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PREFACE 


This book is intended primarily for the student of applied 
mathematics, physics, chemistry or engineering who wishes 
to use the ** special ” functions associated with the names of 
Legendre, Bessel, Hermite and Laguerre, It aims at provid¬ 
ing in a compact form most of the properties of these 
functions which arise most frequently in applications, and 
at establishing these properties in the simplest possible way. 
For that reason the methods it employs should be intelli¬ 
gible to anyone who has completed a first course in calculus 
and has a slight acquaintance with the theory of differential 
equations. Use is made of the theory of functions of a 
complex variable only very sparingly, and most of the book 
should be accessible to a reader who has no knowledge of 
this theory. Throughout the text an attempt is made to 
show how these functions may be used in the discussion 
of problems in classical physics and in quantum theory, 
A brief account is given in an appendix of the main 
properties of the Dirac delta “ function ”, 

I should like to record my debt of gratitude to the late 
Sir John Lennard-Jones, and to my colleagues Mr. B. Noble 
and Dr. J. G. Glume for their generous help in reading the 
first draft of the manuscript and making valuable sugges¬ 
tions for its improvement. I am indebted to Miss Janet 
Burchnall for her assistance in the preparation of the final 
manuscript, to Mr. J. 3. Lowndes for help in correcting 
proof sheets and to Miss Elizabeth Gildart for preparing 
the index. I should also like to thank Dr. D, E, Rutherford, 
general editor of the series, for his advice and criticism 
throughout the preparation of the book. 
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PREFACE 


My debt of gratitude to Professor T. M, MacRobert is 
a much more general one. It was at his lectures that I first 
acquired a taste for the subject* and it will be obvious to 
anyone who knows his published writings how much I have 
been influenced by them, 

Kekle, Staffordshire 
2(M August 1955 
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CHAPTER I 


INTRODUCTION 

§ 1* The origin of special functions. The special functions 
of mathematical physics arise in the solution of partial 
differential equations governing the behaviour of certain 
physical quantities. Probably the most frequently occur¬ 
ring equation of this type in all physics is Laplace’s equation 

VV = 0 (LI) 

satisfied by a certain function ij/ describing the physical 
situation under discussion. The mathematical problem 
consists of finding those functions which satisfy equation 
( 1 . 1 ) and also satisfy certain prescribed conditions on the 
surfaces bounding the region being considered. For 
example, if $ denotes the electrostatic potential of a system, 
i// will be constant over any conducting surface. The shape 
of these boundaries often makes it desirable to work in 
curvilinear coordinates q ly q 2y q 3 instead of in rectangular 
Cartesian coordinates y, z, In this case we have relations 

X = x(q u q 2 , q 3 ), ^ q 2 , q 3 ), z = z(q u q 2l q 3 ) (1.2) 

expressing the Cartesian coordinates in terms of the curvi¬ 
linear coordinates. If equations ( 1 , 2 ) are such that 

dx dx dy + dz dz _ ^ 

dq t dqj dq^qj dq t dq } 

when i 7 ^ j we say that the coordinates q u q 2y are 

orthogonal curvilinear coordinates.! The element of 

t D. E. Rutherford, Vector Methods (Oliver and Boyd, 1939), 
pp. 59-63. 

I 















(1.3) 

(1.4) 


(1.5) 
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length dl is then given by 

dl 2 = h\dq\+h\dq\ + h\dql 

vi/liiprp 

and it can easily be shown that 

= _J_ jjLA*** + _d_ 

h ih 2 h 3 \Sq l \ hi dqj dq 2 \h 2 dq 2 ) 

+ ±.(hhiW\\ 
dq 3 \ h 3 dqjy 

One method of solving Laplace’s equation consists of 
finding solutions of the type 

'A = 61(91)62(92)63(93) 

by substituting from (1.5) into (1.1). We then find that 

LJL( h J t * d Qi\+ 1 a (ihjhSQ^ 

Q i dq 1 \h t dqj Q 2 3q 2 \ h 2 dqj 

+ X±(hh°Qi) = o. 

Q 3 dq 3 \h 3 dq 3 ) 

If, further, it so happens that 

=/i( 9 i)^( 92 . 43) 


etc., then this last equation reduces to the form 
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Now, in certain circumstances, it is possible to find three 
functions g ^ft), g 2 (q 2 ), g 3 (q 3 ) with the property that 

^"1(921 93 ) 0 l( 4 l) + ^ 2 ( 93 > 9 1)02(92) + ^3(9 L> 92)93(43) = 0 - 
When this is so, it follows immediately that the solution 
of Laplace’s equation (1.1) reduces to the solution of three 
self-adjoint ordinary linear differential equations 

It is the study of differential equations of this kind which 
leads to the special functions of mathematical physics. 
The adjective “ special ” is used in this connection because 
here we are not, as in analysis, concerned with the general 
properties of functions, but only with the properties of 
functions which arise in the solution of special problems. 

To take a particular case, consider the cylindrical polar 
coordinates ( 0 , <f>, z) defined by the equations 


x = q cos <j), y = q sin <j>, z = z 

for which h t = 1, h 2 = 5 , h 3 = 1. From equation (1.5) 
we see that, for these coordinates, Laplace’s equation is 
of the form 


d 2 <A l_ _1_ 3V , 3 V 

dg 1 q dg g 2 d<f> 2 dz 2 


(1.7) 


If we now make the substitution 


lA = fi(e)4»(0)Z(z), (1.8) 

we find that equation (1.7) may be written in the form 
I fd 2 R 1 dR\ 1 d 2 0 1 d 2 Z _ 

R\dg 2 g dg) g 2 <S> d <}> 2 X dz 1 

This shows that if d>, Z, R satisfy the equations 


(1.9a) 
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= 0 , 


(1.9b) 


d 2 Ri , 1 dR , 



(1.9c) 


respectively, then the function (1*8) is a solution of Laplace’s 
equation (1.7). The study of these ordinary differential 
equations will lead us to the special functions appropriate 
to this coordinate system* For instance, equation (1,9a) 
may he taken as the equation defining the circular functions. 
In this context sin(fl$) is defined as that solution of 
(1,9a) which has value 0 when <f> = 0 and cos (n<p) as that 
which has value 1 when 0 = 0 and the properties of the 
functions derived therefrom, cf. ex, 4 below. Similarly, 
equation (1.9b) defines the exponential functions. In actual 
practice we do not proceed in this way merely because 
we have already encountered these functions in another 
context and from their familiar properties studied their 
relation to equations (L9a) and (L9b), The situation with 
respect to equation (L9c) is different; we cannot express 
its solution in terms of the elemen tary functions of analysis, 
as we were able to do with the other two equations. In 
this case we define new functions in terms of the solutions 
of tills equation and by investigating the series solutions 
of the equations derive the properties of the functions so 
defined. Equation (1,9c) is called Bessel’s equation and 
solutions of it are called Bessel functions* Bessel functions 
are of great importance in theoretical physics; they are 
discussed in Chapter IV below* 


§2. Ordinary points of a linear differential equation* 

We shall have occasion to discuss ordinary linear differen¬ 
tial equations of the second order with variable coefficients 
whose solutions cannot be obtained In terms of the elemen¬ 
tary functions of mathematical analysis. In such cases 
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one of the standard procedures is to derive a pair of linearly 
independent solutions in the form of infinite series and from 
these series to compute tables of standard solutions. With 
the aid of such tables the solution appropriate to any 
given initial conditions may then be readily found. The 
object of this note is to outline briefly the procedure to 
be followed in these instances; for proofs of the theorems 
quoted the reader is referred to the standard textbooks*! 

A function is called analytic at a point if it is possible 
to expand it in a Taylor series valid in some neighbourhood 
of the point. This is equivalent to saying that the function 
is single-valued and possesses derivatives of all orders at 
the point in question* In the equations we shall consider 
the coefficients will be analytic functions of the independent 
variable except possibly at certain isolated points. 

An ordinary point x - a of the second order differential 
equation 

y"+ce(x)j/+/l(x)y = 0 (2.1) 

is one at which the coefficients oc, /? are analytic functions. 
It can be shown that at any ordinary point every solution of 
the equation is analytic. Furthermore if the Taylor expan¬ 
sions of ol(x) and j}{x) are valid in the range | x—a | < R 
the Taylor expansion of the solution is valid for the same 
range , As a consequence, if a(x) and fi(x) are polynomials 
in x the series solution of (2.1) is valid for all values of x. 

When, as is usually the case, <x(x) and p(x) are poly¬ 
nomials of low degree, the solution is most easily found 
by assuming a power series of the form 

y- £ c r (jc— a)' (2.2) 

r = 0 

for the solution and determining the coefficients c 0 , c u r 2 , 

+ See, for example, E. L, Inee, Ordinary Differentia! Equations* 
(Longmans, 1927), Chap. VII; E. Goursat, A Course m Mathematical 
Analysis, VoL II* Pan II (Girm, 1904), Chap. Ill; j. C* Bufkill, 
Theory of Ordinary Differential Equations Chap. IV (Oliver and Boyd; 
1956). 








6 SPECIAL FUNCTIONS OF PHYSICS AND CHEMISTRY $3 

by direct substitution of (2.2) into (2.1) and equating 
coefficients of successive powers of x to zero. 

The simplest equation of this type is 

y"+y = 0. (2.3) 

Substituting a solution of the type (2.2) with a = 0 into this 
equation we find that, if the equation is to be satisfied, 

£ r(r-l)cX“ ! + £ c t x' = 0. 

r ■ 0 r = 0 

The series on the left is equivalent to 

£ (r+lXr+2)c r+2 * r 

t = o 

so tiiat, equating coefficients of x r , we see that the equation 
is satisfied by a solution of type (2.2) provided that the 
coefficients are connected by the relation 

(r+l)(r+2)c r+2 + c r = 0. (2.4) 

The coefficients c 0 , c 1 are determined by the prescribed 
values of y, y‘ at x = 0, and the others arc determined 
by equation (2.4). From this relation it follows that the 
solution is 

' - c ° ( I_ ir ■ + x i -•)| +Ci (*- ft + r, ~ •••)■ (2 ' 5) 

An equation of the kind (2.4) which determines the 
subsequent coefficients in terms of the first two is called a 
recurrence relation. 

§ 3. Regular singular points. If either of the functions 
a(x), is not analytic at the point x — a> we say 
that this point is a singular point of the differential equation. 
When the functions ct(x), fi(x) are of such a nature that 
the differential equation may be written in the form 

(x - a) 2 y"+(x - a)p(x)y'+q(x)y --- 0 


(3.1) 
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where p(x) and q(x) are analytic at the point x = a we 
say that this point is a regular singular point of the differ¬ 
ential equation. 

If x = a is a regular singular point of the equation (3.1) 
it can be shown that there exists at least one solution of 
the form 

£ cAx-a)**' (3.2) 

r = 0 

which is valid in some neighbourhood of x = a. More 
specifically, if the Taylor expansions for p(x), q(x) are 
valid for |je-o|<2?, the solution (3.2) is valid in the 
same range. 

Putting 

p(x) = £ p r (x-aY, q(x)= £ q t (x~a) T (3-3) 

r = 0 t = 0 

and substituting the expansions (3.2) and (3.3) into (3.1) 
we see that for the equation (3.1) to be satisfied we must 
have 

£ c r (p+ r )(p +r - l)(x - a) e+f 

t = 0 

+ £ PJix-ay £ c ,(e + r)(x-flf + r 

m ° 0 r » 0 


+ I <h(x-ay £ cAx-ay+' = 0. (3.4) 

j “ 0 r = 0 

Equating to zero the coefficient of (*—a)# we have the 
relation 

c 0 q(q ~ i) + Po c o P + <Iq € q = 0 
so that if # 0 we have the quadratic equation 

a 2 -KPo-l)g + So ^0 (3.5) 

for the determination of q. This is known as the indicia! 
equation. Similarly if we equate to zero the coefficient of 
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{x-a) e+r we obtain the relation 

r 

c£e + r )(o + r -1) + £ {p£e+r - s) + q s )c r ., = 0 

S =» 0 

which may be written in the form 

c r {(Q+r)(e + r— I)+i? 0 (£ + r)+ q 0 } 

+ £ {Ps(Q±r-s) + q !C }c F _ x = Q. (3,6) 

s = 1 

Equation (3.5) gives the two possible values q 2 of q. 
If we take one of these values, Q t say, and substitute it 
in the recurrence relation (3.6) we obtain the corresponding 
value of the coefficients c r and hence the solution 

J'iW - £ cXx-a)'**'. 

r=0 

In a similar way the root g 2 of the indicia! equation leads 
to the solution 

y 2 (x) = £ c r (x-ay + <\ 

T w 0 

Three distinct cases arise according to the nature of the 
roots of the indicial equation. 

Case (i) Qi™Q 2 neither zero nor an integer 

In these circumstances the solutions ^i(x) and y 2 (x) 
are linearly independent and the general solution of 
equation (3.1) is of the form 

y= i c£x~ay + “+ £ c'£x- a y + ‘\ (3.7) 

r 0 r » 0 

Case (ii) q x = q 2 . 

If q 1 — q 2 the solutions ^(je) and y 2 (x) arc identical 
(except, possibly, for a multiplicative constant). The 
general solution of the equation can be shown to be 
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yi(x)+y 2 (x) where 
yi(x) = (x-df t 2 c '( x ~ a Y> 

r = 0 

yi«-.Ki<*Mog(x-a)+(x-<>r £ (^r) (*-# 

t - O \ uQ/q = 


Case (iii) q 2 — Qi~n where n is a positive integer . 

In this case all the coefficients in one of the solutions 
from some point onwards are either infinite or indeter¬ 
minate. It can be shown that the appropriate solutions are 


£ c,(x-ay, 

r o o 

00 


+ X b£x-a) r , 

T - 0 j 

where g„ is the coefficient of x* in the expansion of 


(3.9) 


B+I 


{y t (x + n)} : 


exp 


J up(«)dtij. 


It may happen that g„ = 0 in which case y 2 (x) does not 
contain a logarithmic term. 


§ 4. The point at infinity. In many problems we wish 
to find solutions of differential equations of the type (2.1) 
which are valid for large values of x. We seek solutions in 

the form of infinite series with variable -• If we make 

x 

the transformation 


1 



the “ point at infinity ” is taken into the origin on the 
£-axis. With this change of variable equation (2.1) becomes 


d 2 y 

de 


+ 




^ + -4/^V=o. 
« z* \d 


(4.1) 
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rV(r‘) are both 0(1) as £-0 and 
analytic in < then £ = 0 is an ordinary point of equation 
(4.1) and we say that x = co is an ordinary point of 
equation (2.1). Returning to the original independent 
variable we see that the conditions foi the point at infinity 
to be an ordinary point of equation (2.1) are that 

«(*) = - + 0(.v -2 ), fi(x) = 0(x -4 ) as x— qo. (4.2) 
x 

The corresponding solutions are of the form 

r - 0 

Similarly if, as x—co 

*(*) = +0(x~ 1 ), p(x) = & +0 (x- 3 ) (4.3) 

X Xj 

where ct 0 , fi 0 are constants we say that the point at infinity 
is a regular singular point of the equation (2.1). The 
corresponding indicia! equation is 

e 2 +(i-a)e+(} = o. 

If the roots of litis equation are q u q 2 the solutions of (2,1) 
valid for large values of x are of the form 

y t (x) = £ c r x-‘*-', y 2 (x)= £ (4.4) 

r " 0 r - 0 

§ 5, The ganima function and related functions. In 
developing series solutions of differential equations and 
in other formal calculations it is often convenient to make 
use of properties of gamma and beta functions. The 
integral 
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converges if n >0 and defines the gamma function. Similarly 
if m>0, n >0 the beta function is defined by the equation 

B(m, n) = j 1 x* -1 (l-je)" _1 dx. (5.2) 

It is then easily shown that t 

w no = i, 

(ii) T(n + 1) = nT(n), 

(iii) r(n +1) = n ! if n is a positive integer, 

(iv) B(m, n) = 2 J sin 2 " 1-1 0 cos 2 " -1 0 dO, 

(*)«■».») = 2 ^. 

r(m + n) 

(vi) T(i) = Jn, 

(vii) T(p)T(I ~p) = n cosec (pn), 0<p<l, 

(viii) r(i)T(2n) = 2 2 '-T(«)r(« + i)—the duplication 
formula, 

(ix) T(z + 1) = lim -—-, (z>0). 

V .-to (z+ l)(z + 2)...(z + n) 

When n is a negative fraction T(«) is defined by means 
of equation (ii); for example 


rx-fl-gri)- m? 

-1 (-1X-1) 


41X1) 


By means of the result (ix) we can derive an interesting 
expression for Euler’s constant, y t which is defined by the 
equation 

y- lim ( 1 + - -log n) =0.5772... (5.3) 

oo \ 2 n / 

+ For proofs of these results the reader is referred to R. P. Gillespie, 
Integration (Oliver and Boydl, 1951, pp, 90-95. 
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From (ix) we have 

^-{logr(z+l)}= iim flogn---—-—^ 

\ z+1 z-p2 z+n/ 

so that letting z-^0 we obtain the result 


and from (5.1) we find 

r 


log f(z + l) 




' log ( iIt 


-J> 

Integrating by parts we see that 

J e _, log (dt = +e~ x log z + 


(5.4) 


(5.5) 


f ?* 


so that 


- 7 -lim(| ^-dr+logz\ (5.6) 

Closely related to the gamma function are the 
exponential-integral ei(x) defined by the equation 


ei(.x) = — du (jc>0), 

Jx « 

and the logarithmic-integral li(jc) defined by 
«(*)- 

Jo log a 


(5.7) 


(5.8) 


which are themselves connected by the relation 

e i( x ) = (5,9) 

Other integrals of importance are the sine and cosine 
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integrals Ci(jc), Si(.v), which are defined by the equations 

Ci(x) - - P” — du, Si(x) - P — du (5,10) 

J* « Jo U 

and whose variation with x is shown in Fig, 1. 

In heat conduction problems solutions can often be 
expressed in terms of the error-function 

erf (x) = -^-\ e~ ul du (5.11) 

x/njo 



whose variation with x h exhibited graphically in Fig, 2.f 
Similarly in problems of wave motion the Fresnel 
integrals 

C(x) = I cos 2 )du, S(x) = r sin Qnu 2 )dtt (5,12) 
Jo Jo 

occur. The variation of these functions with x is shown 
in Fig. 3, 

t A. C. Aitken, Statistical Mathematics fOliver and Boyd, Seventh 
Edition, 1952), p. 62, gives a short table of values of erf (a). 










erf(x)-»- 




Fig. 3. Variation of the Fresnel integrals, C(x) and S(x\ with x. 
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The importance of these functions lies in the fact that 
it is often possible to express solutions of physical problems 
in terms of them. The corresponding numerical values can 
then be obtained from works such as E. Jahnke and 
F. Emde, “ Funktionentafeln ” (Teubner, Leipzig, 1933) in 
which l hey are tabulated. 


Examples l 

(1) Show that, in spherical polar coordinates r, 6 f 4> 
defined by 

x — r sin 0 cos y = r sin 0 sin z — r cos 0, 


Laplace’s equation becomes 


dr \ dr) sin 0 d0\ 80/ sin 2 


8 V 


8 dip 2 


= 0 , 


and prove that it possesses solutions of the form r a e tm *& 
(cos 0), where Q(ft) satisfies the ordinary differential equation 

O-V)£5 tU+D- J^je-o. 

dr dp l 1 -/i 2 J 

(2) Show that if 

x — a cosh £ cos y — a sinh £ sin z = z 

Laplace's equation assumes the form 

d 2 ijf d 2 tff 2 ( u2 - 2 a 

-rzz (cosh 2 c-cos 2 *;) —^ = 0. 

&T %dz 

Deduce that it has solutions of the form 
where Jlq) satisfies the equation 

4 +(C+l&f cos 2t])f = 0 
dtj 2 

in which G is a constant of separation and q = — a 2 y 2 /32. 
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(3) Parabolic coordinates <f> arc defined by 
* = \A?'/> cos tf>, y = V(^|/) sin <f>,z = 

Show that in these coordinates Laplace’s equation becomes 


_4_ a 

s+vk 




4 a 


i+i? ^ \ 8ii) d4> 2 

Prove that if F„(x) is a solution of the equation 

d 2 F , dF , / m 2 \ _ n 

dx 2 dx V 4 x) 

then F„(£)F_ n (q)e ±lm * is a solution of Laplace’s equation. 

(4) Defining cos a, sin x to be the solutions of 

3 + '- 

which respectively are 1,0 when x = 0, prove 


(i) cos (-a) = cos .v, sin (-a) = -sin x; 

(ii) cos( a+ a') = cos xcos jc'-sin a- sin x 

(iii) sin (*+*') = sin x cos x' + cos x sin a' ; 

(iv) cos 3 A+sin z A = 1; 

(v) — (sin a) = cos a, ~ (cos a) - - sin x. 

dx dx 


(5) The only singularities of the differential equation 
y”+p(x)y’+q(x)y = 0 

are regular singularities at a = I of exponents a, a' and 
at a = - 1 of exponents /}, ft' the point at infinity being 
an ordinary point. Prove that = -a, = -a and 

that the differential equation is 

(A J -l)V + 2(x-lX*-a-a')/+4«a'y =0. 
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Show that the solution is 




where c, and c 2 are constants. 

(6) Apply the method of solution in series to the 
equation 


dx 1 


dx 


showing that, near a = 0, y = Au+Bv where « is a 
Maclaurin series and a = a 1 " °e*. (a is not an integer.) 

(7) Find two solutions of the equation 

(x 2 + 2x)^ + ^-k(k+l)y = 0 
dx 2 dx 


in the form 


y= £ 

It = 0 


Show that* if A is a positive integer, one of these solutions 
is the polynomial 


n~i (fc-n + l)! (2n)l 
(8) Prove that if s is an integer and a is fractional 

,,r(o)r(i-a) 


r(a-s) = (-l)‘ 
(9) Show that 

(i) (a-lX«)„-, 

(ii) («)„_,«(~Xf (a) " 


ni-a+s) 
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(iii) 




(n-s)! 

where 

(«)r = «{«+ l)(a + 2)...fa + r- 1). 

(10) Prove that 

V 2 v 3 

ei(x)= -y-logx+jc - + -. 


2.2! 3.3! 


and deduce that 


Ci(x) - y+log X- + —j 


CHAPTER U 


HYPERGEOMETRIC FUNCTIONS 


§ 6. The hypergeometric series. The series 

\+^l x + ^±iM±i) x 2 + 

l.y 1 * 2y(y +1) 


(64) 


is of great importance in mathematics* Since it is an obvious 
generalisation of the geometric series 

l-tx+x 2 + 

it is called the hypcrgeometric series* It is readily shown 
that, provided y is not zero or a negative integer, the 
series is absolutely convergent if [a|>I, divergent if 
| a- |> 1 , while if | a | = I the series converges absolutely 
if y>E*f$.t I* is convergent when x = —I, provided 
that y>Q£+jj— 1. 

If we introduce the notation 


(a) r = «(« + l)...(a + r-l) = B^ (6.2) 

r(«) 

we may write the series (6.1) in the form 


iF t ( a, p; y, *) = 


y («m A , 
r = 0 r !(y), ’ 


(6.3) 


the suffixes 2 and 1 denoting that there are two parameters 
of the type a and one of the type y* We shall generalise 


+ See J. M. Hyslop, Infinite Series t Fifth Edition (Oliver and Boyd, 
1954), p. 50. 
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this concept at a later stage (§ 12 below) but it is advisable 
at this stage to denote the “ ordinary ” hypcrgeometric 
function by the symbol 2 F t instead of simply F, if we are 
to avoid confusion later. From the definition (6.3) it is 
obvious that 


*; r> *) =!?,(«, /?; r, *)■ (6.4) 


A significant property of the hypergeometric series 
follows immediately from the definition (6.3). We have 


dx 


2 F t (a, 0 ;y; x) = 


r -i (r—l)!(y). 


_ y 1 1 r 

r!(y) r+ , * 

Now (a) r+l = a(a+l) r so the right-hand side of the last 
equation becomes 


afi y (oe+l) r (^-t-l) f 

y r = 0 r!(y + l) r 


showing that 

~ 2 F,(a, 0; y\ x) = ^ iF 1 (a+1, 0+1; y + 1; x). (6.5) 

It should also be observed that 


so that 


2 F,(a, fi;y; 0)= 1 


— P; y;x) 



( 6 . 6 ) 

(6.7) 


Several well-known elementary functions can be 
expressed as hypergeometric series; examples of them are 
given in ex. 1 below. 

It should be noted that, if we adopt a certain convention, 
a hypergcometric series can stop and start again after a 
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$7 


number of zero terms. For example, consider the hyper¬ 
geometric series 2 F,(-«; b ; — n-m; x) where both m 
and n are positive integers and b is neither zero nor a 
negative integer. Because of the occurrence of (-») r in 
the numerator in the expansion in powers of x it is obvious 
that the (n + I)tli term of the expansion will be zero, and 
we are tempted to think that every subsequent term is 
also zero. If we note that, as a result of ex. 9 (iii) of 
Chapter I, 


ni 


(-n), _ 

(-fl-tn), (n + m)! 


(n + m-r)(n + m-r— l)...(it-r+1) 

( 6 . 8 ) 

when the form on the left is not of type 0/0, and if, further, 
we assume that it still has the value on the left when it 
is indeterminate, we see that we may write 

iF,(—n, fc; -ti-m; x) 


- £ fi—=-Yi—r_A.A—^(^(6.9) 

r = o\ n+mj\ n + m-lj \ n+\J rl 


so that although the series stops at the «th term it starts 
up again at the (n+m + l)th term. For instance. 


2F.<-2, l; -5; x) 

= 1+ l x+ 1 x i__L x 6_2 x 7 -jc8 + 

5 10 10 5 

According to a different convention, however, the 
hypcrgeometric function does not restart after a set of 
zero terms. 


§ 7. An integral formula for the hypergcometric series. 
In order to derive some further properties of the hyper- 
geometric series we shall first of all establish an expression 
for the series in the form of an integral. It is readily shown 
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that 


(Ph = W + r, y-fl _I_ p _ ,_ t p+r „ 

(y)r B(fi,y-P) B(fi,y-P) J 0 ( 0 " 


from wliich it follows that 
2 F,(a, /?; y; .v) 


= —-— f ^x r I' (i-tY~ f ~'t r * r ~ l dt 

Interchanging the order in which the operations of sum¬ 
mation and integration are performed we see that 

2 F t («, P‘, y; .y) 

Using the fact that 


1 ( 4 (* 0 ' = (i-y 0 -*, 

t = o r I 

we have the integral formula 
2 Fi(of s ft; y; x) 

valid if | x | < J, y>/?>0. The results hold if x is complex 
provided that we choose the branch of (I in such a 

way that (1 as /-*0 and 0. 

The first application of (7.1) is the derivation of the value 
of the hypergeometric series with unit argument. Putting 
x = 1 in (7.1) we have 

,fm,P\t, l? = - 1 ^ f'(i 

B(p,y-p)j o 
_ m y—a—P) 

B(P, y-p) 
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if y—a—/j> 0,/?>()* Ifwe express the beta function in terms 
of gamma functions we have Gauss’s Theorem 


z F l (<x, p; y; 1) = 


r(y)r (y-a-Pi 

T(y-a)r(y-p) 


Now if of = —n, a negative integer, we have 


(7.2) 


r(y-«-/?) 
r(y-« 


(y-flu 


r(y-«) 

Hr) 


“( 7 ). 


so that equation (7.2) reduces to 


2 F l {-n t p-,r, 1)= ( -^= 

(v)» 

which is known, in elementary mathematics, as Vander¬ 
monde's theorem* 

Again, if we put a- = — 1 and a = 1 + /?—y we have, 
from equation (7J) 


2 Fi(M; ft —ct+1; -1) = 


m+p-ct) f 1 

rw(i-«)J 0 




If we write £ = t 1 in this integral we sec that its value is 
iB($p, I - a). Using this result and the relation iT(£/J)/n/J) 
= f(l +iP)IW +p) we have Kummer’s theorem 


2 Fi(«, p, p-ot+l ; -1) = 


ni+p-«)ni+ffl 

m+»m+ij 8 -«)‘ 


(7.3) 


Further, we can deduce from the formula (7.1) relations 
between hypergeometric series of argument x and those 
of argument xj(x—\). Putting t = 1— t in equation (7.1), 
and noting that 

{l-x(l~T)}- a = (l-x)-‘|l- 
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wc see that 


2 ^i( a > Pi y; *) 

my-P)h l x-i J 


2 fi(«, Pi y;x) = (1 y—pi y; (7.4) 

and* by symmetry, the relation 

Pi y; x) = ( 1 -a) _ * 2 F 1 ^y-a, /?; 7 ; ( 7 *5) 

Using the symmetry relation (6.4) and equation (7.4) with 
a- replaced by a/(a— I) we see that 

= (1-X y~ r zFi{y~p, y-a; y; x), 

so that 


2 ^, P; y; x) = (1 -a)’-*" V^y-a, y-p; y; a). (7.6) 
If we put x = i in equation (7.4) we obtain the relation 

2 ^i(«, Pi y; = 2* 2 F t (a, y ~P; y; -1). 

The series on the right-hand side of this equation can be 
derived from equation (7.3) provided either that 

y = y~P~ a + 1, i.e. p = 1— a, 

or that 

y = a-(y-/I)+l, i.e. y = $(« + /?+!). 


§8 HYPERGEOMETRIC FUNCTIONS 

We then obtain the formulae 

r(iy)r(*y + i) 


^(a, 1-a; y; i); 


r(i«+iy)r(i-i«+iy) 


rr- ft- t r , iff, I.n- mm+fr+W 

lF,< ’’ "■ i “ +w+i ' J) - r ( 1 + J«)r(i+iffi- 
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(7.7) 

(7.8) 


§ 8. The hypergeometric equation. In certain problems 
it is possible to reduce the solution to that of solving the 
second order linear differential equation 

x(l-x)p^ +{y-(l+a+P)x} d ^ - a p y = 0 (8.1) 

dx* dx 

in which a, fi and y are constants. For instance, the 
Schrodinger equation for a symmetrical-top molecule, 
which is of importance in the theory of molecular spectra,! 
can, by simple transformations, be reduced to this type. An 
equation of this type also arises in the study of the flow of 
compressible fluids. In addition certain other differential 
equations (such as that occurring in ex, 1 of Chapter I) 
which arise in the solution of boundary value problems in 
mathematical physics can, by a simple change of variable, 
be transformed to an equation of type (8.1). Indeed it can 
be shown that any ordinary linear differential equation 
of the second order whose only singular points are regular 
singular points, one of which may be the point at infinity, 
can be transformed to the form (8.1). For that reason it 
is desirable to investigate the nature of the solutions of 
this equation, which is called the hypergeometdc equation. 

We may write the hypergcometric equation in the form 

x 2 y"+x(l+x+x 2 +...){y~(a+P+l)x}y' 

-apx(l+x+x 2 +...)y = 0, 

t See, for example, L, Pauling and E. B. Wilson, Introduction to 
Quantum Mechanics> with Applications to Chemistry (McGraw-Hill, 
New York, 1£35), pp. 275-280, and ex. 10 below, 
c 
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so that, in the notation of § 3, we see that near x = 0 

Po = y, <?o = o, 

and the indicial equation is 

0 2 +(y-l )0 = O 

with roots g = 0 and q = 1 -y. 

Similarly, the equation can be put in the form 

(* -l)V—(x-l){y-a—jS-l—y(x—l)+...}y' 

+ a/J(x-l){l~(x—l) + ...}y = 0, 

with indicial equation 

£ 2 -Kfl£ + /?-y)p = 0, 

of which the roots are q = 0, q - y-a-f}. 

Finally in the notation of § 4 we have for large values 
of x 

t \ + 1) v Of 8 

X x 2 


and so the indicial equation appropriate to the point at 
infinity is 

q 2 *- (at+/?)£ +(*/? = {), 

with roots a, f} t 

Thus the regular singular points of the hypergeometric 
equation are:— 

(i) x — 0 with exponents 0, 1 -y; 

(ii) x = oo with exponents a, /?; 

(iii) x — I with exponents 0, y—ot — fi< 

These facts are exhibited symbolically by denoting the 
most general solution of the hypergeometric equation by 
a scheme of the form 


oo 

a 

P 


0 

y = P\ 0 

U-y 


l 

o 

y-a-p 


X 


( 8 . 2 ) 
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The symbol on the right is called the Riemann-J’-function 
of the equation. 

We shall now consider the form of the solutions in the 
neighbourhood of the regular singular points. 

(a) x = 0: Corresponding to the root q — 0 we have 
a solution of the form 

y= I <vc r . 

w = 0 

Substituting this series into equation (8.1) we obtain the 
relation 

(l-x) t vfr-iy 1 

r *= 0 

+ {y-(«+jS + l)x} £ c r r* r_1 -a0 £ c,x r = 0 

r'tO F ** 0 

which is readily seen to be equivalent to 


Z {c r+ i[r(r + l) + (r+l)y]-c r (r+a)(r+^)}x r = 0, 

r — 0 

so that 


(r + oc)(r+0) „ 

Crtl (r+l)(r+rt " 


(8.3) 


from which it follows that 


c r = 


(y) f r! 


c o* 


(8.4) 


It follows that the solution which reduces to unity when 
x = 0 is 


y = i+ y, x+ 

yl! 


a(a+1)0(0+1) , 

V(J + 1)2! 


i.e. y = zFi(a, Pi y; x). (8.5) 

Similarly, if 1 — y is not zero nor a positive nor negative 
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integer, the solution corresponding to the root q = 1 —y is 

y = 

where 

(l-x) £ C,(r+l-yXr-7K‘ 7 + 

r = 0 

+ {y-(a+j3+l)*} £ c t (r+l~y)x'-’-ap £ c^ l " T+r = 0, 

r ~ 0 r-0 

which is equivalent to 
£ c T {(r +1 - y)(r - y)+y(r+1 - y)}x r_ 

r e 0 


- £ c r {(T+l-7)(r-7)+(«+^+i)(r+l-)»)+ fi(J 5}x r -^ 1 = 0, 

r “ 0 


implying that 

C r+l 


(r+«-y + l)(r + /?-y+l) 
0 +l)(r+ 2 -y) 


Comparing this relation with (8.3) and taking c 0 = 1 
we see that this solution is 


jci-r^a-y+l, 0 -y+l; 2 -y; x>. ( 8 . 6 ) 

Combining equations (8.5) and ( 8 . 6 ) we see that the 
general solution valid in the neighbourhood of the origin is 

y = A 2 F t (a, jS; y; x) 

+Bx 1 -> 2 F 1 (a-y+1, p-y+l; 2-y; x ), (8.7) 

provided that I — y is not zero or a positive integer. 

If y = 1, the solutions (8.5) and ( 8 . 6 ) are identical. If 
we write 

7iW = P‘, y; x ) 


ya(*) “ J'lW log x -F 


to 


Z vf 

t - I 


and put 
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|8 -- 

we find on substituting in (8J), with Y — 1, that 

(r+ l) 1 c, + 1 —r(a + ^+ l)c,+ =0 

r!(r+l)I 

from which the coefficients c r may be determined. 

A similar procedure holds when I — y is a positive integer, 

(b) x = 1 : If we let £ = I - x, equation (8.1) reduces to 

«i-8 S +{“+/*-)■+1 -(«+*+D0 % -*Py = o 

at at 

wiiich is identical with equation (8.1) with y replaced by 
a+/J—y+1, and x by % = 1-x. Hence it follows from 
equation (8.7) that the required solution is 

y = A 2 F l (a, p; a+jS-y + 1; 1-x) 

+ B(l-xy-'-' , 2 F l (y-«,y-p-,y-a-p+l; 1-x). ( 8 . 8 ) 

(c) x = oo: Corresponding to the root g = a, we put 


y= £ 


r = 0 


wiiich gives 


(l-x) £ cXr + aXr + a+lJx-'-*-' 

r-0 

-{y-(a + )3+l)x} £ (r+ccfC'X'''*" 1 ~aP £ c,*"'"* = 0, 

r *■ 0 r = 0 

i.e., 

£ cXr+aXr+a-y+l)*-'-*- 1 - £ c r r(r+x-p)x~ r -‘, 

t « 0 r ** 0 

whence it follows that 


- O'+ <*)(*•+«-y + i) 

C,+1 (r+ l)(r + «—/?+!) ' 
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which in turn is equivalent to 


c = («) r (tf-y+ l) r 

' r\(ct~p+1 ), 

Taking r 0 = ! we obtain the solution 


c 0 . 


* *2*1 a-y+ 1 ; u-p+l;~J. 

From the symmetry we see that the other solution is 
x~ e 2 F i fi—y+1; p~a+U~j, 

so that the required solution is 

y ~ Ax~* 2 F t (a, a-y+1; a-p+i ; ^ 

+Bx-“ i F i (p,p-y+lip-a+l;^. 


(8.9) 


§ 9. Linear relations between the solutions of the hyper- 
geometric equation. The series in the solution (8.7) are 
convergent if | x |<I, i.e. in the interval (-1, 1) whereas 
those in the solution (8,8) are convergent in (0, 2). 
There is therefore an interval, namely (0, 1), in which all 
four series converge, and since only two solutions of the 
differential equation are linearly independent it follows 
that there must be a linear relation, valid if 0<x<l, 
between solutions of type (8.7) and those of type (8.8). 

Let 


P’> T, x)=*A 2 F i {a, P: tx+p—y+l; L~x) 

+ B(l~xy~*~ fi 2 Fi(y~a, y-/J; y-ot-/? + l; l~x). 
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then putting x = 0 we have 

1 =A 1 F l (a, p\ <x+P-y+l; 1) 

+ B 2 F i(y-*’ y~P'> y-a-p+l; 1), 

and putting x = 1 we have 

2 F t («> P\ y; i) = 

if we assume that 

l>y>a + p. (9.1) 

Substituting for the series with unit argument from equa¬ 
tion (7.2) we see that 

A r(y)T(y—«—j?) 

r(y- a )T (y-p)’ 

and that 


J = A !"(«+/?—y+l)r(l—y) +B r(y-*-fi+mi-y) 


np-y+l)r(a-y+l) 


r(i-/?)r(i-a) 


so that 


B = 


T(y)r(g+/? — y) 

mm : 


whence we find that 
p\ y; x) 


= fr r( L g P‘’ a+ P- T+ 1 ; 1 -*) 

T(y-<x)T(y-p) 
r(y)T(a+fi—y) .. 

mm 

1 F 1 (y-ot,y-jJ;y-a-0+l; 1-x), (9.2) 
provided that the condition (9.1) is satisfied and 0<*<1. 
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$9 




I f wc replace a - by ^ in equation (9,2) we have 
(a, p ; y; 

i^i 7“0; y-a-jS+1; 1- 0* 


r(y)r(a + /? — y), 
iW(/i) 


and from equation (7.4) 


(«,*,; 1-i) = 


so that 
F 


x* 2 F i(a, y-/J; y; 1-x), 


r(y)r(y-g-P) 

r(y-a)r(y~p) 


x’xFiia, «—y+1; a+p-y+l; 1-x) 


w<» 


iFi(y-a, l-«;y-a~/3+l; 1-x), (9.3) 
where 1<jc<2 and I>y>a+/i, 

These relations are typical of a larger number which 
exist between the solutions of the hypcrgeometric equation 
(8.1). If we change the independent variable in this 
equation to any one of 

1 1 x~l x 
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§ 10 

the equation transforms to one of the same type (but, of 
course, with different parameters). The equation (8.1) 
therefore has twelve solutions of the types (8.5) and (8.6) 
—two for each independent variable—each convergent 
within the unit circle. Any one of these can be expressed 
in terms of two fundamental solutions. In addition twelve 
more solutions of the kinds 

(1 - x )y-^ 1 F 1 (y-a, y~p;y;x), 

x^l-xr^^l-a, I—/?; 2 -y; x) 
can be derived. The relations between these twenty-four 
solutions of the hypcrgeometric equation are of the types 
(9.2) and (9.3); for a full discussion of them the reader is 
referred to T, M. MacRobert, Functions of a Complex 
Variable (Macmillan, 2nd edition), pp. 298-301. 


g 10, Relations of contiguity. Certain simple relations 
exist between hypergcometric functions whose parameters 
differ by ±1. For example if the parameters a and /J 
remain fixed and y is varied we can prove that 

y{y-l-(2y-l-«-jS)A} 2 F,(iJt, /?; y; a) 

+(y-<*)(y-P)xiF 1 (a, p; y + 1; x) 

- 7(7 - 1X1 -x) P\ y — 1 ; 2 () = 0. (10.1) 

The proof follows from the definition (6.3). For the co¬ 
efficient of x? in the expansion of the function on the 
left of (10.1) is 

7(7-1) yfj ~(27 -«- P) (g) "~\ 

(7) n «! (7) n -i(n-l)! 


+ (y-a)(y-P) 


(?+D.- i(n-l)! 


-7(7-1) 


(7-l).«! 


+ 7(7-1) 


CoO-itf)-. 

(y-!).->(«-i)! 


and it is not difficult to show that this is zero. 
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In another kind /J and y are kept constant and a is 
varied. One such is 

{y—u—f}+(f}—a)(l—x)} 2 F l (cc, p; y; x) 

+ «(I-Jc) 2 / r i(c< + l, Pi yi x)-(y~a) 2 F 1 ( a-1, p,y;x)=* 0 

( 10 . 2 ) 

the proof of which is similarly direct. 

In the third type of relation y is kept constant and <x 
and P vary. One of the simplest among these relations is 

(«-/J) 2 F,(a, pi yi x ) = tt 2 F t (a+l, Pi y; x) 

-pzFfa P+lir,x) (10.3) 

The proof of these relations is left to the reader; further 
examples are given below (exs. 3, 4). 


§ 11. The confluent hypergeonietric function. If we 

replace x by x//J in equation ( 8 . 1 ) we see that the hyper- 
geometric function 

2 ^i(«> Pi r, x/p) 

is a solution of the differential equation 




so that letting /?-*oo we see that the function 
lirn 2 F l (a, Pi y; x/p) 

is a solution of the differential equation 

d 2 y ,,*dy ^ 

x~ +(y— a) ~ —ay = 0. 

From the definition of (/?),, we see that 


( 11 . 1 ) 


Jim (Hi = j 

fl-KO p' 


( 11 . 2 ) 
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so that the function (11.1) is the series 
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f (*), 

r -0 (y),' r! 


(11.3) 


and this series we denote by the symbol ,Fj(«; y; x). This 
function is called a confluent hypergeonietric function, and 
the equation ( 11 . 2 ) is the confluent hypergeonietric equation. 

Equations of the type (11.2) occur in mathematical 
physics in the discussion of boundary value problems in 
potential theory, and in the theory of atomic collisions 
(see examples 13, 14 below). 

It is readily verified that the point x = 0 is a regular 
point of the differential equation ( 11 . 2 ) and that, in the 
notation of § 3, p 0 = y and q 0 = 0. The indicial equation 
is therefore 

ete+y-i) = o 


with roots q — 0 and Q = I -y. 

Corresponding to the root q = 0 there is a solution of 
the form 

co 

yi= E c r x r i 

r a 0 

substituting this solution in equation ( 11 - 2 ) and equating 
to zero the coefficient of a*' we find that 


_ (g + r)c r 
,+1 (y+r)(p+l)* 

Putting c 0 = 1 we see that 

c = Wr 1 

' trX'rf 

and if y is neither zero nor a negative integer the solution is 

M*)= iffair, x). (U.4) 

Similarly, the root g = I —y, leads* if is neither zero 
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nor a positive integer, to a solution of the type 

= £ Cr x r . 

r » 0 

If we write 

and substitute in ( 11 . 2 ) we find that w(jc) satisfies the 
equation 

+ ^ 2 ~ 7 ~ x ^Jx -(*-?+ l ) u = 0 ' 

which is the same as equation ( 11 . 2 ) with y replaced by 
2-y and a replaced by et-y + I. We know from equation 
(11.4) that the solution of this equation which has value 
unity when x = 0 is u = ,F,(a-y +1 ; 2 —y; x) so that 

J'aW = * 1 * T i-P|(«-y+l; 2-y; *). (11.5) 

Thus if y is neither 0 nor an integer the general solution 
of equation (31.2) is 

y(x) = A y; x)+Bx l ~ r iF^a-y+1; 2-y; x ), (11.6) 

where A and B are arbitrary constants. 

In the exceptional case y — 1 we have 

yi(x) = iFi(a; 1; x), (11.7) 

obtained simply by putting y — I in equation (11.4). For 
the second solution we write 

y 2 (x) = y t (x) logx+ £ cX- ( 11 . 8 ) 

r = 1 

Substituting this expression in equation (11.2) we find that 
the unknown coefficients c r must be such that 

X i -3'i+ £ {(r + l) 2 c, + t -rc,}x r +Cj =0. 

dx r = l 

Inserting the value of y t (jc) from equation (11,7) we see 
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that these coefficients are determined by the recurrence 
relation 

c lS =l-«,(r+l) 2 c r+ 1 -rc r = (l-«) (11.9) 

r!(r + l)! 

The complete solution is therefore given by y = /ty 1 (-v) + 
By>i{x) where A and B are arbitrary constants and the 
functions y L (x), y 2 (x) are defined by equations (11.7), 
(11.8) and (11.9). The complete solution when y is an 
integer may be found by a similar method. 

If in equation (11.2) we put 

y(x) = x-t’e*’ W(x) (11.10) 

we find that the function W(x) satisfies the differential 
equation 

7? + {- i+ - x + i -if} wix) - 0 ' (1U1) 

where we have written k for iy—a and m for (|—iy). 
The solutions of this equation are known as Whittakers 
confluent hypergeometdc functions. 

If is neither 1 nor an integer the solutions of the 
confluent hypergeometric equation corresponding to equa¬ 
tion (11.11) are given by equation (11.6) with y — 1+2 m 
and a = \—k+nu Thus the solutions of equation (11.11) 
are the Whittaker functions 

M kf m (*) = x*+ m e~* x t FxQ-k+m; l+2m; x% (11.12a) 

Mjt ( - Jx) — x*- m e~* x | —k - m; 1 — 2 m; x). (11,12b) 

Several of the properties of 1 F 1 functions have analogues 
for the yFi functions. Corresponding to equation (7.1) 
there is the integral formula 

tF^jy;*) - —2—- ['(I_,)>-«-»/«-V**, (11.13) 

B(<x, y-«) Jo 
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from which Kummer’s relation 

y; x) = e* y; -x) (11,14) 

may be obtained by a simple change of variable. The 
analogue of equation (6.5) is 

f 71 x)} = - iFi(«+l; y + 1; x), (11.15) 

ax y 

while corresponding to the contiguity relations of § 10 we 
have relations of the type 

K 1 i r 1 (a+ 1 ; y+ 1 ; *)+(?-«) iFi(a; y+ 1 ; x) 

-Yi^i(«; 7 ;*) = o, (11.I6) 

(x+a) 1 F,(a+l; y+ 1 ; x) + (y-ec) ,Fi(a; y+ 1 ; x) 

~y i F 1 (a+l;y;x) = 0, (11.17) 

«iFi(«+l; y; x) + (y~ 2 a-x) ^(a; y; x) 

+(«-y) 1 Fi(a-l;y;x) = 0 , ( 11 . 18 ) 

(a-y)x [F^a; y+ 1 ; x) + y(x+y- 1 ) jF^a; y; x) 

+y(y-l),F 1 (a;y-l;x) = 0. (11.19) 

§ 12. Generalised hypergeo metric series. There are 
two ways by which we may approach the problem of 
generalising the idea of a hypergeometric function. We 
may think of such a function as being the solution of a 
linear differential equation which is an immediate general¬ 
isation of the equation (8,1) or we can define the function 
by a series which is analogous to the series (6,1), 

At first sight it is difficult to see how the differential 
equation (8,1) can be generalised immediately, but if we 
introduce the operator 
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and notice that (8,1) is equivalent to 

{9($+y-l)-x(B+*){9+P))y = 0, (111) 

an obvious generalisation is 

{S(S + 0i — l),..(3 + 0 p —1)~ = 0* 

( 112 ) 


where m u a 2f q ,, g p are constants. Further¬ 

more it is readily shown that this equation is satisfied by 
the series 

y (^ : l)n(^ 2 )fl^ - ^ (12 3 ) 

(Q 1 UQ 2 X •%)« nf 

which is, itself, a generalisation of the series (6.1). Such 
a series is called a generalised hypergeometric series and 
is denoted by the symbol p+ x FJp u ,; q u q p ; x). 
It is left as an exercise to the reader to show that, if no 
two of the numbers 1, q u q 2 , Q P differ by an integer 
(or zero), the other p linearly independent solutions of 
equation (12.2) are 

* 1 ~ e *p+i-f’(i(i+ a i“Ci-"» i+a P+ i-ei; 2-e„ 

1 + 61 -^ 2 . •••> l+Qf-Qi', x), (I = 1 , 2 , ..., n). 

As it stands (12.3) is a generalisation of the series (6.1) 
but it is not sufficiently wide to cover a simple series of 
the type (11.3). To cover such cases we generalise, not the 
differential equation, but the series defining the function. 
The generalisation of (6,1) which includes (12.3) is the series 


y> (« l)> 2 )n---(«p)n ^ 

n = 0 (Ci)fi(02)ti-■•(£',)« a ! 

which we denote by the symbol 


(12.4) 


pFg(CCij ...0f p , "-Qqt X), 

or, if wc wish particularly to throw into relief the difference 
between the numerator and the denominator parameters. 
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by the symbol 

F f a i’ V 

” S Ui. Q'< 

The suffix p in front of the F denotes that there are p 
numerator parameters » lt a p ; similarly the suffix q 

indicates the number of denominator parameters. 

Generalised hypergeometric series do not usually arise 
in mathematical physics because we have to solve equations 
of the type (12.2). Their use is more indirect. Such series 
occur normally only in the evaluation of integrals involving 
special functions. In certain cases these series reduce to 
series of the type 

F h. a ' ; 1 

p xci.e,; _ 

which have unit argument. For this reason it is desirable 
to have information about sums of this type. An account 
of the theory of such sums is given in W. N. Bailey, 
Generalised Hypergeometric Series (Cambridge University 
Press, 1935), Here we shall consider only one such calcula¬ 
tion because it illustrates the use of the theorems of Gauss 
and Kummer proved above (equations (7.2) and (7.3) 
respectively). Other results of this kind are given in 
examples 18 and 20 below. 

By expanding the 3 F 2 series involved we see that 

s _ mmm F r « a n r 

r(i+ a -0)r(i+«-y) 3 2 Li+<*-& i+«~y; J 

_ y r(a+n)r(f} + n)r(y + n) 

» = q n!Rl+a — /J+n)r(l+a—y + n) 
f r(a + n)r(p + n)r(y + n) 

B = o n!r(l+a+2n)r(l+a-^-y) 

r(l+g + 2n)r(l+«-l-y) | 

Rl+a-jS+n)Rl+«-y + /ij 
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Now by Gauss’s theorem (7.2) the expression inside the 
curly brackets is equal to 2 F|(j8+n, y+»; l+a + 2 a; I) 
which may be written 

y R/f+H + ro)Ry + n + ffl)Rl+«+2«) . 
m = o Rj3+n)Ry + n)Rl + a + 2n + m)m! 
whence we find that 

s= y y T(« + b)TQS + b + m)Ry + n + m) 

* = o m- o r(l + a+2« + m)r(l+ce-j?-y)n!m!* 

Interchanging the order of summation and putting p - /i+wi 
we see that 

s= f R/?+p)Ry + p) f Hat + w) 

p- o Rl+<x-/I-y) « = o (t!(p-n)!Rl+a + n + p)' 

Now by example 9 (iii) of Chapter I 

.— ^ «(-!)" (j “. 

(p-w)! p! 

so that the inner sum is equal to 

Ret) F f a, -P‘, -ll 

pIRl+a+p) 1 '[l + a+p J’ 

which by Kummer’s theorem (7.3) is equal to 

Rconi+fr) 

p!Rl + ct)Rl + i« + p) 

Therefore 

„ y Rtt)r(^4-p)Ry + p)Rl + ^t) 
p = o p!(l + a— p— y)Rl + ot)Rl +ia + p) 

. RoQRflRy) F r P, r; i] 

Rl+a)Rl + <x—P—y) 2 J 


D 
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This 2 F, series with unit argument can be summed by 
Gauss’s formula (7.2) and the expression for S found. 

It then follows that 

F f «> P> r, i ] 

3 2 [l+a-p, I+a-y;J 

= rg+i«)r(i+i«-/?-y)r(i+a-jS)r(i+«-y) 

T(1+a)T(l +a-^-y)T(l+ia-^)r(l +i«-y)’ 
a result which is known as Dixon’s theorem. 

Examples II 

(1) Show that 

(0 P\P\z) = 0-z)"*; 

(ii) i«+i;i; z) = i((l-2T)“-+(I+r)-“>; 

(ui) M«+h i* + U \\ z 2 ) = -L {(1 —z) - *—(1+z) - *}; 

2az 

(iv) 2 F 1 (l,l;2;z)=-il0g(l-z); 

z 

(v) zFi(b Ui>z 2 )=~- log 

2 z 1-z 

(vi) 2 Fi(i, i; i; z 2 ) = 

2 

(vii) jF^i, 1 ; - z 2 )= 

(viii) 2 F,(i, i; 1; k 2 )=~ K(k); 

71 

(ix) jF.t-i,*; l;* 2 ) = -£(*). 
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(2) By transforming the equation y"+n 2 v = 0 to 
hypergeometric form by the substitution i; = sin 2 z, prove 
that, if — $ 7 r^z^ jjr, 

(i) cos (nz) = ^(in, -in; i; sin 2 z); 

(ii) sin («z) = n sin z 2 Fi(i—in, i + in; i; sin 2 z); 
and that, if O^zgn, 

(iii) cos (nz) = cos (inn) 2 F,(i«, —in; i; cos 2 z) 

+ n sin (inn) cos (z) jFi(i“in, i + in; i; cos 2 z); 

(iv) sin («z) = sin (i«Jt) 2 Fi({n, — in; i; cos 2 z) 

-« cos (inn) cos (z) 2 Fi(i-in, i + i>»; i; cos 2 z); 

(3) Prove the relations: 

(i) i»-m-x) 2 F l (a,p;r,x) 

= (y-fizFiia, p- 1; y; x)-(y-*) 2 F,(a-l,0; y;x); 

(ii) (y—P~\) 2 F 1 (a t P\ y; x) 

=(y-a-/?-l) 2 Fi(a, j5+l;y; x)+a(l-x) 2 F 1 (a+l, £+1; y; x); 
=(«-^-l)(l- 3 c) 2 F 1 (a,^+l;y;*)+(y-a) 2 F 1 («-l,^+l;y; jc); 

(iii) (y—a —^) 2 F 1 (a, p; y; x) 

= (y-«)jFi(a- 1, P; y; x)-0(l-*)iFi(«, p+ 1 ; y; x); 

(iv) « 2 F,(oe+1; P; y; .v)-(y- l) 2 F,(a, P; y-1; .v) 

= (a+1—y) 2 F|(a, p; y; x); 

(v) (l-x) 2 F 1 (a, /?; y; x)- 2 Fi(a-l, p- 1; y; x) 

= a+ ^ - —- x 2 Fi(a, p; y+1; x); 

7 

(vi) QzP)* 2 F^a, p; y+1; x) 

7 

= 2 F,(a-l, p-l; y; x)- 2 F,(a, /?— 1; y; x); 
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(vii) (l-jc) 2 fi(«, P;y;x) 

= fi-U y; x)+ 2 F 1 (g, jS; y+l; x), 

y 

4. Prove that 

(0 2 ^iC a > P+ 1; y+i; z)- 2 J p 1 (a, /?; y; z) 
y(y+i) 

(ii) 2 F l (« J 0; y; z) = jF^a + l, jS-1; y; z) 

+ g "’^ +1 z 2 Pi(g+l, 0 ; y+l; z). 
y 

Deduce a simple expression for the hypergeometric series 

jFjta, P;P-l,z). 

(5) If n is a positive integer, prove that 

2 F 1 (-n,a+n;y;x) 

= x^^-xr-rty) d" f + m) 

r(y+n) dx" 1 U } 1’ 


and deduce that 


2^1 (-«» «+n; i«+i; 


= (ft a -l)*-*T(^ + i) #_ (2 _ jy+f.-* 
2T(i« + i + n) dp"^ 

(6) If /i is a positive integer, and | x j > 1, prove that 
r /n + 1 « + 2 , l\ {-l)"x" +1 d" f 1 

if ‘ It- t- '• - 

(7) Establish the following formula: 

0) 2 ^ 1 ( 0 *; /?> ^)x 2 Fi( 7i S; y+5 — Q; x ) 

“ 2^1 (»+e»^+e;«+^+e;jf) x 2^1(7 —^jr) - 


W' 
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(ii) D x {x* % F x {a y p; y; kx)} = ax*~ 1 2 Fj(a+1; JJ; y; **), 
where a#0; 

(iii) B(X, y-X) 1 F i (a, 0; y; x) 

- 0;X;xt)dt, 


where | x j<I, 2>0, y—2>0. 
(8) Prove that if 0>O, 

2 F t (a, 0; 20-, z ) 


-.fe ta:: p‘ ( .M)»-r < 1+ f“ s «"'L, 

f) ]„ L+{l-£<»a #>—J 9 

where ( — z/(2-z). 

Deduce that 

iFM 0; 20; z) = (1 -\z)~* ■§«+*; 0+h C a ). 

(9) Prove that 


r 


cos mO cos" Odd = 


jrr(fH-l) 


2 B+1 r(iM + im + 1)r(±n-4m + l) 


and evaluate 


r 


cos 2m0 sin" 0d8 


where m is a positive integer. 

(10) Schrodinger’s equation for the rotation of a sym¬ 
metrical-top molecule is 

sin 0 38 \ dOJ sin 2 8 88 2 

. (cot 2 6+~\^-~ — 005 S ^ ■ % nlAW 

V C) 8y 2 si 


sin 2 8 8yd 4> 


*- 0 . 
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where A, C, W, h are constants. Show that it possesses 
solutions of the form 


tfr = /J; y; x), 


where n^m, x = £(l -cos 0)> 7 = n—m + 1, and «, fi are 
the roots of the equation 


z 2 —(2n + l)z + - n 2 + n- 

c 


%n 2 AW 

it 2 


(II) Prove that: 

(i) ,F,(«; «; x) = e'; 

(ii) t F l (a+ 1; a; x) = ^1 + ^ e*; 

(Hi) iFi(i;|i -x 2 )=^erf(x); 

2x 

(iv) 1 F 1 («+l; 7 ;x)- 1 F 1 (a;y;x)= ~iF,(a+l;y + l;x); 

V 

(v) —x 2 ) = e~ xl ~yjTi x erf (x); 

(vi) *" 1 ^( 11 ; n+1; -x) = « j* t* -1 

(12) Prove that the equation 

d 2 F _ 

5x 2 fc 5f 

possesses solutions of the type 

rn where and C are constants. 
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(13) Show that the Schrddinger equation 
Vty+ fft 2 -= 0 
possesses a solution of the form 


t F i (“ ^ Uikr-ikz'j. 


(14) The Schrodinger equation governing the radial 
wave functions for positive energy states in a Coulomb 
field is 

Show that it possesses a solution 

L = r" e ikT i F x (f# + n + 1; 2?i + 2 ; - 2ikr) 
where k 1 = %n 2 mWjh 2 , a — An 2 mzz'e 2 jk. 

(15) Show that the equation 

2 2 nip 


d2 y . 1 _dy 
dx 2 x dx 




possesses a solution 

y = FMn+i-ip-, 2imx) 

and hence that a solution of equation (1.9c) is 
R = gi , 'e~ it F(in+i; n + 1; 2 img), 

(16) Show that 

(i) Px-a-xr-vJ* 1 - KX "U 

Jo Lpi» J 

-b(1, . 
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(ii) 


f i x , ' 1 0-3c)’"" l pF, 

* 1-x" 

a i> * **> ^ 

Jo 

_0t. P«l 


dx 




(iii) 


[ l u-*r'A 

l-x“ 

Jo 

_Pu ft 


dx 


(17) Prove that 




e ~ax x tt~i J x 


_ _ I~{/0 p f"^li ‘"i ^pj hr ^/^”1 

p+i a U>-.j m J ; 

®JT ^;::Z ±bV y~ x ‘'' ix 

T(ji) „ fa i, ip, in+i\ ±4b 2 /a 2 l 

./>.; J' 




e plxl x p 1 dx 


= M 

2 P" 


p+i 


F r«i.« P . \hi 

q L^i* Pqi 


±b 2 lp 2 


(18) By equating coefficients of x in the relation 
(1 -x)*" fl ~ T 2 F i («, p; y ; x) = jF^y-a, y~Pl T> *) 
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prove Saalschutz’s theorem 

p ["«, p, -n; 1 "1 

3 2 ly, 1+«+jS-v-«J" 


(r-«)„(r-A, 


Hence prove that 


If | x | <3—2^2. 


ia, {+bu-P; 
l+x—p; 


4x ' 

( i-xy 


(19) Show that 


a-Fit#, Pi l+a r p; z) 

= (l-z) _ * 2 Fi(ia, i+i*-P‘> 1 +a —P‘, 0 

where £ = -4z(l — z)~ 2 . Hence deduce the value of 
2 F 1 («, p; I + a— P; — 1) from Gauss’s theorem. 

(20) Show that 

= r(^)r(e)r(g) F V5-a, s-a, <r; 1 

r(a)r((T+^)r(o+y) 3 2 |_ <r+p, a+y 


where a = 5 + e-a — p—y. 

Hence, using Dixon’s theorem, prove Watson’s theorem: 


3 * 2 


«> P, y ; l 

^(l+a +p), 2y; 


r(i>r(j+ y)m+ja+mrg-ia-w+y) 

r(i+*x)r(*+ mm -i« +y)m-iP+y) 


and, using Watson’s theorem, deduce Whipple’s theorem 
that, if a+j? = I, and e+5 = 2y +1, 

“K 1 ] 

= _ «r(g)r(«) _ 

i 2 '- 1 m *+i mi* + mw+map +w. 
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LEGENDRE FUNCTIONS 


§ 13. Legendre polynomials. If A is a fixed point with 
coordinates (a, /J, y) and P is the variable point (*, y, z), 
then if we denote the distance AP by R, we have 

R 2 = (x-a) 2 +(y~P) 2 +(z-y) 2 . 

Furthermore, we know from elementary considerations that 



is the gravitational potential at the point P due to a unit 
mass situated at the point A, and that this must be a 
particular solution of Laplace’s equation. 

In some circumstances it is desirable to expand ip in 
powers of r or r~ l where r = (x 2 +y 2 +z 2 ) i is the distance 
of P from O, the origin of coordinates. This expansion 
can be obtained by the use of Taylor’s theorem for func¬ 
tions of three variables but it is much more suitable to 
introduce the angle 0 between the directions CM, OP (cf. 
fig. 4) and write 

R 2 = r 2 +a 2 —2ar cos 6, 

The expression for \p then becomes 


*- 


yf(a 2 —2arfi-i-r 2 ) 


(13.1) 


where /i denotes cos 0, and this can be expanded in powers 
of rja when r<a and in powers of ajr when r>a. If we 


so 
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denote by PJji) the coefficient of h" in the expansion of 
(I - 2}ih+h 2 )~ i in ascending powers of h, i.e. if 


1 

V(l-2 nh+h 2 ) 


£ PMh\ 

B = 0 


( 13 . 2 ) 


then the potential function (13,1) can be expanded in 


the forms 



(13.3a) 

(13.3b) 


It is clear from the definition (13.2) that the coefficients 
P n (ft) are polynomials in ft. The first one or two can 
readily be calculated directly from the definition. By the 
binomial theorem we have 


(1-2^1 h+h 2 )'* 

= l+(-iX-2 nh+h 2 )+ (-*K-*> (_2;ifi + fe J ) 2 + ... 


= 1 +Hh+ K3 n 2 -1 )h 2 + M5/i 3 - 3 n)h 3 + ■■■ 
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SO that 


Po(ji) = 1. Piin) - /»> 

= i<3/* a -l)» i*30t) = KV-3/t). (13.4a) 

We sliall show below that, in the general case P„(ji) is a 
polynomial in n of degree «; it is called the Legendre 
polynomial of order n. 

The expression for the general polynomial can 
be derived by the method employed to obtain the simple 
expressions (13.4a). 

Expanding (1 — 2;iA + /i 2 ) -i by the binomial theorem 
we have 

(1-2 fih+h 2 r*= Z — (2/ift-A*) r 

f-o r! 

and the coefficient of If in this expansion is the coefficient 
of If in the expansion 

S ( -^(2 t th-h 2 y= t »h~hy-* 

= f ( tr (*>■ 

«-o V (i-H), (ii-c)! ’ 

since by example 9 (ii) of Chapter 1 

Now ihe coefficient of If in the expansion of 

t ~ 1) ° (2fih — h 2 y~* 

(n-e)! 

is 

w a> 

e!(«-2e)!' 
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and, by the duplication formula for the gamma function, 
n! 


(« —2g)! 


so that 


. 2~ 2e 


rQn+i-e)r(in+i-e) 


n! e = o e!(i-n) e \jl 2 J 
a result which may be written in the form 

P.(M) = ( -^P 2 ,F, (i-in, -in; i~n; A). (13.4b) 

Putting p = 1 in equation (13.2) and equating co¬ 
efficients of If we find that 


PJL D-l (13.5a) 

for ali values of n. Similarly if we put p = — 1 in (13.2) 
we derive the result 

PJt-D = (-D B 03.5b) 

which is a particular case of the result 

PA-n) = (-iyPM- (13.6) 

Equation (13.4) gives P„(cos 0) as a polynomial in cos 0 
of degree n so that it should be possible to express 
P„(cos 0) in terms of cosines of multiples of 0. Instead 
of attempting to do this by substituting the appropriate 
expression for cos r 0 in (13,4) we begin afresh with the 
definition (13.2). Writing (1 — 2 cos Oh+h 1 ) in the form 
(1 —e ,s A)(l —e~ ,6 h) we find that 

Z A(cos 0)h n = 

n = Q 

= f f nr+m*+$ +» e t(r-j 9) 

',-*.-0 r(i)r(i)r!.s! 
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Equating the coefficients of h" we find that 


Pn (cos 9 ) = 


y r(HQr(H»-0 c K2r-n)il 

rfo TO)r(i)r!(n-r)! 


Using the duplication formula, we see that 


rq + r)r(l+n-r) = 1 (2n-2r)!(2r)! 
rXi)r(i) 2 2 " r!(n —r)! ’ 


M3 


so that 


P n (cos 0) = 


l y (2n—2r)l(2r)l e ( ( 2 r-n)o 
2*" r -o (rl) 2 {(n—r)!} 2 


from which it follows immediately that 


^2«(COS fl) 

... {(2»)1} 2 

2 4 n (/ i !) 4 

and 


t 1,-1 

— V 

tB-l *->. 


(4)i —2r)!(2r)! 
'o (r!) 2 {(2rt —r)!} 2 


cos(2n—2r)0, 
(13.7) 


/> 2n+t (cos 6) 



2 4n+1 rt'o 


(4n + 2—2r)!(2r)! 
(r!) 2 {(2« +1 —r)!} 2 


cos (2)1 — 2r +1)0. (13.8) 


From these last two equations we may derive a general 
result of some importance. We may write 

P B (cos0)= y c r cos(n-2r)0, (13.9) 

fa 0 

where p = in or in—i according as n is even or odd. In 
particular 

P*( 1 >= f 

r = 0 
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I p n (cos0) ] £ 1, (18-10) 

10 


0-5 - 
P n (cm 0) 


Fig* 6, Variation of P n (cos 0) with 9. 
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so that from (13.5) 


Now, from equation (13.9) we have 

| ^(COS 0) | ^ £ c r . 

r * 0 

and therefore 

| ^(cos fl) | g 1. (13.10) 

The variation of P„(ji) with ft for a few values of n is 
shown in fig. 5. Since, in most physical problems, the 
Legendre polynomial involved is usually P„(cos 0) we have 
shown in fig. 6 the variation of this function with 6. 
Numerical values may be obtained front Tables of Associated 
Legendre Functions (Columbia University Press, 1945). 

§ 14. Recurrence relations for the Legendre polynomials. 
If we differentiate both sides of equation (13.2) with respect 
to b we have 


fi-h 




(1 —2/jft + ft 2 ) 3 
which may be written in the form 


nh"~ t P„{]£) t |/<!<], 


(ft-h) £ VPfn) = (l-2yh+h i ) £ nh'-'Pffi). (14.1) 

ll — 0 ff as 0 

Equating coefficients of /)" we have 

i&M-P.-M 

= (n +1 )P„ + 1 00 ~ 2/u<P„(/0 + (n -1 )P„ _! 00. 

which reduces to 


<" + DP. + iW- (2n+l) f iP H ( fl )+nP^ l ( t t) = 0. (14.2) 

This relation has been proved to hold for | fi j < 1 but 
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since the left-hand side is a polynomial in /t, it must hold 
for all values of ft. 

On the other hand, if we differentiate both sides of 
equation (13.2) with respect to ft we obtain the relation 


h 

{l-2fth+hy 


1 h"PM- 

n « 0 


(14.3) 


Combining equations (14.1) and (14.3), we have 
(n-h) £ k'PXfi)m £ nVPJM, 

It - O It SB Q 

so that equating the coefficients of h n we obtain the relation 
PP'M-P' n - iG») = nPfn), (14.4) 

and since each side is a polynomial in (i this relation holds 
for all values of / i. 

If now we differentiate equation (14.2) with respect to 
H we obtain the relation 

(n + l)P; +l 0i)-(2« + l)P,(Ai) 

-(2n+i)tiP;Q t }+nP^ t (jt) = 0. (14.5) 
Eliminating P'„(}i) from (14.4) and (14.5) we see that 

p; + i 00 -p;-i 0 o = ( 2 «+i)p n (^. (14.Q 

Subtracting (14,4) from (14,6) we obtain the recurrence 
relation 

PU i(M) - - (n +1 )P a (ji). (14.7) 

The differentiations with respect to h and /i under the 
summation sign is justified by the fact that the series on 
the right-hand side of equation (13,2) is uniformly con¬ 
vergent for ail real or complex values of h and pi which 
satisfy the relation j /i | £ 1, [ // | c^/2—L 

§ 15. The formulae of Murphy and Rodrigues, From 

E 
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the expansion (13.2) it follows immediately that 


I A"£;MW = £(l-2/ifc +ft 2 )"* 

n=eo dp dp 

= 2 'A r r<r+ ^ (l- 2 /ift+ft J )~ , ~*. 

Hi) 

Substituting the value /i = I we see that 

i P' r) (l)ft" = 2' fc'(l - h )- (2r+11 
»= 0 r(?) 

_ 2 r // y f~(l+2r+s) ^ 

TQ) .-o I~(l+2r)s! 
where Pi r) (p) denotes d'P 0 {^)(dn T . 

Equating coefficients of h n we see that ^'(l) = 0 if r>/i, 
as is obvious from the fact that P„(/j) is a polynomial of 
degree « in fi, and that 


_ 2 r HI + H + r) 

m r(l + 2r)(n-r)!’ 

From the duplication formula for the gamma function 


2T(r+j) = J_ = 1 

r(i)r(l+2r) r!2 r (1),2'’ 

and from example 9 (iii) of Chapter I, 


so that 


r(l + n+ r) 

(n-r)} 


= (-!)'(«+ lU-n) r . 




(n + l) f (-«) r 


(lUr 


(15.1) 
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Now, by Taylor’s theorem 


mw = 


GO 



(jLjr p(ry {l) 

r! 


Substituting the expression (15.1) in this expansion we 
obtain the relation 

rM - (lj r r! (T 

which gives Murphy’s formula, 

MW = 2F1 (-». «+i; i; (is. 2 ) 

for the Legendre polynomial P„(p). 

If now we put « = I in example 5 of Chapter II we see 
that equation (15.2) is equivalent to 


MW = 


1 d" 
2 V dff 




(15.3) 


which is Rodrigues’ formula for the Legendre polynomial. 

Rodrigues’ formula is of great use in the evaluation 
of definite integrals involving Legendre polynomials. 
Consider, for instance, the integral 


-L 

ila 1 

-r 

2 n! J -, 


f(x)P„(x)dx. 


(15.4) 


By Rodrigues’ formula we may write this integral as 

d" ^ 


/(x)f-(x 2 -iKx t 
dx 


and an integration by parts gives 


1 f d n 




it 
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The square bracket vanishes at both limits so that we have 

d”~ l 


-JLf 1 /', 

2"n! J-i 


dx 1 


Continuing this process we find that 

/ - J l (* 2 - m M {x)dx. (15.5) 

For example if j{x) — P m {x\ m<n f — 0 and so 
/ = 0. In other words 


L 


Pm{x)P„(x)dx = 0, (m ¥=n). 


(15.6) 


If Ax) = P n (x) then 


/ w ( jc> = 1 dl ] (jc 2 - ir 

2 n it ! dx 2 " v 


_ (2n)l 
" 2"n! 

Hence 

J_ t {«*)}’<** = 2>!)a J _ t (1 ~ x2)Vx 

_ (2»)i r(j)r(n + l) 

"" 2 2 "(n!) 2 n» + J) 

Making use of the duplication formula for the gamma 
function we can reduce this to the form 

f { P r,( X )Y dx = (15.7) 

J-i 2 h+ 1 

A convenient way of combining the results (15.6) and 
(15.7) is to write 

I P m (x)Pn(x)dx= (15.8) 

J_t 2n + l 
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Sis 


where S m§ „ is the Kronecker delta which takes the value 
0 if m £ n and the value 1 if m = n. 

Similarly ifyfje) = x m , where m is a positive integer, then 


f M (x) = 

and hence, if m>n, 


r(m+l) 


T(m —n+1) 

0, 


x m ", if m'^n, 
if m<n. 


L 


x m P n (x)dx = 


2 n r(m 


r( m + l) f 1 
-n + l)n! J_ t 


x m -"(l-x 2 rdx. 


If m—n is an odd integer the integral on the right is zero 
while if m—n is an even integer it has the value 




2 xr-\i-x 2 rdx = 


j Wv _ r(^n-jfi+i)rQi+i) 


f(im + in + ^) 


so that, if m is an integer. 


Jl 


x m P„(x)dx 


0, 


if men. 


m T(^m —jp m _ f |^o is even, (15.9} 


2"(m - n )! f(} m + \n +1) 

0, if m — n>0 is odd. 

If m — n the result is 


j 1 x”P„(x)4x = i J‘ (1 -x 2 T 


dx 


1 r(i)T(n+l) 

T r(n+}) 
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which, on account of the duplication formula, is equivalent 


to 



x a P„(x)dx = 


2 " + 1 ( h !) 2 

(2n+l)l 


(15.10) 


§ 16. Series of Legendre polynomials. In certain 
problems of potential theory it is desirable to be able to 
express a given function in the form of a series of Legendre 
polynomials. We can readily show that this is possible 
in the case in which the given function is a simple poly¬ 
nomial. For example, from the equations (13.4a) we have 

i = PM, 
a = PM, 

v 1 = i+iPi GO = iFo00+^ 2 (/i), 

n 3 = ifi+iP 3 (p) = jPM+jPM, 
so that any cubic c 0 ft 3 +CiH 2 + c t ft+c 3 can be written as 
the series 

^fPM+ 2 -^P 2 (M)+ (y+ c *)' , i(mH (j 1 +c^PM- 

It is obvious that we could proceed in this way for a 
polynomial of any given degree n though if n were large 
the arithmetic involved might become cumbersome. Since 
P„(n) is a polynomial of degree n in it emerges as a 
result of an extension of the above argument that any 
polynomial of degree n in p can be expressed as a series 
of the type 

i c r P r (n), {-IZfiZl )• (16.1) 

T ** 0 

The problem which now arises is that of expressing any 
function/(/*), defined in the interval — as a series 

of Legendre functions of the form 

Z CJM- 

r “ 0 


(16.2) 


S1< 
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If it is assumed that the infinite series (16.2) converges 
uniformly in the range (—1, 1) to the sum /(ft), we may 
multiply the terms of the series by P„(p) and integrate 
term by term with respect to ft over the range (—1, 1) 
to obtain the relation 

| f{fi)P„(ti)dn = c r | PMPn(p)dn 

and by equation (15.8) the sum on the right-hand side is 
equal to 

y C - $ = ?Sl 

r?o '2n+l ,,H 2« + 1 

which shows that the series 

r Z o («+i)^(P)J l /(v)P.(vMv (16.3) 

converges uniformly to the sumJlp) in the range (—1, I)* 
The series (16.3) is called the Legendre series of the 
function /(/0« We shall not discuss here the conditions 
which must be satisfied by the function f[p) if this series 
is to be uniformly convergent; for such a discussion the 
reader is referred to Chapter VII of E« W. Hobson, The 
Theory of Spherical and Ellipsoidal Harmonics (Cambridge 
University Press, 1931). 

The possibility of expanding a function in the form of a 
series of type (16.2) is a consequence of the relation (15.6). 
In the theory of special functions we frequently encounter 
sequences of functions $,(*), ^X Y )> **■ which 

have the property 

| 4>Mtf> K (x)dx = 0, (m*n). (16.4) 

We then say that the functions (r - 1,2, ...), form 
an orthogonal sequence for the interval (a, b). If, in 
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addition, the functions are such that 

J {Ux)} 2 dx = 1, (16.5) 

for all values of /r, we say that the functions of the sequence 
arc normalised, and form an orthonormai set. Given a 
set of orthogonal functions it is obviously a simple matter to 
construct a normalised sequence. For example we see from 
equation (15*8) that the sequence of functions P n (x) t 
( n — 0* U 2, ,,,) is orthogonal but not normalised. By 
multiplying each function by we find that the 

functions 

(16,6) 

form a sequence of normalised orthogonal functions in 
the interval (—1, 1), 

There is another property of importance which such a 
sequence of functions may possess. If there is no integrable 
function ^(x), different from zero, such that 

| Hx)<j>„(x)dx = 0 (16.7) 

for all values of n we say that the sequence is a complete 
orthogonal sequence. It may be shown in the case of 
the functions (16.6) f by considering the Fourier coefficients 
of \p{xjn)P„(xjn) in ( — te, n) that if (16.7) holds this function 
is a null-function and hence that tf/(x) is a null-function 
in (-1, 1). In other words it can be shown that the 
functions (16.6) form a complete sequence of normalised 
orthogonal functions. 

§ 17. Legendre’s differential equation. If we write 
v = ( n 2 ~~l) n 

then it is readily shown that 

(1—;< 2 ) ^ +2finv — 0 
dfi 

+ E, W, Hobsora, op , eit.* p. 40. 
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and if Wd differentiate this equation «+ f times using 
Leibnitz’s theorem we find that 




d” +2 o 


-2(i 


d n+, u 
dpi** 1 


d”v 


which when written in the form 


{a-*>£-*i«u + D }(£)- 0 

shows that d"(ji 2 — \) n ld(i n is a solution of the differential 
equation 


,.2s dl y 


(1-/0 ~i -2#t / +n(M + l)y = ° (17.1) 

a\k dfi 


so that we conclude from Rodrigues" formula (15,3) that 
when n is an integer P n (ji) is one solution of the equation 
(17,1). This equation, which we shall now consider in a 
little more detail, is called Legendre’s differential equation. 
We saw in example 1 of Chapter I how such an equation 
arises in the solution of Laplace’s equation when solutions 
of the type i?(r)0(cos 0), Le. m = 0, are sought. 

It is obvious by inspection that the point fi = 0 is an 
m dinary point of the equation (17 J)* Writing the equation 
in the form 


^ d 2 y 


(tt-lVji +(#*-!) 

dft 


2ji_dy 

fi + i d(i 


n(« + l) 
/r+1 


Qt-l)y = Q 


and observing that in the notation of equation (3.1) with 
a=l, 

m= 

1+iOi-l) 

/I—1 

1+*(/<-!> 


= -Hn+iKOt-D- 


90*)= -$«(« + !) 
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we see that ft - I is a regular singular point with indicia! 
equation q 2 = 0. 

Furthermore in the notation of equation (4.2) 


1 -fT 1 -ft 1 

so that as fi-*co 

«0*)~ -» P(m)~ - 

ft n 

showing that the point ft = oo is a regular singular point 
with indicia) equation {g-(n+l)} (g + n) = 0, We thus 
see that the equation is defined by the scheme 


y = P 


-1 

0 

0 


oo 1 

n + 1 0 ft . 
-n 0 j 


(17.2) 


If, however, we put 

* = K1-/0 

in equation (17.1) we find that it reduces to the form 

*0 -+(1 - 2jc > T- + "(« + % - 0 . ( 17 . 3 ) 

ax 1 dx 


which is equation (8.1) with a = « + l, fi = —n and 
7 = 1, so that the scheme (17.2) is equivalent to the scheme 


y = P 


0 oo 1 

0 n + 1 0 -J — \ft 
0 -n 0 


(17.4) 


It should be noticed that the values along the top row are 
those assumed by not by ft. 

We consider first the solution corresponding to the 
singular point at infinity. We write 

yt(p) = v° X <y*~ v 

v « 0 
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which on substitution into (17.1) leads to the recurrence 
relation 


(n — v + 2)(n - v + l)c v _ 2 = — v(2n +1 — v)c v . 
On taking c Q = 1 we obtain the solution 


yiO) = v n - 


n(H-l) 

2. (2n —1) 


w ”2 


n(n-l)(H-2)(H-3) -_ 4 + 
2.4.(2n-l)(2n-3) 


, H y., (~inxi-i«) i 

y 1 1 ■(!—«) P 2 

■ (-*»X“I" + 1)(W")(W« + 1) 1 1 

l-2(i-»Xi-» + l) f " J 


which may be written in the form 

JhOi) = i-«; pj. (17.5) 

Also, if we write for the second solution 


yi(n) = h " 1 I d,fi “ 

v = 0 


we find that 
FiOO = /<""' 


, (n + !)(« + 2) -„-3 
2.(2n + 3) 


( n+l)(« + 2Kn + 3)(n+4) 
2.4.(2n + 3)(2n + 5) 


= • (*«+*> In + Un + bpj, 


- s + ... 
(17.6) 


provided n is any number other than a negative integer 
or half a negative integer. 

These solutions are valid for all values of n for which 
the 1 F i series have a meaning, not only for integral values 
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of n. If n is an integer the series y L (tt) terminates—in 
other words, y t (ji) is a polynomial of degree n in p. If we 
multiply this polynomial by 


( 2 «)! 

2 "( n !) 2 


we obtain the Legendre polynomial of degree « (equation 
(13.4b) above). 

On the other hand the series for y 2 {pt) does not terminate 
when n> — 1 so there is no point in restricting n to be 
an integer. This series solution when multiplied by a 
factor 


IW( H + 1) 

2" +1 r(n + i) 

gives the function 

QM 

- (*" + !• + n + i; I). (17.7) 

The function 
p M . 4SS+JL 


2“{r< 




is a solution of the Legendre equation (17.4) even when 
«is not an integer and it reduces to the Legendre polynomial 
when n is an integer. We shall continue to denote it by 
PJjt) but when n is not an integer shall refer to it as the 
Legendre function of the first kind of degree «. The 
function QM defined by equation (17.7) will be referred 
to as the Legendre function of the second kind of degree n; 
even when n is an integer it is not a polynomial. 

With these definitions we may write the solution of 
Legendre’s equation (17.1) as 


y~APM+BQM 


(17.9) 
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(j n |>1). In some problems we know that the solution 
should be a polynomial in jt; in that case we must take the 
solution to the form y = AP„(ft). 

The variation of Q„ (/t) with n may be computed easily 
from equation (17.7) when /o 1. Tables calculated in this 
way are contained in the volume quoted at the end of 
§ 13. The following fig. 7, which was prepared from these 



12 3 4 


Fig. 7. Variation of (LOO with /i. 

tables, shows the variation of Q„(ft) with ft for a few values 
of n. 

Since Legendre’s equation lias a regular singular point 
at ft = I we may on the basis of equation (3.8) take the 
second solution of Legendre’s equation to be proportional to 

PM !ogOi-l)+Ic r (/< —l) r . 

The coefficients c r can now be obtained by substituting 
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this expression into the differential equation (17.1) and 
equating to zero coefficients of successive powers of 
(/t-1). Instead of proceeding in this way we shall derive 
this second solution by means of a method due to F. E. 
Neumann. 


§ 18. Neumann’s formula for the Legendre functions. 
Let us now consider the integral 


r 

J-i p-Z ’ 

where |/r|>l, and n is a positive integer. Expanding 
the denominator by the binomial theorem, we have, for 
the value of the integral, the series 

From equation (15.9), it follows that the integrals in this 
series are zero if or if s = n + 2r where r is a positive 
integer, so that the above series is equivalent to 

Z f 

r=0 ft" 1 " J_i 

which, by the same formula, is equal to 

J_ 5. (n + 2r)!r(r + j) /lY 

fi n+i r«o 2"(2r)!r(n+r+*) W' 

Now from the duplication formula 

(n + 2r)!r(^ + r) _ r($n44 + r)r($n + l +r) 

2\2r)\ r! 

so that the series reduces to 


1 f r (jH + i + r)r(in + l+r) /l V 

/i" +1 * ™ o r(n + i+r)rl WJ 
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which is equal to 


r(in+i)rq n +i) 

T(n + i)/ + 1 


2 Fi in+1; n+3; 



Noticing that 


m«+i)r(±»+i)/r(n+*) = r(i)r(«+i)/2T(n+i), 

and comparing the result with equation (17.7) we see 
that this series is merely 2Q„(p). Hence we have shown 
that if | n | > 1 

e,oo=if (i8.i) 

J-i p-z 

a result which is known as Neumann’s formula. Equation 
(18,1) holds not only for real values of n greater than 1, 
but for all values of n which are not real. For this reason 
equation (18.1) may be regarded as defining the second 
solution, Q„(p), of Legendre’s equation. 

Certain related formulae, due to MacRobert, follow 
readily from this result. If | ;i |>1 and m is a positive 
integer then 

r r 

J-i /t-f J-1 p~Z 

and the integral on the right is equivalent to the finite sum 

*lV _I_r P modi. (i8.2) 

•■ = 0 J-l 

If it follows from equation (15.9) that each term of 
this series vanishes so that we have 


P m QM = i f‘ (18.3) 

J~. p-Z 

provided m, n are integers and 

On the other hand, if m = w +1 the series (18,2) reduces 
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to the single term 


so that 


>J> 


■ffW 


2"(n !) 2 
(2n+ 1)!’ 




= i f 1 
J-i P~( 


(18.4) 


Z (2n + l)! 

Now if m is an integer P„{p) is a polynomial of degree 
m in fi so that it follows from (18.3) by finite summation 
that, if m£«, 

PMQM -i r ^ )P ”P dt (18.5) 

J-I p-z 

Similarly from equation (18.4) and the definition of P„ +1 (p) 
we have the formula 

^lOO&GO -i (*' — . (18.6) 

J -1 P~Z n + 1 

ff we replace n in equation (18.5) by /i-l-l and m by n 
and subtract from equation (18.6) we obtain the relation 


p n + i(p)QM-P n (ft)Q n+ 100 = -L (18.7) 

» + l 

Other formulae of a similar nature are contained in ex. 24 
below. 

We shall now make use of Neumann’s formula to derive 
the form of the second solution of Legendre’s equation in 
the neighbourhood of the points fi = +1. From (18.J) 
we have the result that if fi> 1, 

Qn 00 = iW log ~ - W n _ ,04 ( 18 . 8 ) 

p-\ 

where W n _ ,(//} denotes the integral 

i f m=m d( 

J-i 
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§ 18 


Now when n is an integer, P n (fi) is a polynomial of degree 
n in /t so that {P n (p) — )}/0* “ £) is a polynomial of 

degree h— 1 in /i. Hence W m -i(pL) is a polynomial of 
degree h— 1 in fi. 

If we substitute from (18.8) into Legendre’s equation 
(17.1) we find that satisfies the differential equation 

(1 - p 2 )w;. I - 2 ^'_ 1 + / I (/ I +I)w;. l = 2P&1). (18.9) 

Now since W„-!(p) is a polynomial of degree n-1 in p 
we may write (of. § 16 above) 

w;-iG0» 1 C'PM- 

r = 0 

Using the fact that 

Pi 00= i (2n-4r~l)F„. Jr „ II 

r = 0 

where p = i( n— 1) or in—1 according as n is odd or 
even (which follows from equation (14.6)), and the result 


(1 - p 2 )P?(ii) - 2/(/>;0i) + „(/, +1 )P r (p) 

= (n-r){n+r+l)P r (p), 
we find, from equation (18.9), that c n _ 2s = 0 t (s — 0, ...,/?) 
and that 

_ 2n —4s—1 
C " _2l_l (2s+lXn-s)' 

Substituting these values in (18.8) wc obtain the formula 


QM = ±p„(p) log'll 
/<-i 



j- 0 


2>t—4s- 1 
(2s + l)(n — s) 


P*-2m-1<J*)‘ 

(18.10) 


Another expression for Q n (ti) may be derived from the 
fact that both PJjt) and QJji) are solutions of Legendre’s 
equation (17,1) so that 


<m f {a -f)pm-Pn(p) t {(i-V)Q,;oo} 

dpi dp 

F 
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which is equivalent to 


~ [v-P 2 ){QMPfo)-p n mm] - o, 

showing that 

(i - -PMW} - c, (i 8 .ii) 

where C is a constant. Now from equations (17.7) and 
(17.8), we can readily show that for large values of ft 


P'MQ&i)~ 


n j_ 
2m +1 ft 2 ' 


Q' n (n)P n (fi )— 


M + l 

2 m 4- 1 


1 


ft 


z 


so that as fi-*<x> the left-hand side of (18.11) tends to 
— 1, showing that C = — I. Writing (18.11) in the form 


i fe„(»)i -i 

dn\pm (fi z -\){P n (p)¥ 
and noting from (17.7) that as ^-*oo we have 


QM = P„(ft) 


r 

Ju 


dt, 




(18.12) 


§ 19. Recurrence relations for the function Q„(ft). 
Recurrence relations for the Legendre function of the 
second kind can be derived from Neumann’s formula (18.1) 
and the corresponding recurrence relations for the Legendre 
polynomials P„(ft), From the recurrence relation (14.2) 
and Neumann’s formula we have 


(»+i)e, + i0t)+«Q (1 - 1 0i)=(«+i) T 

J-i p-z 


Now 


i J* ^ dC - nQMO-i P„(m- 
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If we write the second term on the right as J P 0 (!;)P„(Z)cU; 

we see from (15.6) that it vanishes if n =4 0. Hence we 
have 

in +D&.+tO*M2«+DtfWMH i(M) = 0 O 9 - 5 ) 

showing that the functions Q„(f 0 for three consecutive 
values of n satisfy a relation of the same form as that for 
the functions P„(ft) (equation (14.2) above). 

From Neumann’s formula (18.1) we have 

Pn ^\ dt 


and if we integrate by parts on the right-hand side we 


find that 

e»=i{^- + 



Hence 


=i r i K+M-rua # 

J-t p-z 

f 

i-tfi-z 


= i(2n + \) 


by virtue of equation (14.6). The integral on the right is 
2 Q n (p) by Neumann’s formula so that, finally, 

ql+m-ql-iW = ( 211 + 1 ) 6 ^ ( i 9 . 2 ) 


§ 20* The use of Legendre functions in potential theory* 
In potential theory we have frequently to determine 
solutions of Laplace's equation V 2 ij; — 0 which satisfy 
certain prescribed boundary conditions. If we have a 
problem in which the natural boundaries are spheres with 
centre at the origin of coordinates it is natural to employ 
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polar coordinates r, 0, 0. In cases in wliich there is sym¬ 
metry about the polar axis ip will not depend on <p so we 
may write ip = \f/(r, 0 ). It then follows from example 1 
of Chapter I that 

will be a solution of Laplace’s equation provided that v a 
is a solution of Legendre’s equation (17.1). Taking to 
be P„(cos 0)+C K QJpo$ 0 ) we see that we may write 

Hr, 6)= £ (A n r>+B„r-’>- l )P n {cosO) 

n - 0 

+ t (C„r" + D„r-"- 1 )Q„(cos0) (20.1) 

n = 0 

where the quantities A„, B„, C„, D„(n = 0, 1, 2, ...)are all 
constants. 

Now it is obvious from equation (18.8) that Q B (cos 0) 
is infinite when 0 = 0°, and we know on physical grounds 
that in the case of spherical boundaries ip remains finite 
along the axis 0 = 0. Hence we must take C„ = D n = 0 
for all values of n and obtain the potential function 

*= £ (A„ r "+ B„r~ n ~ 1 )P „(cos 0) (20.2) 

A » O 

which Is valid as long as r is neither zero nor infinite, i.e, 
if a£r£b where a and b are finite and non-zero. If the 
region under discussion is the interior of a sphere, i.e. if 
then to avoid ^ becoming infinite we must take 
B n to be zero to give 

</'= £ A n r"P „(cos 6). (20.3) 

n ^ 0 

On the other hand, if the region being considered lies entirely 
outside this sphere, we must take 

£ B n r' l '- l P n ( cosfl). 

n ^ 0 


(20.4) 


S20 LEGENDRE FUNCTIONS 77 

Examples of the use of the solutions (20.2, 3, 4) in potential 
theory are given in Coulson’s Electricity (Oliver & Boyd. 
1948) 75-78, 80; a further example is given below. 

The Legendre functions of the second kind, Q„(cos 0), 
which are absent from problems involving spherical 
boundaries, enter into the expressions for potential func¬ 
tions appropriate to the space between two coaxial cones. 
lfO<a<0</?<7t we must take a solution of the form (20.1). 
Suppose, for example, that \p = 0 on 0 = «, and \p = Iv" 
on 0 — ft then we must have 

A„P„(cos <x)+C n Q a (cos a) = 0 

and 

A„P n (cos P) +C n Q n (cos p) = B„ — D n = 0, 

the latter results following from the fact tliat if a # (S, 
Q n (cos a)P B (cos p)-P n (cos a)£?„(co5 //) does not vanish. 
Solving these equations for A„ and C„ and inserting the 
solutions in equation (20,!) we find that in the space 
between the two cones 

” j Q„(cos a)P B (cos 0)-P n {cos «)Q n (cos 9) j > 

«-o " 1 q„(cos a)P„(oos P) - P n (cos a)S„(cos /?) j 

To illustrate the use of the solution (20.1) and of some 
of the properties of Legendre functions we shall now 
consider the problem in which an insulated conducting 
sphere of radius a is placed with its centre at the origin 
of coordinates in an electric field whose potential is known 
to be 

£ a nr n P »(cos (?) (20.6) 

it = i 

and we wish to determine the force on the sphere. The 
conditions to be satisfied by the potential functions i p are 
(i) that ^ is a solution of Laplace's equation; (ii) that \j; 
has the form (20,6) for large values of r; (ill) \j/ = 0 on 
r — a. 
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The conditions (i) and (ii) are satisfied if we take 

0 = /*-(«» fl) 


We 


and (iii) is satisfied if we write B n - —« H a 2n * t , 
therefore have 

*>= n E -^rpj P„(cos 0). 

The surface density of charge on the conductor is 

ff “ - = - j- £ (2n + l)a"" 1 « B / > „(cos0) 

4?t \OP/, = a 4)1 b = I 

and since the force per unit area on the conductor is 27 r<r 2 
the resultant force on the sphere is in the 0 = 0 direction, 
and is of magnitude 


F= P Ina 1 
Jo 


Ina 1 sin 0 cos Odd 


= i« 2 I E a m a « a " 


+n-2i 


(20.7) 


n ~ 1 m — 1 

where I mn denotes the integral 


J, 


(2n + l)(2m +1) cos 0 sin 0 P m (cos 0)P „(cos 0)d0. 


Changing the variable of integration to fi = cos 0 and 
using the recurrence relation (14.2) we find that 


f mb = (2m +1) 




PMP a+ [GO d}l 


■"I. 




and by the orthogonality property (15.8) this reduces to 


5 21 

the form 
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7™ = 2(n +1 )<i B , „+i + 2n<5«i. «-i- (20.8) 

Substituting from (20.8) into (20.7) we find that the total 
force on the sphere is 

P= £ (« + l)ov* B+1 a 2 " +1 . 

fi “ l 

§ 21. Legendre’s associated functions. We saw in 
example 1 of Chapter I that the solution of Laplace’s 
equation in spherical polar coordinates reduces to the 
solution of the ordinary differential equation, 

(1 -/< 2 ) 7^ -: 2 o ~ + {«(« + 0- r~y| © = 0 (21.1) 

dfi dft (. l~fi) 

which reduces to Legendre’s equation when m = 0. This 
equation is known as Legendre’s associated equation. To 
solve this equation we may write 

0 = Oi 2 -l ( 21 . 2 ) 

Substituting this expression in the differential equation we 
find that y satisfies the equation 

- 2 ( 1 +(n + m)0t-m + l)y = 0 
dfr dft 

and differentiating this equation m times with respect to fi 
by Leibnitz’s theorem we find that 

{»-*>£-*.£ + «i. + l>}g- 0 (2.-3) 

showing that if d m yjdp m is a solution of Legendre’s equation 
(17.1) the function 0, defined by equation (21.2), is a 
solution of Legendre’s associated equation (21.1). 

Similarly if we put 0 = (n 2 — \)* m y in equation (21.1) 
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we find (hat 

(1—/t 2 ) 3-7 -2(1 +m)n ^ +(n — ;nX»+m + \)y = 0. (21.4) 
tf/r dft 

If now we differentiate equation (17.1) m times with 
respect to /* we obtain the equation 

j(l -/i 2 ) fs - 2(1 +nt)fi 4- +(M ~ m)(« + m +1)1^ = 0 
( d/i d/i J d/i ffl 

showing that if \{n) is a solution of Legendre’s equation 
then 

(21-5) 

tift" 

is a solution of Legendre’s associated equation ( 21 . 1 ). 
Taking the two solutions of Legendre’s equation to be 
and it follows from (21.5) that the functions 

™ = (/i 2 -l) im ^^, 

d(i m 

oroo= ( 2 i. 6 ) 

a/i 

are solutions of Legendre's associated equation. As a 
consequence of equation (21,3) we see that so also are the 
functions 

K m (n) - (/< 2 -i) _im J" j v ... j' PMmr (2i.7) 

and 

Qn m (n)=(^-iy im r p... r (uotdcr. ( 2 i. 8 ) 

J oo J 00 J oo 

It is an immediate generalisation of (6.5) that 

P‘, T, *) = a F,(a + m, J 3+m; y + m; x) 

dx m (y) m 
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so that using Murphy’s form (15.2) for P n (ft) we see that 


dft m (— 2)“f n! 


jf, pn-n, n + m+I; m + 1 ; y^. 


Hence P%(n) as defined in (21.6) may be written in the 
form 

™ “ + ° i2 “ l)im 

x 2 F, ^m-n, n+m+1; m+1; y^. (21-9) 

Other expressions for the first of the two functions (21.6) 
can be easily derived. If we make use of the result (7.4) 
we see that 

™ = ~5^— t T Tv^ + ir lm (/i-1) 

2 n mir(n-m +1) 


z Fj |m-n, -n; m + 1 ; ( 21 . 


10 ) 


and similarly, if we make use of relation (7.6) we may 
derive the expression 

= r(w+m+l) //*-iV M 

" m tr(n — m + 1 ) \n +1 / 

x 2 F, pi +1, -n; m + 1; y^. (21H) 

From Rodrigues' formula (15.3) we derive the simple 
expression 

•fTOO * —t (m 2— i)*" 


Jm+B 


2 "n! 


dft' 


m + n 


(F 2 -l )"• (21.12) 
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The simple differentiation 
d” -«-t-ar ( ]yH I~(m + n + l+2r) 

r(n + l+2r) 

may (because of the duplication formula) be written in 
the form 


rv, | n <*" + 

v n"+' +2r j 

= (- l)"T(m +n +1) (i»» + i«+l) r gm + in+ 1 ) r 

v * x 7 u m+n+l+2r 

showing that, by term-by-term differentiation of the solution 
(17.7) of Legendre’s equation, we obtain the solution 




r(j)r(m + n + l) 1, 2_ 

2 " +1 r(n+|) # W ’ 


x 2 -fi im+iu + 1; n+$; 4^ 


(21.13) 


of Legendre’s associated equation. 

Solutions of the type (2L7) and (21.8) can be derived 
in a similar fashion. Using the result 




r! 

(m+l)r 



in equation (21*7), with Murphy’s expression (15.2) for 
P n (fi), we derive the expression 


“(M-ira*i (-m»+l;m + li (21. 

n i n 

-I 


14) 


for 
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™ = (■"• n+l1 m+u ¥)■ (2U5) 

The solution provided by Rodrigues’ formula (15.3) can 
obviously be written as 

K m (v) = ~ J; -r^O* 2 -!)". (21.16) 

2 ill dft m 


In a similar way the solution 


e™=(-ir 


rtOTi-m+i) 

2' , + T(n + i) ‘ fr m + l 


x 2 F, 




}m + i, in~im +1 


;n+ *y) 


is derived from equation (21.8) and the result 


(21.17) 


r<»+i)P 

J 00 J co J o& ^ 

- (- 1)T(„ - m+1) 
used in equation (17*7). 

The four functions Pn(p)> Qn(v)> K m (v\ Qn m (p) defined 
equations (21.9), (21,13), (21.15) and (21,17) respectively 
are therefore solutions of Legendre’s associated equation. 
They are known as Legendre’s associated functions- Al¬ 
though the expressions above have been found by assuming 
m and n to be integers it is readily shown that the solutions 
quoted are valid even when m and n are not integers. 
Since Legendre’s associated equation is of the second 
degree it follows that only two of these four functions are 
linearly independent, and that the other two may be expressed 
simply in terms of them. It follows immediately from 
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equations (21.11) and (21.15) that if m, n are integers 

K m (M) = (21.18) 

r(«+ m+ 1 ) 

Furthermore, if we apply the result (7.6) to the hyper¬ 
geometric series on the right-hand side of equation (21.17) 
we find that 


2,™ = (-!)" 


IWQi-m+l) 

2 " +l r(n + S) 


/i"'" _ "“ 1 0i 2 _i) im 


x 2 ^i im+in+1; n + i; 


which, on comparison with equation (21.13), reveals the 
relation 


Qr<ji) = 


F(h — m + 1) 
r(n + m+1) 




(21.19) 


It is now a simple matter to prove that when m and n 
are integers, and m is fixed, the polynomials P%(jt) form 
an orthogonal sequence for the interval (-1, 1). Making 
use of the results (21.18), (21.12) and (21.16) we find that 





F( n + m+1) I_ 

r(n —m + 1 ) 2" +n ’/i!/j'! 


L 

and after inlegrating by 
on the right reduces to 


J rt-m jn' + m 

parts n—m times the expression 


F(n + m + 1) (-l)"" 1 " f* 
F(n-m + l) 2" + "'n!n'! 




<T + "‘ 

d/i n+ "' 


(M 2 ~iy'dn. 


This integral is evaluated by the method used at the end 
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of 5 15 and we find that 

f 1 PWPttM* = C-0 m z~7 K ( 2I - 2 °) 

J_, r(n-m + l) 2n + l 

In many physical problems /.i = cos 0 so that — 1 Sp= 1- 
It is then not always convenient to have a factor of the 
form ()i 2 — l) lm . We use instead Ferrer’s function 


™ = ( (21.21) 

With this notation we may write (21.20) in the form 


i: 


COOTCOMi = p / K (21 -22) 
i r(n-m + l) 2n + l 


The other formulae are amended similarly. 


§ 22, Integral expression for the associated Legendre 
function. If we assume Cauchy’s theorem in the form | 


1_ 

2ni 


f 

J cC-ti 


where J[Q is an analytic function of the complex variable 
C in a certain domain R which includes the point ( = /i 
and where the integral is taken along a closed contour C 
which includes £ = fi and lies wholly within the domain R , 
then by differentiating both sides of the equation r times 
with respect to fi we obtain the result 


_ si r m dc 
dn r 2ai J c (C-p) r+1 ‘ 


( 22 . 1 ) 


Substituting m+n for r and (f 2 — 1)" for AO in this equation 


+ E* G* Phillips, Functions of a Complex Variable (Oliver and Boyd, 
1940), p. 93. 
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we find that, as a result of equation (21.1), 


2". n! 2ni 


' f (C 3 -!)" 
J c (c-A»r + " H 


df. (22.2) 


If /i>0 we may take the contour C to be the circle 

Integrating round this contour we obtain from equation 
(22.2) the equation 

1 P {/i+V^-l) cos w-!//)}" C0S {m<f>)d<j> 

2tcJ 0 sin 


cos 


(j#V?0i) 


(n+m)! stn 

from which follows immediately the Fourier expansion 
{n + y/in 2 - 1) cos (^-^)} n 
n! 


= P n (^)+2 £ 


: P”(/i)cos»i(i^-0), (22.3) 


m = i (n+m)! 

Changing/i to — (ji+ 1) we obtain the expansion 

0*' +VCm'* - 1) COS 

= W)+2 t P^/Ocos (22.4) 

»t = 1 tt! 

Applying ParsevaPs theorem for Fourier series to the 
scries (22.3) and (22.4) we find that the series 

W W) + 2 2 (- I) ni P; n (;i)P™0(') COS m<p 

m = I (/!-}• »l)! 

converges to the sum 


2* J-« O'' 


{z'+VV-l)cos 


+ - s /(/i' 2 -l) cos V>} 


n+l V 
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This integral may be evaluated by means of Cauchy's 
theorem.t Its value is found to be 

^O^'+vTO' 2 - Ofp' 2 - 0] cos 

Writing /i = cos 0, ;t' = cos O' and 

cos 0 = cos 0 cos O' + sin 0 sin O' cos tj> 

we obtain the result 

P„(cos 0) = P/cos 0)P„(cos O’) 

+ 2 t (_l)"t!lZ^p;( C o S 0)P;(cos0) ( cos(ffl^) (22.5) 

m = t (H + Wj! 

which is often of value in the solution of problems in wave 
mechanics. 

§ 23. Surface spherical harmonics. From the two sets 
of orthogonal functions 7"”(cos 0), cos wc can form 
a third set of functions 

X m , n (0, 4>) 

= {o7Tm)l}* 7 " ,(cos0)COS m4> ' ( '"- n) (23 '° 

which is an orthogonal set of functions on the unit sphere, 
i.e. the functions of the set satisfy the normalisation relation 

f sin OdO f X„ t m X„,_ m -d<j> = 5 nn .5 mm : (23.2) 

Jo Jo 

In a similar way we can construct a set 

(m£n) (23.3) 


+ For details, sec E. W. Hobson, op. eit., pp. 365-71, 
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which satisfies the relations 

I * sin OdO [ 2 ‘ y„, m y„, M .d<t> = 5 nn .5 am ., (23.4) 

Jo Jo 

psinOf/fl I** X„.„Y n .. m .d<l> = 0 (23.5) 

Jo Jo 

for all integral values of n, n', m and tri with m£/>, 

Because of these orthogonality relationships we can 
establish an expansion theorem which is a straightforward 
generalisation of the Legendre series (16.3). It is readily 
shown that for a large class of functions /, the function 
J{0, 4>) can be represented by the series 

£ COS0) 

n = 0 

+ 11 (23.6) 

b = 1 m = l 

where the coefficients c nf x nmi y„ m are given by the ex¬ 
pressions 

c - = ~r~ P d <t> f'M 4>)P n (cos 0 ) sin 0 do, (23.7) 
4k Jo Jo 

** X m (0, (23.8) 

} 

(23.9) 

0 

For any given function f(0 7 <p) the series (23.6) can therefore 
in principle be computed by a series of simple integrations. 

The functions X„ f m and Y„ t m? which are known as 
surface spherical harmonics, can be constructed easily from 
the known expressions for the associated functions T™. 
We find, for instance, that 

X K i(&* $) - — s * n ® cos ^ 
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i(«, «--(■ 

—\ sin 0 cos 0 cos 4>, 

\nj 

x 2 , 2 ( 0 ,4>) = (j 

—) sin 2 0 cos 2 <j>, 

I6ji / 

* 3 . .(0. *) = - (: 

—V sin 0(5 cos 2 0-1) cos 4>, 
\2n) 

x 3 . 2 ( 0 ,<t>)= (■ 

— ) sin 2 0 cos 0 cos 2<j>, 

6 nj 

1 

1 

II 

w* 

m 

n 

* 

— ) sin 3 0 cos 3^, 


and the corresponding expressions for the Y m> m are ob¬ 
tained by replacing cos (m^) by sin (m$) in the expressions 
for the X Hi 

The functions X Ht m and Y Ht m have the important 
property that they are solutions of the partial differential 
equation 


i a 

sin 0 SO 



+ 


1 d 2 X 
sin 2 0 d # 2 


-m(n+lpf = 0 


(23.10) 


so that the function 


(. AP'+Br-‘- l )X n> JLO, 0) + (Cr"+nr-"- 1 )y, 1 , m (0, <f>), 

where A, B, C and D are constants, is a solution of Laplace’s 
equation. It follows immediately from equations (23.7)- 
(23.9) taken with the expansion (23.6) that the function 

Hr, 0 , <M = c« Q ^( cos G ) 

+ i i WKJ,ju)+a,j.«} 

N - i m * I \aj 

G 
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satisfies Laplace’s equation in the region 0is finite 
at r = 0, and takes the value yffl, <j>) on the sphere r = a. 

For example* suppose we wish to find the solution of 
Laplace’s equation which takes on the value x 2 on the 
surface of the sphere r — a. Here we have 


/(0, 4 s ) — sin 2 0 cos 2 4> 
( 3 cos 2 0- 


a 2 a 2 ( 3cos 2 0-l\ ( 2 . 2 „ *, 

—--I- +i*r stir 0 cos 26 

3 3 V 2 J 

= ~ - j ? 2 (cos 0)+ X 2t M <t>)a\ 


Thus the required solution is 

4>(.r, 8, 0) = ^a 2 - 1 r 2 P 2 (cos 0)+ r 2 X 2i 2 {Q, <f>). 

Substituting the values of P 2 and X 2t 2 and transforming 
back to Cartesian coordinates we see that the required 
solution is 

ip = l(a 2 + 2x 2 —y 2 ~z 2 ). 


§ 24. Use of associated Legendre functions in wave 
mechanics. To illustrate the use of associated Legendre 
functions in wave mechanics* we shall consider one of 
the simplest problems in that subject—that of solving 
Schrodinger’s equation 

VV+ (W-VW = 0 (24.1) 

hr 

for the rotator with free axis, that is for a particle moving 
on the surface of a sphere. In equation (24.1), W represents 
the total energy of the system, V the potential energy. 
In the case under consideration V is a constant, V 0 say, 
and the wave function ip will be a function of 0, <$> only. 
If the radius of the sphere is denoted by a then equation 
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(24.1) is of the form 
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J_ dhp cotO # 1 d 2 ip 8n 2 m(W-V 0 ) 

a 2 30 2 a 2 39 a 2 $mQd(f> 2 h 2 

(24.2) 

If we consider solutions of the form 


41 = Qe ±im + 


then 


e-H-e'c* 

h 2 


m 


sin 


2 0 


0 = 0 . 


Substituting 

n %n 2 ma 2 (W-V 0 ) { 1A 

/t = cos 0, ---— = n(/i 4-1) (24.^) 

h 2 

we find that this equation reduces to Legendre’s associated 
equation (21,1) and hence has solution 

0 = 4P“(cos 0) + B(JJ(cos 0). 

However, for the same reason as in the case of potential 
theory (§ 20 above), wc must take B = 0. The solutions 
of equation (24.2) will therefore be made up of com* 
binattons of solutions of the form 

K M <P) = A„ n e ±tm *P"(cos 0 ), (24.4) 

where A mn is a constant. 

The physical conditions imposed on the wave function 
ip are that it should be single-valued and continuous. 
Obviously then the “ physical ” solutions will have m an 
integer, since \pm, $ + 2 ji) must equal ip mt n (9, $i* 

Further, in order that the series for shouid converge 
for the values fi = ±1 it is necessary that it should have 
only a finite number of terms. This is possible only if n 
is a positive integer. If therefore the solution (244) is to 
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be valid for 6 = 0 and 0 = ji we must have n a positive 
integer. The physical conditions on the wave function 
are therefore not satisfied by systems with an arbitrary 
value for the energy W but only by systems for which 

W=V 0 + -%—;nin+i), (24.5) 

o7t fttt 3 

where n is a positive integer. In other words, the energy 
of such a mechanical system does not vary continuously, 
but is capable of assuming values taken from the discrete 
set (24.5). 


Examples III 

(1) Show that, if n is odd P„(0) = 0, and that, if n is 
even, 

(—1) 4 "«! 

" (i«)! 2"{(*n)!}*‘ 

(2) Prove that 


as n n + 1 

I —- = i log 

« = o n +1 

(3) If C = p + J (,<*-]) show that 




(0(1-*£)-*(!-A/0"*“ I h n P„(tx); 

ft — 0 

(ii) PM= -n;i-n;r 2 ). 


«!f(i) 


Deduce that 


«ir(i) 


(hi) iFiih -n;i-n; 1) = 

r<n+*> 

(4) If n is a positive integer prove that 


I 


1 p a (p)(i-2nh+h*r%t= 

-1 2n +1 
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and hence, making use of Rodrigues' formula, deduce that 

fi 72 a+l / n| \2 

(1-P 2 )V -2pft + /i 2 )—*dp = “ 

J-i (2n + l)! 

(5) Prove that 

Pfc) = X (2n-4r-l)P„_ 2r _ l (x), 

r = 0 

where p — 1) or — 1 according as n is odd or even. 

Deduce that for all x in the closed interval (-1, 1) and 
for all positive integers n y the values of the functions 

| P&)\, n -2 | Ph(x) |, n~* | P”(x) |, ... 

can never exceed unity. 

(6) Prove that 


1' 


p„(mW 1 


2n + l 


tt.-i&O-J’.+ iW}, 


and deduce, from example 1, that if n is an odd integer 


I 




What is the value of the integral when n is even? 

(7) Using equation (14,2) and the results of the last 
example show that, if n is even, 


r 


(«— 2)1 


pP,(p)dp = (-i)** 1 - -—. 

o 2 (in + l)!(in —1)! 

and that the integral has the value 2 ero if n is odd. 


(8) If 


w» = 

prove that (n+\)u n + } 


PMPn-M ^ 

J-I P 

+/w rt = 2. Hence evaluate 


««■ 
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(9) If in and n arc positive integers, prove (hat 


1. 


1 . 2" + * +1 f('m + nVf 2 

a+nr* m wto = ~— t, " • 

_i m!(m+2ij + l)! 


(10) If n is even and m> — I, prove that 


r(im+i)r(im+i) 


i: 


U m PJjt)dll = 

J0 2r(i»i + in + l)r(ii?7-in+I) 

Deduce that 

*i 


I 


(l-kft 2 )-”-ip 2n (fi)d t i^ - 1 - 
o 2n + l 




(II) Show that 


and hence that 


It 

(}zj) = t i Q ^ 2,r - 


Deduce that 

P n (cosh 

(12) If f(ji) ~ (^ 2 -“l) B show, by using Rollers theorems 

that /'Qt) must have at least one zero between — I and 1. 
Proceeding in this way deduce that has n zero, 

between —1 and 1, 

Hence show that when n is even the zeros of P n (jt) 
occur in pairs, equal in magnitude but opposite in sign, 
and that when n is odd, /t = 0 is a zero and the others 
occur in equal and opposite pairs, 

(13) If >'(/0 is any solution of the linear differential 
equation 


dfr dii 
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ill which of, /f and y are continuous functions of /( whose 
derivatives of all orders are continuous, prove that y{p) 
cannot have any repeated zeros except possibly for values 
of /( which satisfy the equation a(j.t) — 0, 

Deduce that all the zeros of P„(/<) are distinct, 

(14) Prove that the Legendre polynomial P n (x) has the 
smallest distance in the mean from zero of all polynomials 
of degree n with leading coefficient 2 n (i)Jnl 

(15) If R denotes the operator 

prove that 

| -n(n +1) j* P„(x)f(x)dx 

provided that f(x) and f'(x) are finite at x = +1. 

Prove that if 1 


r 


log (1 —*)P „(x)dx = —-- 

-i ii(n + l) 

(16) Prove that 

= V 2 - 

2n + r 

P n (x)dx _ 2h n 


r 1 p„( X )dx = 

J -. va -*) 
p 

(ll) J _! (l-lhx + h 1 )* - 1 - ft 2 ' 

(17) If R 2 - l—2hx+h 2 , prove that if | h |< 1, I, 


(0 

(ii) 


[ log (1 - hx + R)P„(x)dx = 


2/t" * 1 

(2h + 1)(h + 1)’ 
2ft" 

h(2m +1) 
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(18) If 
/ 

prove that 


f(x, a) = --- — d%, (m > 1), 

J 0 0-2fx+£ 2 )* 

f 1 2a" +ffl 

/(x, a)P„(x)dx = --- --- 7, -f 

J-i (n + m)(2M+1) 


(19) If z is real and ( z J < 1 prove that 
dd> 7t 


r 


J 0 1 + Z COS 0 ,/(l -z 2 ) 

Putting z = +/i^/(/i 2 —1)/(1 -/yt) where A is so small that 

| h{fi± s /(ji i -\)cos <f>} |<I 

expanding both sides in powers of h , and 
equating coefficients of if show that 

PM = - f cos 4>}"d4>, 

n Jo 

Hence evaluate the sum 

t n C r P,(cos 0). 

r =* 0 

(20) Show by making the substitution 
Z = ±hJ(jt 2 -])l(hti- 1) 
in the formula in ex, 19 that 

d<(> 


P 00 = l p 

« Jo 


0< + V0* 2 -l)cos 


(21) Making use of the integral expression for F n (/0 
derived in ex. 19 show that 
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Deduce that 


PM 

PM 

PM 

= (l -/< 2 ) 3 

PM 

0 

PM 

PM 

PM 

PM 


0 

PM 

0 

P 2 (f 0 

PM 

PM 


Pit® 

0 

PM 


(22) Prove that if /j > 1 

and deduce that 

5*00 = | {H-JU* 1 ~ 1) cosh 0}" dO, 

where a = ± tog(/i + l)/(/i- 1). 

Hence find expressions for Q 0 (n) and Qi(fi). 

(23) Determine the simple expressions for Q 0 (n), 
QM< QAd) and Qs(ji) by working out the values of the 
polynomial H^ n _,(/j), occurring in equation (18.8), for these 
values of n. 


(24) Establish the following formula: due to MacRobcrt: 

iPJLMJiQ 


(0 pPMQJj*) 


< 




■ d%, m<n; 


l 


2n + l 


os) T ^ + 

J-i /*-£ 

(iii) Wn 1 -t)Q'M = i [ 1 d{, n>m; 

(iv) (» 2 -l)PMQ'M = l P dt, n>m. 

J-i n-Z 

(25) Prove that, if n is a positive integer and ( — /* + 

s /o< 2 -n, 

| _ r f * P„(cos fl) sin 6 dO 
'"'Jo I —2C cos 0+C 2 


QnU 0 
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Deduce that, if [ ( | > 1, 

ft.GO “ Z C" m j F„(cos 0) sin mO 
"* 1 Jo 

By evaluating this integral show that 


dd. 


- V 




2 F,(i, n+1; n-f-J; C 2 ). 


r(n+j) 

(26) Prove that, if m is a positive integer, 
f ft—T-fo) = (— l) w (2m) !(1 — 

nT m " 2" , m!(l~2/r/i + /i 2 ) m+ *' 


(27) Show that, if | /* [> I, n> — 1 then 

i= r(n+m + l) (ji 2 - 
r(n + l) 2 
Deduce that, if | /r+1 |>2 then 

Q"( /i ) = V ,ir ( n + m + 1 ) O' 2 -!) 1 


o"oo = ‘v—^ v* - 1 )*" f‘ 

; ~. "," +l J_, 0(-O" + ' n+l 


2"“T(» + 3) (/I + I> 


.H + J 7 I+ 1 


.Ft 


^n + 2, n 


+ iw+1; 2n + 2; 


ii+ij 


Find a simple expression for Q„ l (p). 

(28) Prove the following recurrence relations for 
Ferrer’s Associated Legendre Functions: 

(i) 7T+VCP) - IT-V00 = (2«+DO -M 2 ) 1 ™; 

(ii) (n - m + 1)7T + ,00 - (2n +1 )n1?(u) 

+ (n 4- w)7'„"L ,(/<) = 0; 

(iii) IT "00- TT-VW = (« - mXl “ 

(29) Derive the expressions for T7(0, <£) and X "(9, $) 
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for m ** i, 2, 3, 4. Express the functions sin 3 0 sin 2 
sin 2 0 cos 2 <j>, sin 0 cos 3 9 cos <f> in terms of surface spherical 
harmonics. 

(30) Find the function which satisfies Laplace’s equation 
in the interior of the sphere x 2 +y 2 +z 2 = a 2 , remains 
finite at the origin and takes the value ccx 2 +fiy 2 + yz 2 on 
the surface of the sphere. 

(31) The Jacobi polynomials are defined by the equation t 

PJfl, b, x) = jFjC-m, a+rn; b; x). 

If D = d/dx show that 

D"{x l ‘+ m - 1 (l-x)' +m - b }=(b) m x b - l 2 F I (b-a-m, b + m\ b; x) 
Deduce that 


(i) ejfi, b, x) = ZSJLJLlirtf*"- i(i 
(b)„ 

(ii) | ‘ x b -\l-xy- b f(x)Sr m (a, b, x)dx 

= — P (-l)7 ( ' ,,) (.x).x 6+m_1 ( 1 ~x) a+m -Mx; 
(b) m Jo 

(iii) j* x fc-1 (l -x)"'‘# r ffl (fl, b, x)F„(a, b, x)dx = 0, m + n 

(iv) | 1 x i -*(l-x)“- fc {^ w ( < i, b,x)} 2 

__ r(b)r(a + l-b)(a + l-b) m „ ml 

r(a)(fl) m (fe) m a+2m 

(32) The Tchebkhef polynomials of the first and second 

t It should be observed that this name is sometimes applied to 
the polynomial 

if" w = **$**&+p+i. «+i, »-wl 
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kinds may be defined by the equations 

r a (*) = 2 F, 

!/.(*) = (n+l) 2 F,(-H ) n + l;}; ^ 

respectively. 

From (2) on p. 43 show that 

T b (cos 0) = cos nO, (7„(cos 0) = S1 

sin 0 

Prove also that 


(i) 1+2^ T„(x)f = (1-I*X1-2l*+l*r\ 

| * | < 1, |x|<l; 

(ii) 2"(l),T t W = (- 1)"(1 -x 2 )'” f-(l-x 2 y-i; 

dx " 

(ii*) j 1 T m (x)T n (x)(l-x 2 y*dx^{nd m „- 


(iv) £ 0-2x1 +t 7 r , i 

It i* 0 

(V) 2" +1 (i) n + 1 t/„W 

=(-iY(n+w-x 2 r*f-(t-x 2 r +x -, 

dx n 

( vi ) J VmWU'ixXl-x^dx = ind^; 

(vii) T„(x) = V„(x)-xU„-i(x); 

(viii) {l-x 2 )t/ B _,W = xT n (x)~T n+1 (x), 

(33) The function C*(x) defined by 

(l-2xt+t 2 )- y = £ C„ v (x)r B (v>0) 

n 0 

is called the Gegenbauer polynomial of degree n and order v. 
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Prove the relations: 


(i)q<*) = 


f(n + 2v) _ 
h!F(2v) 2 1 


( 


(u)c:(i) = (-i)"c;(-i) = 


n + 2v, -it; v+i; 

r(n + 2v). 
n!F(2v) ’ 



(iii)C;(x) = (-2)" 


r(v+n)r(2v + n) 
n!r(v)r(2v+2n) 


(l-x 2 )* _v 


dx” 


(1—x 2 )" +v- *; 


(iv) p c:(x)c:(x)(i-x 2 y-Mx 


2'- 2v nr(M + 2v)^ 
n!(v+n)[r(v)] 2 ~ 














CHAPTER SV 


BESSEL FUNCTIONS 

§ 25, The origin of Bessel functions. Bessel functions 
were first introduced by Bessel, in 1824, in the discussion 
of a problem in dynamical astronomy, which may be 
described as follows. If P is a planet moving in an ellipse 
whose focus S is the sun and whose centre and major axis 
are C and A'A respectively (cf. Fig. 8), then the angle 



Fig. 8 

ASP is called the true anomaly of the planet. It is found 
that, in astronomical calculations, the true anomaly is 
not a very convenient angle with which to deal. Instead we 
use the mean anomaly £, which is defined to be 2n times 
the ratio of the area of the elliptic sector ASP to the 
area of the ellipse. Another angle of significance is the 
eccentric anomaly, u, of the planet defined to be the angle 
ACQ where Q is the point in which the ordinate through 
P meets the auxiliary circle of the ellipse. 


S2S 
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It is readily shown f by a simple geometrical argument 
that, if e is the eccentricity of the ellipse, the relation 
between the mean anomaly and the eccentric anomaly is 

C=h — esinu. (25.1) 

The problem set by Bessel was that of expressing the 
difference between the mean and eccentric anomalies, 
h—C, as a series of sines of multiples of the mean anomaly, 
i.e. that of determining the coefficients c t (r = I, 2, 3, ..,) 
such that 

«-C= £ c r sin (tO- (25.2) 

t - i 

To obtain the values of the coefficients c r we multiply 
both sides of equation (25.2) by sin (jf) and integrate 
with respect to ( from 0 to 7i. We then obtain 


£ 


(«-0 sin (sQdC 


£ f* • 

= I Cr S' 
r “ 1 Jo 


sin (r£) sin (s0d£. 


Now 


f 


sin (rQ sin (s()^C = 


and an integration by parts shows that 
(»-0sin (sQd( 


1 

= £ (f-«) COS (sC)j + “J - l) COS (sCK. 

From (25.2), £—« is zero when £ = 0 and when £ — ti, 
so that the square bracket vanishes; the integral can be 
written in the form 


ip 

* Jo 


cos s£ du 


t Cf* D. E* Rutherford, Classical Mechanics (Oliver and Boyd, 
1951), § 42. 
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and hence, using equation (25.1) we obtain the result 

c s ——\ cos{s(n-esin u)}dw. (25.3) 
Its Jo 

The integral on the right-hand side of equation (25.3) is 
a function of s and of the eccentricity e of the planet’s 
orbit. If we write 

j n ( x ) = i f cos(x sin d—nff)d0 (25.4) 

*Jo 

it follows from equations (25.3) and (25.2) that 

ti-C = 2 £ J,(er) S ^^-. (25.5) 

r - 1 r 

The function J„(x) so defined is called Bessel’s coefficient 
of order n. 

We shall now show that 7„(x) is equal to the coefficient 
of t” in the expansion of exp ')}; in other words 

we may define J„(x) by means of the expansion 

exp jix (t - ^ t * JJtxY. (25.6) 

To prove this we need only show that the J„(x ) of (25.6) 
can be expressed in the form (25.4). We first of all observe 
that 

£ (-l)V-j(x)f= £ J„(x) (-;)"= £ Ux)t" 

!,= -«, n “ —to \ X} o = -® 

since both expansions are equal to exp {i*(f — 1/f)}. 
Equating coefficients of t " we have 

(-l)V_ n (x) = J B (x). (25.7) 

In the expansion (25.6) we may write t = e' e to obtain 
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the relation 

00 

exp (ix sin 0) — £ J„(x)e inB . 

It “ - 00 

Making use of the result (25.7) we see that the series on 
the right can be put into the form 

JoM+2 £ J 2m (x) cos (2m0)+2i £ J 2|B+I (x)sin(2m+1)0 

m ■ 1 in = L 

so that by equating real and imaginary parts we obtain 
the expansions 

cos (x sin 0) = J 0 (x) + 2 £ / 2m (x) cos (2m0), (25.8) 

ffl “ 1 

sin (x sin 0) = 2 £ J 2m+1 (x) sin (2m+ 1)0. (25.9) 

ffl - 0 

If we now multiply (25.8) by cos nfl, (25.9) by sin nO, 
integrate with respect to 0 from 0 to it, and use the formula 

J* cos (mO) cos (nfl) dO = J sin (m0) sin (nfl) dO = ^ „ 

we obtain the formula: 

7„(x) = - | cos (x sin 0) cos (n0) d9, (« even), (25.10) 
it Jo 

J„(x) = — f sin (x sin 0) sin (nfl) dO, (n odd). (25.11) 
it Jo 

Because of the periodic properties of the trigonometric 
functions we know that the integral on the right of equation 

(25.10) is zero if tt is odd, while that on the right of equation 

(25.11) is zero if n is even. Thus for all integral values of n, 
we have 

J n (x) = - f {cos (x sin 0) cos (nfl)+sin (x sin 0) sin (nO)}dO 
it Jo 

which is identical with the expression (25.4). 

H 
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la particular 

.JoO) = - f’ COS (x sin 0 ) do. (25.12) 

* Jo 

In what follows we shall assume that the Bessel co¬ 
efficients of the first kind are defined by equation (25.6) 
or, which is equivalent, by equation (25.4). 

§ 26. Recurrence relations for the Bessei coefficients. 
If we differentiate the generating equation (25.6) with 
respect to x we obtain the relation 

which is equivalent to 

i W = 0. 

Z It “ — CO I) = —00 

Equating to zero the coefficient of t" we obtain the relation 

2j;(x) = J fl _,(x)-J fl+1 (x). (26.1) 

On the other hand, if we differentiate (25.6) with respect 
to i the resulting equation is 

ix (l + i) exp {ix(l- i)} - _ £ _ 

and this is equivalent to the relation 

ix £ (l»+l"- 2 M.(x)- f nJ£x)f~ l = 0. 

it =* -oo n o ~a> 

Equating the coefficient of t n_i to zero we obtain the 
recurrence relation 

2 Jij n (x) = J n - l (x)+J„ +l (x). 

X 


(26.2) 
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Adding equations (26.1) and (26,2) we find that 

xJ' n {x) = xJ„. i(x) - nJ„(x) (26.3) 

and subtracting equation (26.1) from (26.2) we obtain 

xJ'„(x) = nJ n (x) - xJ„+ i(x). (26.4) 

Putting n = 0 in this last equation we have the important 
special case 

J'oix) = (26.5) 

and putting n = 1 in equation (26.3) we find that 

= (26.6) 

X 

Differentiating both sides of (26.5) with respect to x and 
making use of the result (26.6) we have 

/SC*)=-Jo(*)+~J r i(x) 

X 

which, as a consequence of equation (26.6), may be written 


Jo{x) + - J' 0 (x)+J 0 (x) = 0 (26.7) 

X 

showing that y = J 0 (;r) is a solution of the differential 
equation 




(26.8) 


We can show similarly that the Bessel function J„(x ) 
satisfies the differential equation 

0 + ;l + (‘-?) J,=o ( " in,egril) ' <26 - 9) 

For, from equation (26.4) we find, as a result of differ¬ 
entiating both sides with respect to x, that 

xr H (x)+J'„(x) = nJ' n (x)-J n+ ,(x)~xj' + ,0c). 
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Now, from equation (26,4) t 

nJ'„(x) = n — J„(x)-nJ n+1 (x), 


and putting rt + 1 in place of n in equation (26.3) we see 
that 

xj; +i (jc) + (n + l )J„ + i(jc) = xJ n (x) 

so that 

xJ“(x)+J'„(x) = — J n (x)-xJ„(x) 

X 


which shows that J n (x) is a solution of equation (26.9) 
provided, of course, that n is an integer. 

As we pointed out in g 1, equation (26.9) is known as 
Bessel’s equation. What we have shown is that if the n 
which occurs in Bessel's equation is an integer, one solution 
of the equation is J R (x)> It is because of this fact that 
Bessel coefficients are of such importance in mathematical 
physics, for as we saw in § 1, the equation (26,9) arises 
naturally in boundary value problems in mathematical 
physics* 

g 27* Series expansion for the Bessel coefficients. We 
shall now find die power series expansion for the Bessel 
coefficient If we write 


exp 




and make use of the power series for the exponential 
function we obtain the expansion 


§ 27 
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By our definition (25.6), the Bessel coefficient J„(x) is the 
coefficient of in this expansion. If n is zero or a positive 
integer, we find that 


J*(x) = 


CO 


£ 


(-1) 1 (x\ n+ls 
s!(n+s)! \2/ 


(27.2) 


and when n is a negative integer we can deduce the series 
for ./„(*) from equation (25.7). 

Writing equation (27.2) in the form 


we see that 


z - (—lx 2 )* 

2 nt 1 = 0 s!(n +1), 


^W= + ~ix 2 ). (27.3) 

z n i 


The variation of the Bessel coefficients ./ 0 (x), J t {x), 
J 2 (x) for 0£xg20 is shown graphically in fig. 9. These 
are the Bessel coefficients which occur most frequently in 
physical problems and their behaviour is similar to that 
of the general coefficient 

Simple relations for the Bessel coefficients may be 
derived easily from the series expansion (27.2). For 
example, since this equation is equivalent to 


X M J*(X) = 


CO 



(—l)V "* 21 
s!(n+s)l 



(27.4) 


it follows, as a result of differentiating both sides of this 
equation with respect to x, and making use of the fact 
that (2n+2s)l(n+s)\ = 2/{n-l+j)!, that 


!;{".«} = 



+ /|Vi-l + 2s 
s!(»i —1+s)! U.J 


which, by comparison with (27.4), shows that 


-{.x"J„(.v)}=A-,«. 


(27.5) 
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If we write this result in the form 

x dx 



we see that if m is a'positive integer less than n, then 

(I i-Y xV n (x) = x- m J n - m {x). (27,6) 

\x dxj 

Similarly we can establish that 

i-{*-"/„(*)} = -x-"J B+1 (*) 
dx 

or, which is the same thing, 


(27.7) 
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a result which may be generalised to the form 

(x £)" <*~ v " ( * )} = (-ir* - — j^x). 

In particular we have the relation 

x-Ux) = {-1)" £j Jo(x), (27.8) 

which shows how the Bessel coefficients J n (x) may be 
derived from J 0 (x). 

Another interesting property of the Bessel coefficients 
also follows from the power series expansion (27.3)* This 
concerns their behaviour for small values of the argument x. 
Since 

Uni qFjCm + I, ~ix 2 )= 1 

je-*0 

it follows from equation (27.3) that 

lira = -L. (27.9) 

x-*Q 2 HI 

In other words, for small values of x, the Bessel coefficient 
J n (x) behaves like x f, /2 fl /i!. 

§ 28. Integral expressions for the Bessel coefficients. 
We have already derived one integral expression for the 
Bessel coefficient of order n (equation (25.4) above). In 
this section we shall consider other simple integral ex- 
pressions for these coefficients. 

We shall consider the integral 

/ = j‘ (i -fy-Wdt 

in which n> — 4. If we develop exp (ix/) in ascending 
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powers of ixt we see that the value of this integral is 

£ (l-py-ifdt. 

If s is an odd integer then the corresponding integral 
occurring in this series is zero, and if s is an even integer, 
2 r say, then the integral has the value 

r - tSLtlSM 

Jo r(n + r+l) 


so that 


i= t 


(-I)V' r( n +i)r(r+i) 


r» o (2r)! 

= r(i)r(«+i) t 


r(n + r+l) 

(-l)V 


r = o r!r()i + r+l)2 2r 

since, by the duplication formula for the gamma function, 
r(-£)(2r)! = 2 2 V!r(r+i). It follows immediately from the 
series expansion (27.2) for J„(x) that 


JnW = 


dxy 




t 2 )"~ i e u, di (28.1) 


iwtiH-i), 

and it is easily shown that this is equivalent to the formula 
x" 


J»(x) 2 "- 1 r(i)r(n+i) 

In particular 

Ja(x) 




cos (xt)dt. (28.2) 


- r 

* Jo 


cos (xi) 

Vo-' 2 ) 


dt. 


(28.3) 


The result (28.2) may be expressed in a slightly different 
form by means of a simple change of variable. It we put 
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t - cos 0 we obtain the integral expression 

J„(x)= - ; — - f cos(x cos 0) $in 2n 0d6, (28.4) 

2 n - 1 r(i)r(n+i) Jo 

while if we make the substitution / = sin 0 we get the 
formula 

J„(x ) =-—- 1 1 " cos (x sin 0) cos 2n 0d0. (28.5) 

2"-T(i)r(n+i)J 0 

The particular forms appropriate to n = 0 are 

2 f** 2 

j 0 (x) - - cos (x cos 0)d0 — - \ cos (x sin 0)dO * (28.6) 

Jo 31 Jo 

§ 29, The addition formula for the Bessel coefficients. 

In certain physical problems we have to reduce a Bessel 
coefficient of type J n (x+y) to a form more amenable to 
computation. We shall now derive an addition formula 
which is of great use in these circumstances. From the 
definition (25.6) we have the expansion 

= £ Jxix+y) 1 ”- 

Writing the left-hand side as a product 

and inserting the appropriate series from (25.6) we find 
that 

£ J'(x+y)t*= £ £ JXx)jjb>y + \ 

n — - CO r='-oo»w—co 

Equating coefficients of f we obtain the addition formula 

JXx+y)= £ J r (x)J n - r {y). (29.1) 
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To put this in a form which involves only Bessel co¬ 
efficients of positive order we write the right-hand side in 
the form 


-i 




£ Jr(x)J n -r(y) + 
- 0 


£ JAWn-Ay) 

- n+l 


and note that because of the relation (25.7) the first term 
can be written as 


£ 


— co 


(-ly^-rWn-rO) 


£ (-i yj r (x)J B+ Ay). 

r = 1 


Similarly the third term is equal to 


£ J n +r(x)J-Ay)= £ (-l)V n+ ,(x)J,(>') 

r = 1 t - 1 

so that finally we have 


J*(x+y) = £ J r (x)J n - r (y) 

f- 0 

+ £ (-l)UW^ +r W+^ + ^K0')}. (29.2) 

r =1 


§ 30* Bessel’s differential equation* We showed pre¬ 
viously (§ 26 above) that* if n is an integer, J n (x) is a 
solution of Bessel’s equation (26.9). We shall now examine 
the solutions of that equation when the parameter n is 
not necessarily an integer. To emphasise that this para¬ 
meter is, in general, non-integral, we shall replace it by 
the symbol v, so that we now consider the solutions of 
the second order linear differential equation 




Writing the equation in the form 


— +(-v 2 +x 2 )v = 0, 
ax dx 


(30.2) 


no 
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we see that the point x = 0 is a regular singular point and 
that in the notation of § 3, p 0 = I and q 0 = — v 2 * The 
indicial equation (cf. (3.5) above) is therefore 

g 2 — v 2 == 0 

and this has roots q = ± v. 

First of alt we shall suppose that v is neither zero nor 
an integer. Then the first solution is of the form 

y= £ <V* P+y - (303) 

r = O 

Substituting this series in equation (30.2) we see that the 
coefficients c r must be such that 

£ {(v + r)(v + r — 1)+(v + r) — v 2 }c r x r+v -f £ cX +y * 2 - 0. 
t = o r tr o 

Hence we must have 

c,{(v+ l) 2 ~ v2 } = 0 

and in general 

e r {(r + v) 2 -v 2 } = -c,_ 2 , (r = 2, 3, ...). (30.4) 

We must therefore take c x to be zero and hence, in order 
that (30.4) may be satisfied for all r^2, we must take 

C 2r+ i=0 

and 

c = _ (-Ifro _, 

(2r+2v)(2r+2v - 2),, .(2v+2)2r{2r - 2).. .2* 
an expression which may be put in the form 

c lf =-M-AY. 

r!(v+l),\ Aj 

Taking c 0 = l/2 v v! we see that the basic solution of type 
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(30.3) may be taken as 


** f (-ix a r 

2T(v+l) r - o r!(v + l)/ 


(30.5) 


Comparing this series with the series (27.2) we see that it 
is of precisely the same form as that equation, the only 
difference being that n is replaced here by v. If we take 
the series (27,2) to define the Bessel function of the first 
kind of order n t even when n is not an integer, then we 
may write the solution (30.5) in the form 


y = A(*). 

Similarly, if we substitute a series of type 


y= £ * 

r - 0 

to correspond to the second root of the indicia! equation, 
we find that it must be of the type 


y = 


f (~i* 2 ) r 

2T(—v+l) r = o r!(—v+l) r 


(30.6) 


and with the extension of the definition (27.2) to non- 
integral values of v we may write this solution in the form 

y = J-v(a-). 

Thus when v is not an integer we may write the general 
solution of equation (30J) in the form 


y = AJAx)+BJ„,(x) (30.7) 

where J v {x) is defined by the equation 


j (jc) = _JC- 0 F,(v+ 1 i -|x z ). (30.8) 

2T(v+l) 

It should be observed that the results of §§ 27, 28, with 
the exception of (27.8), are true when n is not an integer. 


5 30 
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since they were derived directly from the definition (27.2), 
which is equivalent to (30.8). The transition is effected 
merely by replacing n\ by T(v+ 1 ). 

When v is zero or an integer we know from equation 
(25.7) that the solutions J*(*) and J- V (x) are not linearly 
independent. We must therefore use the formulae (3,8) to 
calculate the second solution. 

We shall: consider first the case in which v = 0. If we let 


w = c r x F+g 

t « 0 


then in order to satisfy the recurrence relation (30.4) we 
must have 

w = x t £ izMl (30.9) 

r-0Hfe+l) r 

and putting q = 0 we obtain the first solution 

Wq = Jo(x). (30.10) 

Using the result 

l-l -!_(£ JJ 

Oete+i), (e+l), V-ie+sJ 

we see that 

^ =wl 0 g x-x« £ s=MLUsl 

dQ e r -irl(e+l)r\'-iQ+4 

Putting 0 — 0 and substituting the value (30.10) for w Q 
we find that the second solution (dwld(?\ = 0 is 

Yo00 = Jo(x) logx- £ t(30.11) 
rml ( r! ) 

where 

<P(r)= £ (30.12) 

The function Y 0 (;c) so obtained is called Neumann’s 
Bessel function of the second kind of zero order. Obviously 
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if we add to Y 0 (x) a function which is a constant multiple 
of J 0 (x) the resulting function is also a solution of the 
differential equation 


d ll + + v = o. 

</.x 2 x dx 


(30.13) 


In particular the function 

Y 0 (x) = - {Y 0 (x)-(log 2-y)J 0 (x)}, 
% 


where y denotes Eulers constant, will be a second solution 
of the equation. Substituting from equation (30.11) for 
Y 0 (jc) in this equation we obtain the expression 

Y 0 (*)=-{log(ix)+y}J 0 (x)-^ t 4>{r) (30.14) 

% n r-i (r\y 

where <j)(r) is defined by equation (30.12). 

The function Y 0 (x) so defined is known as Weber’s 
Bessel function of the second kind of zero order. 

Thus the complete solution of the equation (30.13) is 

y = AJ Q (x) + BY Q (x) (30.15) 

where A % B are arbitrary constants and J 0 (*X Y 0 (x) are 
given by equations (30.8) and (30.14) respectively. 

It can be shown by an exactly similar process that when 
v is an integer the complete solution of the equation 
(30.1) is 

y = AJ v (x)+BY£x) (30.16) 


where A> B are arbitrary constants, J v (x) is defined by 
equation (30.8) and Y v (x) is given by 


Y„(x) = - {v + log(ix)}J v (x)- - V f 1)1 
n n r «o r! 




(~l)W +2f 

r!(v+r)! 


Wr+vnmb 


(30.17) 
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§30 

The function K v (x) so defined f reduces to the Y 0 (x) of 
equation (30.14) as v-^0 and is known as Weber’s Bessel 
function of the second kind of order t?. 

The variation of Y 0 (x) and Y j(x) for a range of values 
of x is shown graphically in fig. 10. 

The functions J v (x) and Y v (x) are independent solutions 
of the equation (30.1), but in certain circumstances it is 



advantageous to define, in terms of them, two new in¬ 
dependent solutions. If we write 

Hl'Xx) = /,(x) + /Y v (x), (30.18) 

H?\x) = J v (x)-iY v (x), (30.19) 

then it Is obvious that we can take the general solution 
of Bessel’s differential equation (30.1) to be 

y rn A^XxHAi ff‘ 2 >(x), (30.20) 

where A t and A 2 are arbitrary constants. The functions 
Hj°(x), Hj 2) (x) defined by equations (30.18) and (30.19) 
t This function is denoted as M v (a) by Courant and Hilbert. 
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are called Hankers Bessel functions of the third kind of 
order v. 

The Hankel functions and H {2 \x) bear the same 

relation to the Bessel functions JXx), i\(x) as the functions 
exp(±ivx) bear to cos vx and sin vx, and are used in 
analysis for similar reasons. It should also be noted that 
the Bessel functions Y,(x), H[ l \x), H {2 \x) satisfy the 
same differential equations and recurrence relations as the 
function J„(x). 


§ 31. Spherical Bessel functions. A problem which 
arises in mathematical physics is that of the solution of 
the wave equation in spherical polar coordinates 


aV + 23£ + 

dr 2 r dr 


1 


r 1 sin 0 


- (sin fl 
c0 \ dOj 


. 1 
r 2 sin 2 0 d<l> 2 


c 2 ax 2 ' 


(31.1) 


I f we take a solution of this equation of the form 

* = Y*. n(0> mr)e le “, (31.2) 


where Y m _ „(0, tj>) is the surface spherical harmonic defined 
by equation (23.3) and >p(r) is a function of /• alone which 
satisfies the equation 


dr 2 r dr r 2 c 2 


(31.3) 


Now putting 


4> = r~*R 


we see that equation (31.3) becomes 


(31.4) 
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whose general solution is readily seen to be 

R = AJ„ +i ((or/c)+BJ _„_ i ((orlc). (3 l.S) 

Hence the function 


* = r~*Ym, *(0> Mif'+vivrlcW* 1 ( 31 . 6 ) 
is a solution of the equation (31.1). 

The functions J ±(n+i) (x) which occur in the solution 
(31.6) are called spherical Bessel functions. We shall now 
show that they are related simply to the circular functions. 
First of all we consider the Bessel function J+(x). If we 
let v = | in equation (30.6) and make use of the duplication 
formula for the gamma function we obtain the result 


Jj(x) = 


_2 V S (~l) f x 2r+1 
nx) r = o (2r + l)l 


which shows that 


J^x) = 



(31.7) 


Again, if we put v => in equation (30.6) we obtain 
the relation 

./_$(x) = J cos x. (31.8) 


The other functions J m (x) where m is half an odd integer 
may be worked out in a similar fashion. It is left as an 
exercise to the reader to show that 


Jm(x) ■■ 


= ( ~) {/■(*) Sin X ~9m(x) COS x}, 
J-Jx) = (^J sin x+f m {x) cos x). 


where the functions f my g m 

i 


are given in Table I. 
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These functions which arise in the way described have 
been tabulated in Tables of Spherical Bessel Functions 
2 vols. (Columbia Univ. Press, 1947) prepared by the 
Mathematical Tables Project of the National Bureau of 
Standards. 


Tabic 1 


m 

fm 

Q m 

3 

l 

i 

2 

X 

5 


3^ 

2 

X 

7 

15 6 

^-1 

X 2 

2 

x 3 a: 

9 

2 

105 45 , , 

4 2 

r x 

105 10 

x 3 x 

U 

2 

945 420 15 

X s x 3 + x 

945 105 . , 

4 2 

X X * 


§ 32. Integrals involving Bessel functions. In this 
section we shall derive the values of some integrals in¬ 
volving Bessel functions which arise in practical applica¬ 
tions. In the first instance we shall consider definite 
integrals. 

From equation (27.5) we have the relation 


| 32 
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If n> 0, ^J n (x)^() as so that the lower limit is zero 
and we obtain the integral 


j: 


x n J„ „ jOHx = a n J „(a), (n > 0), 


(32.1) 


which, by a simple change of variable, gives the result 


j; 


= (n> 0 ). 


(312) 


A particular case of this result which is of frequent use in 
mathematical physics is obtained by putting /i — 1 in 
equation (32.2). In this way we obtain the integral 


I, 


rJ 0 (^r)dr=^J i (aa 
o Z 


(32.3) 


Further results may be obtained from (32.2) by familiar 
devices such as integration by parts. For example, making 
use of equation (27.5) we may write 

I*' r 3 J 0 (Zr)dr = fr 2 |J {rJ t (Z>1}dr, 

Jo Jo Zdr 

and, integrating by parts, we see that the right-hand side 
of this equation becomes 




-tj> 


(Sr)dr 


which reduces, by virtue of (32.2), to 


a 3 , „ , la 2 


Z <r 

Now by the recurrence relation (26.2) we have the expression 


j 2 m= f j,«<i)-j #({«) 
& 
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so that finally we have the result 

J"r 3 J 0 «r)dr=^|j 0 (^)+ (32.4) 

Combining this result with equation (32.3) we obtain the 
integral 

\‘ r(a 2 -r 2 )J 0 (Zr)dr = ^ Jtfa)- ~ J 0 (£a). (32.5) 

Jo S b 


The most commonly occurring infinite integrals are 
most easily evaluated by means of substituting the formula 
(27.3) in parts (ii) and (iii) of example 17 of Chapter II. 
From part (ii) of that example we see that 


r(ja+v-t-l)q T 

2T(v+l)p" +v+1 

zFi (i^+iv+i, i^+Jv + 1; v+1; - (32.6) 


If we make use of equation (30.8) on the left-hand side of 
this equation and of equation (7.4) on the right-hand side, 
we see that this result is equivalent to the formula 


r 


J v (ax)^ ,l e px dx = 


r(p+v+ l)a’ 


2T(v + l)(« 2 -lV)*' t+ * v+ * 


x zF, i/t+iv + i; i v~ip; v+1; 


a 2 + p 2 


} 


(32.7) 


where p>0, /i+v>0. 

The hypcrgeometric series occurring on the right-hand of 
this equation assumes a particularly simple form if either 
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ft = v or (i — v +1, and we obtain the formulae 

a' 

r(i) V+p 2 r+’ 

2 ,+ T(v+i) pa' 


1 


I, 


J t (ax)x v e~ px dx = 2 I 3 V . + ^ - 
J r (flx)A; v+ 1 e~ px dx — 


m '{a 2 +p 2 y+*‘ 


12S 


(32.8) 


(32.9) 


Two special cases of the formula (32.8) which occur 
frequently are 


J 0 (ax)e px dx = —-i— 
o s/(a 2 + p 2 ) 


xJt(ax)e px dx — —--—. 

o (« 2 + p 2 )* 


(32.10) 


(32.11) 


Integrating both sides of equation (32.11) with respect to 
p from p to oo we find that 


I 


Jy(ax)e px dx = -- lArr~\ 

o a a^f(a £ + p*) 

A special case of (32.9) which is often needed is 


1 


xJ 0 (ax)e px dx = - j -f 
0 (n 2 + pT 


(32.12) 


(32.13) 


If we let p lend to zero on both sides of equation (32.7) 
we find that we can sum the hypergeometric series by 
Gauss's theorem (7.2) provided that ||<| v+1 |. We 
then have the result 


I, 


= yni + ^ + H 

0 a »+*T(i~ip+l V ) 


(32.14) 


Similarly from part (iii) of example 17 of Chapter II 
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we have the equation 


i; 


o^ifv+l; -ifl 1 3c a >"' w ^ +, - , rf* 


_ r(^+jv) 

2f 

which, because of (27.3), is equivalent to 






= oT(^ + jv) 

2 v+ V + T(v + 1) 


V +1 ^1 + tr-/.. , ^ + V +1 


■-$) 


(32.15) 


From parts (i) and (ii) of example 11 of Chapter II we have 
the special cases 

P« ..v„-i! s / 4 p 1 

j o x** 1 J£ax)e~ r ** I dx = T , (32.16) 

and 

r w-*. - ^4^ (»+> - $ 

(32.17) 

of which the most frequently used are 

f" .vJ 0 (ax)e-^dx = ie-W (32.18) 
Jo 2/ri 


and 


[ x 3 J 0 (ax)e p ’ x \ix = -L e -° 2 i+p\ (32.19) 


§ 33. The modified Bessel functions* By an argument 
similar to that employed in § 1 we can readily show that 
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Laplace's equation in cylindrical coordinates 

dfy i #, i ay ay _ 0 

3g 2 g e 2 d<p 2 dz 2 


possesses solutions of the form 

ip = e ±lv4,±im: R(Q) 


where R(g) satisfies the ordinary differential equation 


d^R + IdR 
dQ 2 e dQ 



(33.1) 


Writing x in place of mg we see that this equation is 
equivalent to the equation 


d 2 R ] 1 dR 
dx 2 x dx 



(33.2) 


If we proceed in exactly the same way as in § 30 we 
can show that if v is neither zero nor an integer the solution 
of this equation is 

R = AI v (x)+BI. t (x ) (33.3) 

where A and B are arbitrary constants and the function 
l v (x) is defined by the equation 

t Jx) = *' f 

A ’ 2T(v+1) r = o r!(v+l), 

= 2^M) oFl(V + l;iJf2) ' (334) 

Comparing equation (33*4) with equation (30*8) we see that 

/ v (x) = r'Uix) (33.5) 

a result which might have been conjectured from the 
differentia! equation itself. 
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If v is an integer, n say, then /_ n (x) is a multiple of 
/„(*) so that the solution (33,3) in effect contains only 
one arbitrary constant. By a process similar to that 
outlined in § 30 we can show that in these circumstances 
the general solution of equation (34.2) is 

R = AI n (x)+BK„(x) (33.6) 

where the function KJ(x) is defined by the equation 

W = (-ir +l l n (x){log(ix) + y} 

+i i i-m-'-w (M—* 

r - 0 rl 

+ K-D" Z 1 {4>(r)+4>(n+r)}(ixr +2r . (33.7) 

r = j r!(« + r)! 

The functions /„(x), K„(x) defined by equations (34.4) 
and (34.7) respectively arc known as modified Bessel 
functions of the first and second kinds. 

The result (33.5) is very useful for deducing properties 
of the modified Bessel function I„(x) from those of the 
Bessel function J n (x). For instance, when n is an integer 
it follows from equation (25.7) that 

/_„(x) = I„(x) (33.8) 

and from equations (26.1) to (26.5) respectively that 


2/;w = v.w+/.hW, 

(33.9) 

2 -/„(x) = / n _ I (x)-J n+1 (x), 

X 

(33.10) 

x/;(x) = x/ H _ 1 (x)-n/„(x), 

(33.11) 

xf'(x) = n/„(x)+x/„ + ,(x), 

(33.12) 

l’ 0 (x) = I Ax). 

(33.13) 
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§ 34 

Similarly equations (27.5) and (27.7) imply the relations 

J- (33.14) 

ax 

~{x-”I l ,(x)} = x-I n+l (x). (33.15) 

ax 

All of these relations can* of course, be derived directly 
from the definition (33,4) or I v (x) and it is suggested as an 
exercise to the reader to derive them in this way. 

It should also be observed that K„(x) satisfies the same 
recurrence relations as I n (x). 

The variation of I 0 (x\ I\(x) and J 2 (x) with x is shown 
graphically in fig. 11 and that of K 0 (x\ K v (x) and K 2 (x) 
is shown in fig. 12. 


§ 34. The Ber and Bei functions. If we wish to find 
solutions of the form 


* - 

of the diffusion equation 

+ 1 # _ [ # 

Bq 2 q Bq k dt 

we have to solve the ordinary differential equation 
d z R + 

dg 2 q dg k 

On changing the independent variable to ,v - we 

see that this latter equation is equivalent to the equation 


d 2 R 1 dR 
dx 2 x dx 


-iR — 0. 


(34 J) 


Formally vve may take the independent solutions of 













1-0 






Fig. 12. Variation of e*Ke(x\ ^K^x) and e‘K t <ix) with x. 
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§34 


this equation to be / 9 (f*je) and K 0 (i-x). Kelvin introduced 
two new functions bcr(x) and bei(.v) which are respectively 
the real and imaginary parts of /o(/*x)» > e - 

ber (x) + i bci (x) = / 0 (i*x). (34.2) 


From the definition (33.4) of / 0 (.v) we see that 


ber(x) = 


f (-1M 1 

o (2s!) J * 


(34.3) 



and that 


bei(x)= £ 

s = 0 


(2s+ 1) I 2 


(34.4) 


The variation of the functions ber(;t) and bei(jc) with x 
is shown diagrammatically in fig. 13. 

In a similar way the functions kerf*} and kei(jc) are 
defined to be respectively the real and imaginary parts of 
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the complex function K^x), i.c. 


ker (x) + i kei (x) = A' 0 (f*x). (34.5) 

From the definition (33.7) of A^(.v) we can readily show 
that 

ker(x) = — {log (ix) + y} ber(x)+£7t bei (x) 



(2 r )! 2 


^(2r), 


(34.6) 



and that 


kei (x) — — {log(^x)+y} bei(x)-$jr ber(x) 


(34.7) 


Fig. 14 shows the variation of the functions ker(x) and 
kei(.v) over a range of values of the independent variable x. 
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The four functions ber, bei, ker and kei arc used more 
often in electrical engineering than they are in physics or 
chemistry. For a full account of their properties and those 
of their generalisation to higher order and of their applica¬ 
tion to engineering problems the reader is referred to 
N. W. McLachlan’s Bessel Functions for Engineers (Oxford 
University Press, 1934). 

§ 35. Expansions in series of Bessel functions. We 
know from § 30 that 

j* 2 ~ +:x £ + (A V - m 2 )J J m (Ax) = 0, (35.1) 

jx 2 £ +* £ + Qi 2 x 2 — n J )| J„(»x) = 0, (35.2) 

so that multiplying equation (35.1) by J n (jix)jx, (35.2) 
by J m (Xx)jx, integrating with respect to x from 0 to a 
and subtracting we find that 

(A 2 -p 2 ) j xJ m (Xx)J n (jix)dx+(n 2 - m 2 ) f J m {Xx)J n (nx) — 

Jo Jo * 

= ~XJ n (ga)J’ m {Xa)l (35.3) 

if n> — l s m> — 1* 

Putting m = n in this result we find that if A ^ 

J xJ n {Xx)J n (jix)dx 

= {jUfXa)J’ n {ga)~XJ n {tia)y K {Xa)]. (35.4) 

ir-pr 


The corresponding expression for X — p is obtained by 
putting /( = A+a, where t is small, using Taylor’s theorem 
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and then letting 8 tend to zero. We find that 

x{J n (Xx)} 2 dx 


I 


= +(l-Mfl)} 2 . (35.5) 


Suppose now that /. and ft arc positive roots of the trans¬ 
cendental equation 

hJ H (Xa)+kXaJ'„(Xa) = 0 (35.6) 

where h and k are constants. It follows then that 


i; 


xJ„(Xx)J„(nx)dx = cA, „ 


(35.7) 


where 


Cl 


= + k 2 n 2 }. (35.8) 


If we now suppose that we can expand an arbitrary 
function^} in the form 


fix) = I a.JnO-tX) 


(35.9) 


where the sum Is taken over the positive roots of the 
equation (35.6) then we can determine the coefficients a t 
as follows: Multiply both sides of equation (35.9) by 
xJJiXpc) and integrate with respect to * from 0 to a; then 

j xf(x)J„(AjX)dx = V <3, (' .Y J„{).,x)J n 0.jx)dx 
Jo ■ Jo 

from which it follows that 




xf(x)J„(Xjx)dx. 


(35.10) 
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Because of Us similarity to a Fourier series, a series of the 
type (35.9) is called a Fourier-Bessel series. 

In particular if the sum is taken over the roots of the 
equation 

j;(Aa) = 0 (35.11) 


then the coefficients of the sum (35.9) are given by 




2k) 1 

{uwy'ixy-n 1 ) 


i; 


xf(x)J n (XjX)dx. 


(35.12) 


Similarly if the sum is taken over the positive roots of the 
equation 

J„(Xa) = 0 (35.13) 




we find that the coefficients aj are given by the formula 


a 2 {-W-A} ; 


i: 


x/(x)J n (XjX)dx, 


(35.14) 


In tills section no attempt has been made to discuss 
the very difficult problem of the convergence of Fourier- 
Bessel series. For a very full discussion of this topic the 
reader is referred to Chapter XVIII of G. N. Watson’s 
A Treatise on the Theory of Bessel Functions, 2nd edition 
(Cambridge University Press, 1944). 


§ 36. The use of Bessel functions in potential theory. 
As an example of the use of Bessel functions in potential 
theory we shall consider the problem of determining a 
function z) for the half-space 0, r^O satisfying 

the differential equation 


Q QQ 

and the boundary conditions: 
(i) ^ = JIq), on z = 0; 


lot + d ^t = 0 

6q 2 o co 3z 2 


(36.1) 
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(ii) as 2 ->oo; 

(iii) -- +Ktj/ = 0 on g = o\ 

OQ 

(iv) ^ remains finite as g->0. 

We saw in § l that a function of the form ^ = R(g)Z(z) 
is a solution of equation (36.1) provided that 

-Xfz - 0 (36.2) 


and that 


d*R tdR 
dg 2 q dg 


+ Xf R = 0 


(36.3) 


where X, is a constant of separation. To satisfy the 
boundary condition (ii) we must take solutions of equation 
(36.2) of the form 

Z = 


and to satisfy the condition (iv) we must take as the 
solutions of equation (36.3) functions of the form 

R — J o(Aj£>) 

since the second solutions Y 0 (X,q) would become infinite 
in the region of the axis q = 0. 

The differential equation (36.1) and the boundary con¬ 
ditions (ii) and (iv) are satisfied by any sums of the form 

<ACp. z)= £ a t e~^J 0 (X t o) (36.4) 

i 

where the a { and A f are constants. But if we are to satisfy 
the boundary condition (iii) we must take the sum over 
the positive roots of the equation f 

o(Ajfl) -f- kJ o(Affl) — 0. (36,5) 

+ For properties of the roots of this equation see example 16, p. 144. 
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$36 


The solution is determined therefore if we can find constants 
a t such that condition (i) is satisfied* i.e. such that 

/(c)=lMo(^)- (36.6) 

i 

From equations (35.10) and (35.8) we see that we must take 


a, = 


2 Xf 


a 2 (Xf + K 2 ){J 0 (Xfi)Y 

Hence the required solution is 


j>' 


e'AeVo V-,e'W. (36.7) 




(36.8) 


where the sum is taken over the positive roots of the 
equation (36.5). 

If, instead of the boundary condition (iii), we had the 
condition \jf = 0 on q = a then it is easily seen that the 
solution would have been 


1KP. *)= - Z 1 r e'f(eVoW)d6' (36.9) 

or I {JiMY Jo 

where the sum is taken over the positive roots A s of the 
transcendental equation 

J 0 (la) = Q . (36,10) 


For example suppose that ^ satisfies the conditions 


on 


= 0 on q — a, as 2 ->oo* \j/ - #q(« — G ) 

z — 0* 0then the solution to the problem is given 
by equation (36,9) with ffa') - if/ Q (a 2 —Q t2 )> By equation 
(32,5) we have 


i; 


' q'Kq'VMW = A -VbLj^Xfl)- JoQ.fi) 


a f 

4#o 

A? 


A? 


Jittfi) 
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since X l is a root of equation (36,10). Thus the required 
solution is 


v(p, 

a t 


WtM ' 


(36.11) 


Tables of the first forty zeros Cj of the function J 0 (£) 
with the corresponding values of are available f so 
that it is convenient to express results of the kind (36.11) 
in terms of them. It is readily seen that in this case 


Wp> 2 ) = 8a Vo 


• z?mi) 


where ( — zja and a = Qja . 


(36.12) 


§ 37. Asymptotic expansions of Bessel functions. In 
certain physical problems it is desirable to know the value 
of a Bessel function for large values of its argument. In 
this section we shall derive the asymptotic expansion of 
the Bessel function of the first kind J R ( X ) and merely 
indicate the results for the other Bessel occurring in mathe¬ 
matical physics. 

We take equation (28,1) as our definition of the function 
J n (x). Applying the theory of functions of a complex 
variable it is readily shown that this definition is equivalent to 

'■» - rfej {jV<W+ JV<W} 

(37.1) 

where L x is the straight line ${t) - — 1 in the upper 
half of complex /-plane and L 2 is the corresponding part 
of the straight line $(t) — + L By changing the variable 

t A. Gray, G. B. Mai hews and T. M. MacRobcrt, A Treatise 
art Bessel Functions and Their Applications to Physics, 2nd edition 
(Macmillan, 1931), 


$37 
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from ( to u = ix(l — t) in the first integral and lo 
u — -ix(l-t) in the second we see that 


where 
/„(*) - 



(37.2) 



and j*(x) denotes its complex conjugate. Expanding 
(I + j w/2jr>' 1 ~* by the binomial theorem and integrating 
term by term we find that 

/,(*) = t F 0 (±+n t i- n; (37.3) 

If we adopt Hankel's convention of writing 


(n, r) = (-ir 


U-hUI + h), 

r! 


in equation (37.3) and substitute the result in equation 

(37.2) we find that for large values of * the asymptotic 
expansion of the Bessel function J n (x) is 


jcos(jc—inti—in) £ ^ 

\jnxl r = o (2x)* r 

. , , , , £ (— 1 /( 11 ,2r+1)| 

-.in(x-fr,-M > g a (2x) „„ }• 


(37.4) 


The corresponding expansion for the Bessel function of 
ihe second kind is found to be 


(37 ' 5) 
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Substituting these asymptotic expressions in equations 
(30.18) and (30.19) we find that as x -» oo, 

07.6) 

where j„(x) is given by equation (37.3). 

In certain problems only a very crude approximation 
to the behaviour of the Bessel function is desired. In these 
circumstances the following formula are usually sufficient: 





irni-£jc). 


Y n (x) 



sin (jt — inn — Jir); 


(37.7) 


WKx) 


~ITe» 

V TEX 


-inn?- 


//< 2 >(jc)~ (37.8) 

V TEX 


Similar formulae exist for the modified Bessel functions. 
Proceeding in the same way as in the establishment of 
equation (28,1) we can show that 


/„(*) = 


V*n»+i) 


(i*r 


|* e ±x, (l — 


t 2 y-*dt 


which becomes 


«*) = 


1 

V(2TEx)r(n+i) 




by a simple change of variable. By a method similar to 
that employed above to obtain the asymptotic expansion 
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of J m {x) wc can then show lha* if — |ir<arg x <\jty 

I(x) _ f i-mn,r) e-* + «+** f (a, r) 

V(2 nx) M o (2xX V(2 ttx) ,£o (2x) r 

and that if - f ?r < arg x < Jtt the factor exp {—x+(n -f £)jri} is 
replaced by exp { — x— (n+\)ni}. The corresponding formula 
for the modified Bessel of the second kind is 



as x-mx>. 


Examples IV 

(1) Making use of example 2 of Chapter JI and of 
the expansion (25,8)* expand cos (x sin 0 ) as a power series 
in sin 0 in two ways. Hence by equating powers of 
sin 2 ’ 0 show that if s is a positive integer 


x 2 * = 2 2i+l 



(n+s— 1 ) 
(n-s}\ 


-Jlnix)* 


Derive the corresponding result for x 2J+l and show that 
the two results may be combined into the single formula 


x r = 2' £ (r t? 1)1 (r+ 2 n)J r+ 2 n( x ) 
n *=0 n 1 


(r = I t 2* 3, 

(2) If u and £ denote the eccentric and mean anomalies 
of a planet show that 


« 1 

cos (m<) = n £ - J m _„(me) cos (mQ, 

m = - go ttj 


to j 

sin (nu) = m £ - J„-„(me) sin (mf). 

m ■ -a fB 
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(3) Making use of the expression for P n (cos 0) given 
in example 19 of Chapter HI show that 


E P„(cos 0) = e r ca ' °J 0 (r sin 0). 

- m i " 


n = 0 ill 

(4) Show that 

(i) 8j;(z) <= J„ _ 3 (z) - 3/„ _ ,(z)+ 3J„ + t (z) - J B+3 (z); 

(ii) 4JJ(z) +VJ(z)+/ 3 (r) = 0. 

(5) Prove that 

'cos(AJ + i)u du 


2 r~i n J 0 


COS yj(x 2 — u 2 ) 


and deduce that 


~+ £ (-l)7 0 (rx) = 0. 

i r ** 1 


(6) Show that the curve with freedom equations 

x = t — sin t $ y — cos t/(l — cos /) 

may be represented in the interval 0<t<n by the Fourier 
scries 

y = 2 E 2„(n) cos (nx). 

fl « 1 

(7) Prove that 

fi, * r = (^)* Jo (2 " + ^^wwceos 0 )- 

(8) Prove that 

lim P„ ("cos -) = J 0 (x). 

\ W/ 
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(9) Show that 


(0 E = J 0 { y J(a 2 -2ax)}; 

a ° 0 ft! 

(ii) E ^ = (.v + fl)-i'"7„,{2 v /(.v + n)}. 

rt = o n! 

(10) Prove that if n is a positive integer 

i„(Vx) = (~ 1)" X* £; J 0 (2jx). 

(11) If ,v >o show that 
e“ I “*COS (x sin B)d0 = J 0 {-J(x 2 — fl 2 )}. 


£ 


(12) Prove that, if — 1 <x< 1, 

1 


70 ~* 2 ) 

Deduce that 


— ^ i 0 (»!n) cos(tmix). 


m - i 


M*) = 


sin x 


x l 


(i+j** £ fchwgol 

hi = i x—mn ) 


(13) Prove that 

7;(x)y„ v (x)-7 v (x)j'_ v (x) = 

X 

where A is a constant, and, by considering the series for 
J v {x) and J v (x) when x is small show that A — (2 In) sin (vn). 

(14) Show that the complete solution of Bessel’s 
equation may be written in the form 




■ w f 




${US)} 2 

where A and B are arbitrary constants. 
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(15) Show that the complete solution of lhe differential 
equation 

d 2 y . 1 


is 






where c? — 4x*j21 and A and B are arbitrary constants* 

(16) If a and b are real constants show that the roots 
of the equation 

axJ‘ n (x) + bJ„(x) = 0 

are simple roots except possibly the root x — 0. 

Show also that the equations J„(x ) = 0, J„(x) — 0 have 
no roots in common except possibly at = 0. 

(17) If a > 1 andm + « + l>0, prove that 


i: 


•""WO*" - ** 


-fr- 




where Q„(x) denotes the associated Legendre function of 
the second kind* 


(18) Prove that 

U*+y)= t Ux)in-*,(y) 

m » 0 

+ 1 {Ux}I n+ni (y)+I H+ Jx)IJy)}. 

IH 14 1 

(19) Show that 

J‘ e*>‘PJ(n)dv= 

(20) Prove that 

f Jn{kx)i” = e" f krru x) 

n « ~ oo it™—eo 
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and deduce that 
(0 Ure iB )= t 


m ™ 0 


mi 


(ii)«*)- I 

(21) Using the expansion of the last question prove 
that J r (ae ix +be ifi ) is the coefficient of t r in the expansion of 

By putting R - a cos a+b cos /?, 0 = a sin a + b sin jff, 
f} = ot + fl — n prove Neumann's addition theorem 

UR) = 1 J n+m (a)J m (b)e^ t 

\ Ct~^DB / m ~oo 


where R 2 = a 2 +b 2 — Tab cos 0. 
(22) Prove that 


-If' 

*Jo 


J a (Rx)dO, 


J 0 (ax)J 0 (bx) 
wliere R 2 = a 2 +b 2 —Tab cos 0, and deduce that 

k m 


I, 


J 0 (ax)J 0 (bx)e-"dx - 
o n^ab) 


where 


K{k) = , k 2 = 

Jo sin 2 


4a b 


I, 


Prove, in a similar way, that 

” J 1 (ax)J l (bx)e^dx = — i— 
o iiky/iab) 


(a + b) 2 +c 2 
{(l~ik 2 )K(k)-E(k)} 
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where k and K{k) are as defined above and 


= j N /(t — A: 2 si 


m 

(23) Show that 

J m (x)J n (x) = - 
n 


sin 2 <j>)d<f>. 


J\- m (2xsi 


sin 0)e ltt * mV>i dO 


and hence that 


{•/„(*)} 


f 

n Jo 


J 2 „( 2x sin 0)d0. 


Deduce that 


R(*>} 2 = (±*) In I 

(24) Prove that 


(2n+2s)! 
o .5!(2n + s)!(n + s)t 


<-± * 2 Y- 


r ~ «p (- ° 4 /) '• ( 0 ) 

(25) We define the Bessei-integra! function of order n 
by the equation 

= - f' n = 1,2,3, .... 

Jo » 


Prove that if n is even 


m 1 f" 
* Jo 


cos (ji(J)cj(x sin &)d0, 


where c/(x) denotes the cosine integral, and derive the 
corresponding expression when n is odd. 

Show that: 

(i) !*{«*■**-«} = £ f"J(.(x); 

It = “ <Xt 

(ii) = d,(x)/x; 
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(iii) (II - 1 v, (x) - (ii +1 +, (x) = 2nJ i'„(x) ; 

(iv) ci(x) =■ J/ 0 (x) -2 Ji 2 (x)+ 2J/ 4 (x)— 

(v) si(x) = 2/i,(x)—2di 3 (x)+2Ji s (x) —...; 

(vi) di 0 (x) = y +iog(+x)- ^ 2 f,(l, 1; 2, 2, 2; -x 2 /4). 

O 


(26) If 


= — f‘x 2 “ +1 

JlWJo 


J 0 (Xx)4x 


where X is a zero of d 0 (x), show that 

\nf_ 

;. 2 


I — I 

m X 


, 4m 2 , 4 2 m 2 (m—l) 2 

1 —~ Hr --- 




4“m 2 (»i —1) ...2 .1 
X 2m 


Show that 


1 - X 2 " = t, a i J o(2|*), 0 £ x <j i, 

i - 1 

where X lt 2 2 , ... are the positive zeros of J 0 (x) and 

2 

£1, — 


1 "w (1 u ' 


Putting nt = 1, 2, 3, .... and rearranging the series show 
that 

r = i XfJjiX,) 8 l } ' 






















"M3 
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Deduce that 


_i__ i, * i = jL- y i = _19_ 

i tfJtHi) ~ H S j = i tfJAld ^128’ iVt XlJ^Xt) ~ 460S 

(27) If a it a 2 * *** denote the positive zeros of J r ( x ) 
show that 


£ —r-~—r | x/(x)JJa,x)dx = j [ .v v+ %x)dx. 
*- i «/J.+ i(«i) Jo Jo 

Taking J{x) = * v , x' l+1 show that 

S = _J_. y i = 1 

( =t« ( 2 4(v+l)’ i= i af 16(v+l) 2 (v+2)' 

(28) Putting v = 0, J[x) — in (27), show that 



f = t A; 4 (/h 4-1) 


where 2,, A 2 , ... are the positive zeros of J 0 (x) and /„, is 
defined in (26). 

Substituting the value for I m derived in (26) and using 
the fact (obtained by putting v = 0 in (27)) that £Aj" 1 = 
show that the sum 

Jp 

satisfies the recurrence formula 


m /i\in-r+l e 

s 2m+2 = t , c- ir _r j r , 2 ; 1 , nI 

r -1 V 4 / [(m-r+1)!] 2 

+( - ir (i)‘*W 


Deduce that 


f i = i. f i r i. f JL_ 

t-i Af 32’ j-i Af 192’ ,£■ Af 


!(tn +1)! 
11 


12288 
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r 


(29) Show that 


£ | 0 g x . 

A|Vf(A,) 


Multiplying both sides of this equation by x and integrating 
term by term* show that 

&M<X)- M, - 2,oex) ' 


and hence that 


1 


t = i A fj j(Aj) 


= i 


























CHAPTER V 

THE FUNCTIONS OF HERMITE AND LAGUERRE 


§ 38, The Hermite polynomials. The Hermite poly¬ 
nomial H n (x) is defined for integral values of n and all 
real values of by the identity 


ir -0 n\ 


( 38 . 1 ) 


If we write 

/(*, 0 - e 2 *-** = e*V ( *“ r >* 
then it follows from Taylor's theorem that 

HM-(a) .az.-*-*] 

WJ.. o Lai" J,_. 

Now it is obvious from the form of the function 
exp {— (jc—/) 2 } that 


and so we have the form 


KM = (- 1)V* ~ (e~*) 
ax 


( 38 . 2 ) 


for the calculation of the polynomial 

It follows from this formula that the first eight Hermite 
polynomials are: 

H 0 {x) - 1, 

H,(x) = 2x, 


]» 
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H 2 (x) = 4x 2 - 2. 

H 3 (x) = 8x 3 — 12x, 

H 4 (x) = 16x 4 —48x 2 + 12, 

H s (x) = 32x 5 —160x 3 + 120x, 

H b (x) = 64x 6 —480.v 4 + 720x 2 - i 20, 

H n {x) = 128x 7 -1344x 5 + 3 360x 3 - 1680.x. 

In general we have 

H„(x) = (2xy- 

+ «(« — 1)(H -■ 3)(H - 4) Q^n-4 , 

or, in the notation of Section 12, 

KM = (2x)" 2 F 0 ^-in, i-in;- (38.3) 

Recurrence formulae for the Hermite polynomials follow 
directly from the defining relation (38.1). If we differ¬ 
entiate both sides of that equation with respect to x we 
obtain the relation 

#t = o n l 

from which it follows directly that 

2nH n _ t x = H'(x)’ (38,4) 

On the other hand if we differentiate both sides of the 
identity (38.1) with respect to / we obtain the relation 


2{x-t)e 


2ijc“( 2 _ 


= I 


r -i 


i (n-l)! 


which can be written in the form 

a KM a ff n (x) , _ a 
2x £o — ‘ \h~nT 1 -„h(n-iy. 


KM (n _, 
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to yield the identity 

2 xH„(x) = 2 nH n . 1 (x)+H n+ ,(x) (38.5) 

by the identification of coefficients of t". 

Eliminating j(jr) from equations (38.4) and (38.5) 
we obtain the relation 

fftx) = 2xH„(x)-H n+ t (x). (38.6) 

Differentiating both sides of this identity we find that 

H;(x) = 2xH’ n (x) + 2 H n (x) - HU ,(x) 

and, by equation (38.4), H„ + l (x) = 2(n+ so that 

H n „(x) - 2xH' n {x)+2nH„(x) = 0. (38.7) 

In other words y = H£x) is a solution of the linear differ¬ 
ential equation 

y*—2xy'+2ny = 0. (38.8) 


§ 39. Hermite’s differential equation. We saw in the 
last section that //„(*) is a solution of the differential 
equation (38.8). Replacing the integer n in that equation 
by the parameter v we obtain Hermite’s differential equation 

pL-l, ^+2,,-0. (39.1) 

ax ax 


if we assume a solution of this equation in the form 


y= I 

r 0 

and substitute in the equation (39.1) we obtain the recur¬ 
rence relation 


a r+ 2 — 


2(r+Q-v) 


a r 


(39.2) 


(r + g+2)(r+g4*l) 
on equating to zero the coefficient of Equating to 
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zero the coefficient of we obtain the indicia! equation 


= 0. (39.3) 

Corresponding to the root q — 0 we have the recurrence 


relation 


a r+1 — 


2(r — v) 
(r+lXr+2) 


a r 


(39 .4) 


which gives the solution 


ir fr v , (, 2v , 2 2 v(v—2) 4 2 3 v(v—2)(v—4) 6 \ 

l--* + —jp- x 6 T '- x + -) 


(39.5) 


where a x is a constant. 

Similarly, corresponding to the root g = 1 of the indicial 
equation we have the recurrence relation 

2(r + l~v) 


a r + 2 — 




(r+3)(r+2) 
from which is derived the solution 

2 2 (v-l)(v-3) 


(39.6) 




(, 2 (v — 




5! 




where a 2 is a constant. The general solution of Hermite’s 
differential equation is therefore 

y = yi(x}+Y2{x)* (39.8) 

For general values of the parameter v the two series 
for >*i(a-) and y 2 (x) are infinite. From equations (39.4) 
and (39.6) it follows that for both series 


a r + 2 ~ a r> as 

r 


(39.9) 


If we write 


exp (x z ) = 

L 
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then 2 

b r+2 ~ - b„ as r-»oo. (39.10) 

r 

Suppose that a N jb N is equal to a constant y, which may be 
small or large, then it follows from equations (39.9) and 
(39.10) that, for large enough values of N, % +2m /£> A r +2m --'y. 
In other words the higher terms of the series for y,(x), 
y 2 (jr) differ from those of exp (x 2 ) only by multiplicative 
constants y lt y 2 , so that for large values of | x |, 

y 1 (x)~y 1 e*\ y 2 (x)~y 2 e** 

since for such values the lower terms are unimportant. 

We shall see later (§ 41 below), that in quantum 
mechanics we require solutions of Hermite’s differential 
equation which do not become infinite more rapidly than 
expQx 2 ) as |x|->oo. It follows from the above con¬ 
siderations that such solutions are possible only if either 
y t (x) or y 2 (x) reduce to simply polynomials and it is 
obvious from equations (39.5) and (39.7) that this occurs 
only if v is a positive integer. For example, if v is an 
even integer n we get the solution 

y(x)=c//„(x), (39.11) 

where c is a constant, by taking a 2 = 0, 




= (-!)*" 


n! 

ft") I 


c. 


Similarly, if v is an odd integer n we get the solution 
(39.11) by taking a, = 0 and 


a 2 = (-!)*"-* 


2n! 

ft"-i)i 


c. 


Hermite’s differential equation therefore possesses solutions 
which do not become infinite more rapidly than exp ftx 2 ) 
as | x |-*oo if and only if v is a positive integer n. When 
this is so the required solution of Hermite’s equation is 
given by equation (39.11). 
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§ 40. Hermite functions. A differential equation closely 
related to Hermite’s equation is 

< pL+(X-x 2 )^ = 0. (40.1) 

dx i 

If we transform the dependent variable from ^ to y where 

</, = e ~ ix2 y ( 40 . 2 ) 

and put A = 1 -H 2v then it is readily shown that y satisfies 
Hermite’s equation (39.1). The general solution of equation 
(40.1) is therefore given by equations (40.2) and (39.8) 
with y 1 (x), y 2 (x) given by equations (39.5) and (39.7) 
respectively. 

The argument at the end of the last section shows that 
the equation (40.1) possesses solutions which tend to zero 
as j a: | —► oo if and only if the parameter X is of the form 
1 + 2/t where n is a positive integer. When X is of this 
form the required solution of (40.1) is a constant multiple 
of the function T„(x) defined by (he equation 

TO = e-**'H n (x) (40.3) 

where H„(x) is the Hermite polynomial of degree n. The 
function ¥„(x) is called a Hermite function of order n. 

The recurrence relations for x F„(x) follow immediately 
from those for f/„(x). For instance equation (38.4) is 
equivalent to the relation 

2»n-iW = x'F n (x) + H'^) (40.4) 

and equation (38.5) is unaltered in form so that 

2xTO = 2rt x F n _ ,(x)+T n+1 (x). (40.5) 

Eliminating 2n'¥„- t (x) from equations (40,4) and (40.5) we 
have the relation 

TO = xTO-’W*)* (40.6) 

From the point of view of mathematical physics the 
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most important properties of Hermite functions concern 
integrals involving products of two of them. In establish¬ 
ing most of these properties the starting point is the 
observation that the function *f? m (x) satisfies the relation 

i r a +(2n + l-jc 2 )'F 11 = 0, (40.7) 

as is obvious merely by substituting 2n+ 1 for A in equation 
(40.1). 

Writing down the corresponding relation for T m , 

V*+ (2m +1 ~x 2 )^ m = 0, (40.8) 

multiplying it by *P B and subtracting it from equation 
(40.7) multiplied by f,, we obtain, as a result of integrating 
over (—co, oo), the relation 

2{m-n) f“ '¥ m '¥„dx = f" C¥ m W;~^Jdx. 


Now an integration by parts shows that the right-hand 
side of this equation has the value 


- r (Vn-VWJdx 

L J-co J 

and, if we remember that, for all positive integers n, 
^ B (x)-*-0 as j x j->co, we see that this has the value zero. 
Hence if we let 


».„ = r vjxyrdOdx 

J — to 

L. p = 0, if m ^ It. 


we see that 
In particular 

^n - 1, p+i = 0 

so that from equation (40.5) we have 


J: 


(40.9) 


2xV H (x)'¥ ll . 1 (x)dx = 2 n I,_ 


(40.10) 
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Now if in equation (40.3) we substitute for H h (x) from 
equation (38.2) we have 

^(x) = (-ire* xl ~(e~ xl ) (40.11) 

dx 

so that the left-hand side of equation (40.10) is equal to 


-r 




and an integration by parts shows that this is equal to 

fl4-f n +i i 1 

i.e. to I Hw n . Hence from equation (40.10) we have 

^,. = 2 m n - ltn - t . 

Repeating this operation #i times and noting that 




_ r 

>. 0 - 

J —0 


e **dx = sjn 


we find that 

I n ,» = 2*n\jh. (40.12) 

Combining equations (40.9) and (40.12) we have finally 

= (40.13) 

The evaluation of more complicated integrals can be 
effected by combining this result with the recurrence 
formulae we have already established for the Hermite 
functions. For instance, it follows from equation (40,5) 
that 


r 


x'? m (x)'f'„(x)dx = nl mi , 


showing that 


i: 


x v f m (x) v f , n (x)dx = 0 if m ^ n ± I (40.14) 
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and that 

*TJixyv m +i(x)dx = 2"(it +1) !>■ (40* 15) 

Similarly, making use of equation (40,4) and equations 
(40 A 3*15) we can show that 

0if m ^n±\ 

2"" , jil7^if m = n-l 
— 2 "(n +1) if m = n +1. 

§ 4L The occurrence of Hermite functions in wave 
mechanics* The Hermite functions which we have dis¬ 
cussed in the last section occur in the wave mechanical 
treatment of the harmonic oscillator *t Although this is 
a very simple mechanical system the analysis of its properties 
is of great importance because of its application to the 
quantum theory of radiation* 

The Schrddinger equation corresponding to a harmonic 
oscillator of point mass m with vibrational frequency v is 

pL + (W-htWx 1 )# = o, (41.1) 
dx 2 hr 

where W is the total energy of the oscillator and h is 
Planck’s constant. The problem is to determine the wave 
functions $ which have the property that 

(i) ^-*0 as | x I^oo; 

(ii) f" [ * |Hx = l. 

J — 00 

If we let 

* = r 

2n \j mv 

t N. F. Mott and 1. N. Sneddon, Wave Mechanics and Its 
Applications (Oxford, 1948), p. 50. 


i 


'¥ m (x)'¥'(x)dx = 
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then the equation (41*1) becomes 



and the conditions (t) and (ii) become 
(0 iH0 as | £ |-»oo; 

The argument given at the beginning of § 40 shows 
that equation (41*2) possesses solutions ^ which satisfy 
the condition (ii') if and only if the parameter (2^/Av) 
which occurs in the equation takes one of the values 
1+2 n where n is a positive integer* In other words 
solutions of this type, which are known by the probability 
interpretation of the wave function ^ to correspond to 
stationary states of the system can exist if and only if 

W=hv(n + i) (41*3) 

where n is a positive integer* When this is the case the form 
of the wave function \j/ is known from § 40 to be 


* = C%(0 (41*4) 

where C is a constant. Applying condition (ii 1 ) and 
equation (40*12) we see that 




4umv^* 1 
~h~j 2*«(n\) r 


Thus the wave function corresponding to an admissible 
energy (h+|)/iv is 




WVjwo 

. h ) 2 1 *(n !)*’ 




(41.6) 


In quantum theory the matrix elements (n \ x ] p) 
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defined by the equation 

(« | x | p) = I x<l/ n (x)4/ p (x)dx 
J - <C 

are of considerable importance in the case of the harmonic 
oscillator. In terms of the variable £ we have 

(« M p) = 7T- P ZMQUOK* 

4n 2 mv J_„ 

so that substituting from equation (41,6) we have 


('< * \P) = - 


h f° 


ji V 4jtwv 


J-, 




It follows then from equations (40J4) and (40.15) that 

(«|x|p)«0 ifp#n±l (41.7) 


and that 




(41.8) 

(41.9) 


§ 42. The Laguerre polynomials. The Laguerre poly¬ 
nomials L n {x) are defined for n a positive integer and x 
a positive real number by the equation 


exp 



cc 

0-0 I 

s- 0 


L„(x) fB 
n! 


(42.1) 


Expanding the exponential function we see that 


(l-.)-exp(- j^) 


f. (-1 YxY 

r = 0 r 1(1 ~ty +1 
£ g (-l) r (r+l)s ;:Y+I 
r - o i - o r!s! 
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The coefficient of /" in this expansion is 


y (-mr+a-r ^ 

r = o r !(w — r)! 


Using the relations 


( r +l)n-r= 


r!’(n-r)! »! 

we see that this sum can be written in the form 


("~ n )r v r 

r-o (r!) 2 

Identifying the coefficient of t n with L„(x)jn ! and adopting 
the notation of § 11 we see that 

L^x) = it! I; x). (42.2) 

It should be observed that L n {x) is a polynomial of degree 
n in ,v, and that the coefficient of x? is (— 1) B . 

A useful formula for the polynomial L n (x) can be 
obtained by finding a new representation for the confluent 
hypergeometric function on the right-hand side of equation 
(42.2), By Leibnitz’s theorem for the nth derivative of a 
product of two functions we have 

s'DVe'^e 1 V (-l) r ^^(D n “ r x n )(B r O 

r -o rl 

where D denotes the operator d\dx. Using the relations 
e*D r (0 = (- If, D n ~'x n = n Ix'/rl 

we see that 

e*D"(x"e~*) = nl £ (42.3) 

r-o (r!) 

It follows immediately from equation (42.2) that 

L n (x) = e*~(x n e-*). (42.4) 

dx" 
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The first five Laguerre polynomials can be calculated 
easily from this equation; we find that 

M*) = l, 

Li(x) = 1 -x, 

L 2 (x) = 2—4x+x z , 

L 3 (x) = 6 —18jr+9jc 2 —jc 3 , 

Z- 4 (a) = 24-96a+72at 2 — 16a' 3 + a*. 

Equations (42.4) can be used to show that the functions 

Ux) = ~ e~**L a (x) (42.5) 

it! 


form an orthonormaL system. From (42,4) we have as a 
result of m integrations by parts 

f® f« A* 

e’VIJj t)dx= xT —{x n e-*)dx 

Jo Jo tlx" 

and this is 2 ero if n>m. Since L m (x) is a polynomial of 
degree m in x it follows that 



e x L m (x)L„(x)dx = 0 if m/«. 


(42.6) 


Since the term of degree n in L„(x) is (-l)V it follows 
that, when m = n. 


e~*{LJix)} 2 dx = (-1)"1“ fVLi# 


t'.x a e-*dx 


=>(«!) 2 
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Combining this result with equation (42,6) we find that 


i: 


<t>n,(x)<f>n(x)dx = 


(42.7) 


showing that the $*s form an orthonormal set. 

Recurrence formulas for the Laguerre polynomials may 
be derived directly from the definition (42.1), Differ¬ 
entiating both sides of this equation with respect to t we 
obtain the identity 


x ( xt\ 

- (T^o 3 p \ ~t) 


=(l-o i - £ LXxr 

u = 0 (n — 1)1 n~Q 


nl 


which may he written in the form 

, £ ^£ L -^g 

n * 0 u ! n = 0 (/I — 1 )! 


-(l-o t ^^ = 0. 


n! 


Equating to zero the coefficient of t n in the expansion on 
the left we obtain the recurrence relation 

W*)+(*“^-l)U*) + » a ^-i(x) = 0- (42-8) 

Similarly if we differentiate both sides of (42.1) with respect 
to x we obtain the identity 

, l AM im r .o 

it *> 0 ft! ii *0 fl! 

which yields the recurrence relation 

Itfx)-ii£i. 1 (x) + »I* M Cx) = 0. (42.9) 

Differentiating equation (42,8) twice with respect to x and 
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replacing n by n +1 we find that 

W*)+(*-2«-3)£;+i(j0 

+(n+l) 2 i;(x)+2i; +1 (x) = 0. (42.10) 

Now from (42,9) 
and hence 

U+iW-(» + «{!«*)-W}. 

A similar expression for Z*+ 2 (x) in terms of L£x) and its 
derivatives can be readily obtained. Substituting these 
values of L" + l (x) and jE4+*(x) in equation (42.10) we find 
that 

xia*)+(l - x)mx)+nLJtx) = o, (42. 11) 

§ 43, Laguerre’s differential equation. Equation (42.8) 
shows that y — AL n {x) is a solution of Laguerre’s differ¬ 
ential equation 

X +(1—x) ^ +vy = 0 (43.1) 

dx ax 

in the case in which v is a positive integer n . If we put 
y — a — —v in equation (11.2) we see that it takes 
the form (43.1) so that it follows from equation (IJ.4) 
that one solution of equation (43 J) is 

yi{x)= l F l (-v;l;x) t (43.2) 

Similarly we see from equation (11.8) that the second 
solution is 

. log JC+ £ CrX' (43.3) 

r ■ l 

where the coefficients c r are defined by equations of the 
type (11.9). The general solution of Laguerre’s differential 
equation may therefore be written in the form 

y = Ay t (x)+By z (x) 
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where A and B are constants and y^x), y 2 (x) are defined 
by equations (43,2) and (43.3) respectively. 

If we are interested only in solutions which remain 
finite at x — 0 it is obvious from equation (43.3) that we 
must take the constant B to be zero. Further if a r is the 
coefficient of x r in the series expansion for yi(x) it is easily 
shown that a r +ila t ~r~ L . As the result of a discussion 
similar to that advanced in § 39, it follows that if v is 
not an integer 

y 1 (x)'^e x as x-*co. 

if, therefore, we are looking for solutions which increase 
less rapidly than this we must take v to a positive integer, 
in which case y t (x) reduces to a polynomial. 

The equation (43.1) possesses a solution which increases 
less rapidly than e x as x-+ oo if and only if the parameter v 
occurring in it is a positive integer, n say. If it is also required 
that the solution shall remain finite at x = 0, the solution 
is of the form 

y - (43.4) 

where A is a constant and L n {x) is the Laguerre polynomial 
of degree rc. 


§ 44. The associated Laguerre polynomials and func¬ 
tions, If we differentiate Laguerre’s differential equation 
m times with respect to x we find that it becomes 


im + 2 „ jm+ 1 Jm 

x + 1 -•*) + (”■ m ) 


dx 


dx” 


which shows that Z£(x) defined by 

K(x) = ~L n (x), (nSm) 
dx m 

is a solution of the differential equation 
„ d 2 y , , v dy 




(44.1) 
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The polynomial L™(x) defined by equation (44.1) is called 

the associated Laguerre polynomial. It follows from 

equation (42.2) that it may be presented as a series by 
the equation 

L%X) = ^ !F,(-n +m; m +1; x) (n £ m). (44,3) 

m l(/i — m)} 

Similarly equation (42.4) leads to the formula 
d m ( d n 1 

L m n { x ) = _ \e x — (x n e~ x ) L (44.4) 

dx m \ dx n J 

The simplest associated Laguerre polynomials are: 

*.!(*)=-!. 

Lj(x) = -4 + 2.X, L\(x) = 2, 

L\(x) = —18 + 1 Sjs — 3jc z , Lfrc) = 18-6*, Ll(x) = -6, 

Lj(x) = —96+ 144x —48* 2 + 4x 3 , LjJ(x) = 144-96.v + 12x 2 

Ll(x) = - 96 + 24*, Li(x) = 24. 

The definition (44.!) for the associated Laguerre poly¬ 
nomial is the one usually taken in applied mathematics. 
In pure mathematics the function 

L m (x) _ , Fl (_„ ; m +1 ; *) (44.5) 

mini 

which is a solution of the differential equation 

* + (m +1 —*) — +«y = 0 

dx dy 

is often taken as the definition of the associated Laguerre 
polynomial t so that care must be taken in reading the 

t See. for instance. B. T. Copson, Functions of a Complex Variable 
(Oxford University Press 1935). p. 269. 
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Literature to ensure that the particular convention being 
followed is understood. 

It is readily shown from equation (42.1) which defines 
the generating function for Laguerre polynomials that the 
associated Laguerre polynomials may be defined by the 
equation 

(-l) m rexpf - =(l-0'" +1 y ^l t \ (44.6) 

V 1 — 1/ FI ^ ™ ni 


This identity can then be used to derive recurrence relations 
for the associated Laguerre polynomials similar to those of 
equations (42.8) and (42.9) (cf. examples 6 (it), (iii) below). 
The Laguerre functions R B t(x) are defined by the equation 

RJx) = e'VC/W (»£/+!)• (44,7) 


If in equation (44.2) we replace »i by 2/+1, a by « + l 
and y by e ix x~ l R we see that the function R n! (x) is a 
solution of the ordinary linear differential equation 


(FR + 2dR _ fl 
dx 2 x dx (4 



«±»j. 


0 . 


(44.8) 


In most physical problems in winch this equation arises 
it is known that / is an integer. By reasoning similar to 
that outlined in § 43 it can readily be shown that the 
equation (44.8) possesses a solution which is finite at 
x = 0 and tends to zero as x-.oo if, and only if, the para¬ 
meter n which occurs in it assumes integral values. When 
this does occur the solution is AR nt {x) where the function 
^ni(*) is given by equation (44.7) and A is an arbitrary 
constant. 

We shall now evaluate the integral 



which arises in wave mechanics. From equation (44.6) we 
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have the identity 

y y ^+iH*Hv + |(*) jiH-y-ff 

fi ^ J41 n' ^i+1 (« 4- /) !(n' + /)! 

f xt xx 1 
exp ^-—> 

l 1 — / 1 —tJ ^21+ I—214-1 

-- ^ - TTrrrr * t 


(1-0 3,+2 (1-t) 


21 + 2 


from which it follows that 


<o oa Pco 

.L,,JL.o&fco iJ. •^"‘■'es’wwew^ 

j:— P {— 


_0t) 

\2l + 2 


(l-0 (1 —t) 

This last integration is elementary and gives for the expres¬ 
sion on the right 

( 2/+2)i(tr)* l + 1 (l —Q(l—t) 

(1-(t) z,+3 

and by means of the binomial theorem we may expand 
this function in the form 

(1 — t —t — ft) £ ( 2 /+r+2)l ( ft )H+ r+l 
t = o r! 

Now the coefficient of (tr)’ ,+ t in this expansion is 

2»{Q2 + Q1} 

(n—1—1)! 

and this is equal to /„,/{(«+/)!} 2 . Hence 

Jo (n — i — 1)! 


(44.9) 
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§ 45, The wave functions for the hydrogen atom. We 
shall conclude this chapter by discussing the motion of a 
single electron of mass m and charge — e in the Coulomb 
field of force with potential 


V(r) = ~ 


Ze 2 


r 


due to a nucleus of charge Ze. In a first approximation we 
may treat the mass of the nucleus as', infinite and in this 
case the wave function i// of the system is governed by the 
Schrodinger equation 

VV+^^+^<M 0 (45.1) 

and the conditions: 

(i) ^{r, 0, <f>+2n) = 0, <f») for all r, 9, <j>; 

(ii) i jt must be bounded in the range O£0 £ji for all 

r, $; 

(iii) ijf—tO as r-> oo; 

(Iv) if/ remains finite as r-*0; 

(v) if/ is normalised to unity, i.e, J| ^ | 2 dt = l where 

the integral is taken throughout the whole of space. 
W is the total energy of the system. 

The equation (45J) may be solved by setting 

ip = (45,2) 

where, by the method of separation of variables we have 


4-u 2 «D = 0, 

d<j> 2 


(453) 


1 d 
sin 6 d9 


(s i " e f) + {' ( ' +1) - ^} 0 -°- (45 4) 


M 
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R = 0. 


(45.5) 


To satisfy condition (i) we must choose as a solution of 
(45.3) a function <t>(0) such that <t>(<t>+2n) = for 
all <j>. Thus u occurring in equation (45.3) must be an 
integer and a convenient solution will be 


0 = Ae ,u * 


(45.6) 


where A is an arbitrary constant. Equation (45.5) is the 
well-known equation of which the solution is the associated 
Legendre polynomial Pftcosfl). If I is integral and 
| u |, then /'“(cos 0 ) is the only solution which is 
bounded in the range and is therefore the only 

one leading to a wave function ij/ which satisfies condition 
(ii) above; if /is not an integer no bounded solution exists. 
To solve equation (45.5) we write 


a 


2 


8 nhnW _ 4n 2 mZe 2 

h 1 ’ V h 2 a 


(45.7) 


and change the independent variable from r to x where 
x = 2ir. We then find that 

+ (45.8, 

ax x ax (4 x x J 


where, in order that conditions (iii) and (iv) should be 
satisfied, R must be such that J?-»0, as x-*ao and as 
x->0. From the arguments of § 44 it foUows that this 
is possible only if / is a positive integer and only if v is 
an integer, n say, which is greater than /+1. When this 
is so we may write the solution of equation (45.8) in the 
form (44.7) so that the solution of (45.5) is proportional 
to R„ t {2ar). Now by the second of the equations (45,7) 
we have 

2K 2 Z 2 e*m 

ftV 


w 


(45.9) 
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for the possible values of the total energy W. The wave 
functions corresponding to the value of energy given by 
the integer n are 

0, 4>) = C„, u fi nl (2ar)/>“(cos 0)e ,u + (45.10) 


where C„ iu is a constant determined by the condition (v). 
In polar coordinates dx = r 2 sin 0 drdO d<f> so that this 
condition gives 

1 = Cni» | d<j> f sin 0{/’“(cos 0)} 2 d0 f r 2 R 2 t (2xr)dr 
Jo Jo Jo 


so that it foUows from equations (21.20) and (44.9) that 


(2ai) i f (-lfti-«)l(ii-I-l)!(2I+l) ]» 

1 8* h (/ + h)!{(b + I)!} 3 J ' 


(45.11) 


Now from equations (45.7) and (45.9) we find that 
a = An 2 Ze 2 m}{h 2 n) so that if we introduce the Bohr radius 
a by means of the equation 


h* 

47i J fne 2 


(45.12) 


we find that a = Zjan. Introducing this result into (45.11) 
and substituting the value obtained for QtJa into equation 
(45.10) we find that 


'I'nluir, 0, 4>) 



(/-u)!(«-l-i)!(2f+l) j* 

(/+ti)!{(n+J)l} 3 8ffii j 


/?„ 1 (2ar)Pj'(cos 0)e ,u + 


(45.13) 


with w < /g « — 1 are the wave-functions corresponding to 
the energy (45.9) of the hydrogen atom. 

If we write = —2n 1 Z 2 e 4 mjh 2 then corresponding 
to the energy W 0 we have the wave function 


'l'ioo(r, 0, <t>)= e * 


(0 = Zrja), 
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and to the energy level W 0 j4 we have the three wave 
functions 

0.« - (f ) 1 e*-*" cos 0. 

feu {r, ft *) “ g^j(f J &e ~ ie sin de± ‘*- 

From the last of these functions we can construct two 
functions 

1 {Z\* . „ sin , 

—= [ — 1 Qe * e sin 0 <f>. 

Ayjln W cos 

Similarly corresponding to the energy W 0 j9 we have the 
six wave functions 

*300(1V e , *) = (27-18 e+ 2 e 2 ) e -* 

*310 (r, 0 , *) = JL (?Y (6 -e)Qe~*° cos 0 , 

81 V«W 

*3ii(' - . 0 , *) = Q (6-eJee - ** sin Qe ±l *, 

*320(r, fl, *) = gi^ (a) e**"*^ 3 cos * e-1 >> 
* 32 i(f, 0, *) = Q 1 e~ i<l sin 0 cos Ge*'*, 

*^ r - *♦>-1 ^{-j e ' e ~' ,s> °‘ ee± ‘ u - 
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in wave mechanics j *(r, 0, *) j 2 dz represents the 
probability that the electron whose wave-function is 
*(t, 0, *) is to be found in a small volume dx centred at 
the point whose polar coordinates are (r, 0, *). To make 
the total probability unity we must have 

J|*(r, 0, *)| 2 dr = 1, 

where the integral is taken throughout the whole space. 
Since, in polar coordinates, dx = r 2 sin 0 drd 0 d<j> it follows 
that the probability that the electron is at a distance 
between r—\dr and r 4- \dr from the nucleus is <b(r)dr 
where 


*(r) 


= r 2 pd* r 

Jo Jo 


sin OdO. | *(r, 0, 


Hence for an electron in the state defined by quantum 
numbers », /, u (equation (45.13) above) we have 




/2ZV (n—1—1)1 

\an) 2n{(n + /)!} 3 


r 2 {^ Bl (2ctr)} 2 . 


The mean value of any function / which depends on r 
alone is given by 





f(r)<f>(r)dr. 


that is, by the formula 


/= 


/ 2 ZV (»—/—!)! 

\ an) 2n{(n + /)!} 3 


r. 


r 2 f( r ){Jt n (2ar)} 2 dr. 


For examples of the use of this formula see example 17, 
p. 178. 
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Examples V 

(1) Prove that if m<n 


(2) If 


A m 2 m n 1 

— {//„(*)} = 
dx mX ' (n-m)! 


y, 0 = £ WJfW'nO')'". 

n ” 0 


where the tp n (x) are the orthonormal set of Hermite functions 
defined by the relation 

prove that 

e~ ixl+Ux ~ xl K(x, y, t)dx = 


i: 


Assuming that A'tx, y, t) is of the form 
A exp (Bx 2 + Cxy + Dy 2 ) 

where A, B, C and D arc functions of t prove that 


f 4xyt-(x 2 + /)(l + t 2 ) 


(3) Show that 


} 


B O 0 




(4) The Schrodinger equation for the three-dimensional 
oscillator in cylindrical coordinates (g, <j>, z ) is 

VV+ -p- {W—'ln 2 m{v 2 Q 1 +o> 2 z 2 )}\li = 0 . 


Show that the solutions ij/ of this equation such that 
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^-+0 as g and z~+oo and i j> is finite at the origin are of 
the form 

C exp[ iu-ittY-lfi 2 ! 2 }/,, „ 
where C is a constant, I, n, | u | are positive integers, 
* = W(wv/A), p = 2 ityJ(motlh), and is a poly¬ 

nomial of degree 21 in £ which satisfies the equation 


€L 
de 



Show also that the corresponding values of W are given 
by the equation 


W - (2/+| u | + l)Av+(/i + £)Aa>. 


(5) Prove that 

L„(2x) = nl 

(6) Prove that 


K 





Ln-mM. 


(*) fLyx) = L” +, ( x y, 

dx 


(ii) £-: +1 (x)+(x-2/>-l)L»+mL':- 1 (x)+B 2 i;_ 1 (x) = 0 

(iii) + = 0. 

(7) If wi£n prove that 

[ o JmV-Jx)e-* / ‘x B -* H dx = (-ir a” + i e-°L”(a). 


(8) Prove that if m£n, a> 0, 

f* ,-^iaxvi, - 

Jo ’ 

and deduce that if n^l+l 



x‘ +l RJx)dx~ 


{(n + Z). 1 } 2 / <y+fy 2 t +2 

(n-l-tV. 
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(9) Show that 

j 1 x”(\-xyn;(ax)dx = (-ir l !(<»)■ 

(10) Prove that the function 

x»e-* x LU,(x) 

is a solution of the differential equation 

<Py ct -y + 1 dy , f 1 , 2/?+«+1 _ y(2a-y) | q_ 

dx 2 x dx \ 4 2x 4x 2 J 

(11) The potential energy for the nuclear motion of a 
diatomic molecule is closely approximated by the Morse 
function 

p'(r) = De~ 2 "-2De~". 

Show that the spherically symmetrical solutions of the 
Schrodinger equation with this potential are 

4>(r) = ^ {exp(-be’“ r )}(2hr)‘-"-’ i L^:^" t 1 (26e- w ), 

(0gn£fc-i), 

where b = 2n(2mD)j(ah), C„ is a normalisation constant, 
and the corresponding values of the energy are 

<s±ft+Qg£|. 

(12) Prove that the normalisation constant C„ in the 
previous example has the value 



N „ = { r ( 2f >-«)} 2 


n 



T(2h—2n + r—1) 
H ‘ 


where 
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(13) Show that in parabolic coordinates i, tj, <f> defined 
by the equations 

x = (in)* cos <t>, y = (in) * sin <f>* 2 - K£-if) 
Schrodinger’s equation for a hydrogen-like atom of nuclear 
charge Ze takes the form 


5_ 

si 



S4> 2 


+ 


2te 2 hi 

ft 2 


{W(i+n)+2Ze 2 }<f, = 0, 


Show that this equation has solutions of the form 

■«*■(£) «*■(£) 

where a is the Bohr radius, n = k + l+a+l and W = 
-2 K 2 me*z 2 !{n 2 h z ). 

Determine the constant C so that ^ is normalised to 
unity, 

(14) In the theory of the rotation and vibration spectrum 
of a diatomic molecule there arises the problem of solving 
the Schrodinger equation with potential energy 


Show that the energy levels arc given by 
W = ~ 2 n 2 mZ 2 e 4 l(h 2 <7 2 ) 

where a - 11 +£+<</{£+(/+i) 2 } with /#* / integers and 
b — 8 n 2 niBjh 2 , and that the corresponding wave functions 
are 

c„ ( -Y"* e-"°'LU' (z) P %cos 6)e?*, 

\uaj \aa/ 

where a is the Bohr radius, sc = 2 N /{h+(/+^) 2 } and C„ a is a 
normalisation factor. 
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(15) Show that the constant C„ a occurring in the last 
example has the value 

f 2 ]* f (2/+l)«!(1 -u)1 )» 

(£rjr(n + a+l)j {47r{2n + a+1)(1 + m)!J ' 

(16) If 

/(/; n; s) = J x ,+ 1 {/i Bl 0e)} 2 <fx, 

prove that /(/; n\ s) -i- {(/i + /)!} 2 is the coefficient of (rr) n + 1 
in the expansion of 

(I-l) i+ i(i- T y+> y (2Hs + r+2)! (fT)2l+f+l 
r r! 


Deduce that 

/(/; it; 


D=iM, 

C»—/—1)1 


/(/; n; 1) 


= {(« + Q!} 3 

(„_/_!)! 


{6n 2 —2/(/+1)}. 


(17) Show that in a hydrogen-like atom of nuclear 
charge Ze the average distance of the electron from the 
nucleus, in the state described by quantum numbers /, 
n, is 



Find the average value of I jr and show that the total 
energy of a hydrogen atom is just one-half of the average 
potential energy. 


APPENDIX 


THE DIRAC DELTA FUNCTION 


§ 46. The Dirac delta function. In mathematical 
physics we often encounter functions which have non-zero 
values in very short intervals. For example, an impulsive 
force is envisaged as acting for only a very short interval 
of time. The Dirac delta function, which is used extensively 
in quantum mechanics and classical applied mathematics, 
may be thought of as a generalisation of this concept. 

If we consider the function 


5 . 00 = 2 a 


0, 


x\>a. 


then it is readily shown that 

S a (x)dx = 1 . 


r 


(46.1) 


(46.2) 


Also, if J{x ) is any function which is integrable in the 
interval {—a, a) then, by using the mean value theorem 
of the integral calculus, we see that 

j* f(x)S/pc)dx = ~ [" f(x)dx = f{0a\ j 6 ) £ I. (46.3) 


We now define 


<S(x) = lim 5 a (x). 


(46.4) 


Letting a tend to zero in equations (46.1) and (46.2) we 
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see that £(x) satisfies the relations 

5(x) = 0, if x # 0, 



5(x)dx — 1. 


(46.5) 

(46.6) 


The “ function ” <5(x), defined by equations (46.5) and 
(46.6), is known as the Dirac delta function. It is unlike 
the functions we normally encounter in mathematics; 
the latter are defined to have a definite value (or values) 
at each point of a certain domain. For this reason Dirac 
has called the delta function an *' improper function ” and 
has emphasised that it may be used in mathematical 
analysis oniy when no inconsistency can possibly arise 
from its use. The delta function could be dispensed with 
entirely by using a limiting procedure involving ordinary 
functions of the kind £„(*), but the “ function ” <5(x) and 
its “ derivatives ” play such a useful role in the formulation 
and solution of boundary value problems in classical 
mathematical physics as well as in quantum mechanics 
that it is important to derive the formal properties of the 
Dirac delta function. It should be emphasised, however, 
that these properties are purely formal. 

First of all it should be observed that the precise variation 
of <5(x) in the neighbourhood of the origin is not important 
provided that its oscillations, if it lias any, are not too 
violent. For instance, the function 


JOO-lfan 

u-*oo 71X 

satisfies the equations (46.5) and (46.6) and has the same 
formal properties as the function defined by equation 
(46.4). 

If we let a lend to zero in equation (46.3) we obtain 
the relation 



f(x)S(x)dx =/{0), 


(46.7) 
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which a simple change of variable transforms to 



f(x)5(x-a)dx=f(a). 


(46.8) 


In other words the operation of multiplying f(x) by <5(jr-«) 
and integrating over all x is merely equivalent to sub¬ 
stituting a for x in the original function. Symbolically 
we may write 

f{x)5(x-a)=f(a)6(x~a) (46.9) 


if we remember that this equation has a meaning only in 
the sense that its two sides give equivalent results when 
used as factors in an integrand. As a special case we have 

= 0. (46.10) 

In a similar way we can prove the relations 

5(-je) = a(jc), (46.11) 

S(ax) = i 5(x), a>0, (46.12) 

a 


<5(a*-x 2 ) 


2a 


{3(jc-a) + <5(x + a)}, 


n>0. 


(46.13) 


Let us now consider the interpretation we must put 
upon the “ derivatives ” of <5(.v). If we assume that <5'(x) 
exists and that both it and <5(x) can be regarded as ordinary 
functions in the rule for integrating by parts we see that 

f" /(x)5'(x)dx = r/(x)5(x)T 

J - 00 L _J— DO 

-I: 

Repeating this process we find that 

f(x)5 M (x)dx = (-lT/^m. 


i: 


nx)S{x)dx = -/'(0). 


(46.14) 
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The statement is often made that the Dirac delta 
function is the derivative of the Heaviside unit function 
//(x) defined by the equations 





if x>0; 
if x<0; 


and it is easy to see on geometrical grounds that there are 
reasons for conjecturing such a relationship. To make it 
precise we note that if, in the definition | of the Stieltjes 
integral 



f(x)dF(x), 


we lake F(x ) to be the Heaviside function H(x), we find 
immediately that, for any integrable function f(x). 



f(x)dH(x)=f( 0). 


(46,15) 


Comparing equation (46.15) with equation (46*7) we see 
the relation between H(x) and 8(x), It may be seen from 
these equations that <5(x) is not a function but a Stieltjes 
measure, and that the use of the Dirac delta function could 
be avoided entirely by a systematic use of Stieltjes in¬ 
tegration. 

t The simplest definition of a Stieltjes integral J * f(x)dF(x) is as 

the limit of approximative sums Ef(£ r )[F(x T )— where the 
x ri the points of subdivision of (a, b) and £ rt lies in the interval (x ri jr f ). 
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